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ABSTRACT 


The  NPS  computer  simulation  model  was  modified  to  study 
the  angular  dispersion  and  ring  of  sputtering  sites 
formation  of  sputtered  copper  isotopes.  The  sputtering  ratio 
as  a  function  of  the  ion  mass  and  the  microcrystallite  size 
needed  to  contain  low-energy  events  were  also   investigated. 

A  shift  in  the  ejection  angles  of  sputtered  Cu  (63) 
isotopes  toward  the  surface  normal  is  observed.  When  the 
"primary  knock-on  atom"  was  a  light  isotope,  all  of  the 
sputtered  atoms  were  ejected  more  nearly  normal  to  the 
surface.  The  predominance  of  sputtered  atoms  in  the  <100> 
and  < 1 0 1  >  directions  from  the  impact  atom  indicated  an  ease 
of  momentum  transfer  in  these  directions.  Results  obtained 
from  studies  based  upon  a  grid  of  impact  points  compares 
very  favorably  with  the  results  obtained  when  randomly 
selected  impact  points  were  used.  The  rings  of  sputtered 
atoms  present  a  "liquid  splash"  effect. 
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I.   INTRODUCTION 


The  phenomena  of  physical  sputtering  have  been 
extensively  studied  since  it  was  first  mentioned  in 
Philosophical  Transactions  by  Grove  in  1852.  Sputtering  is 
the  ejection  of  atomic  material  from  a  target  resulting  from 
bombardment  by  an  energetic  particle.  Most  scientific 
investigation  has  been  restricted  to  ionic  bombardment 
because  the  ion  can  easily  be  accelerated  to  any  desired 
velocity.  Sputtering  has  defied  complete  physical 
description  because  of  complicated  mechanisms,  difficulty  in 
treating  the  behavior  of  a  many-body  system,  and  an 
incomplete  understanding  of  the  interacting  forces  and 
potentials  in  crystal  structures. 

3arly  experiments  typically  determined   the   "sputtering 
ratio"    of    single   and   polycrystalline   materials.    The 
"sputtering  ratio"  (interchangeable  with  sputtering  yield 
Y)  is  defined  as  the  number  of  target  atoms  ejected  from  the 
surface  per  incident  ion. 

In  1923,  Kingdon  and  Langmuir  [1]  bombarded  thoriated 
tungsten  with  various  ions  in  a  glow  discharge  tube.  This 
was  a  special  case  of  sputtering  since  the  thin  surface  film 
of  thorium  on  a  tungsten  substrate  was  sputtered  rather  than 
the  tungsten  itself.  The  results  of  this  experiment 
gualitatively  suggested  an  ejection  mechanism. 

A  few  years  later,  Blechschmidt  and  von  Hippie  [2,3] 
proposed  a  theory  which  described  sputtering  as  an 
evaporation  of  surface  atoms  (an  emission  mechanism). 
Earlier,  von  Hippie  [U]  had  found  by  spectroscopic:   analysis 
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that  at  least  some  sputtered  atoms  were  in  an  excited  state. 
The  sputtering  theory  showed  that  atoms  in  the  region  of 
impact  could  rapidly  acquire  thermal  energy  if  the  kinetic 
energy  of  the  incident  ion  was  converted  to  thermal  energy 
at  the  target  surface.  If  an  atom  acquired  a  sufficient 
amount  of  thermal  energy,  it  would  then  evaporate  from  the 
surface,  some  atoms  evaporating  while  in  an  excited  state. 

Direct  application  of  statistical  methods  to  sputtering 
was  made  by  Harrison  [5]  who  envisioned  the  interaction  of 
two  distribution  functions  (the  crystal  lattice  and  the  ion 
beam) .  These  models,  based  on  statistical  methods, 
implicitly  accept  ejection-type  mechanisms. 

One  of  the  most  important  contributions  to  the  study  of 
sputtering  was  made  by  Wehner  [6]  in  1953.  He  discounted 
the  evaporation  theory  and  presented  strong  evidence  for  a 
momentum  transfer  process.  Shortly  after  Wehner's  findings 
were  reported,  Henschke  [7]  proposed  a  theory  of  sputtering 
based  solely  on  classical  theory,  treating  normal  and 
oblique  incidence  sputtering  separately.  His  concept  of 
normal  incidence  sputtering  required  many-body  collisions  in 
which  the  ion  was  eventually  reflected  outward  to  sputter 
surface  atoms.  This  theory  was  plausible  for  oblique 
incidence  sputtering;  however,  the  case  for  normal  incidence 
sputtering  required  that  the  ion  be  reflected  inside  the 
crystal. 

Wehner  and  Rosenberg  [8]  reported  that  in  the  sputtering 
of  polycrystalline  metals  at  low  ion  energy  (normal 
incidence)  more  material  is  ejected  obliquely  than  in  a 
direction  normal  to  the  target  surface.  This  is  due  to  the 
fact  that  less  than  a  full  (180°)  overall  direction  change 
of  momentum  is  more  favorable  for  sputtering  because  this 
requires  fewer  number  of  successive  collisions  for  atom 
ejection.     In   1960,   Wehner   and   Anderson   [9]   reported 
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sputtering  patterns  from  single  crystal  targets  which  showed 
an  anisotropic  spread  of  energy  from  a  collision  center  and 
atoms  which  are  preferentially  ejected  in  close- packed 
directions. 

When  an  atom  collides  with  another  which  has  a  nigher 
mass,  the  scattering  angla  can  vary  between  0°  and  180°.  In 
a  collision  with  a  lower  mass  atom  the  scattering  is 
confined  to  between  0°  and  90°  (thus  backward  reflection 
cannot  occur) .  In  low-energy  sputtering  of  multi-element 
targets,  Winters  and  Sigmund  [10]  showed  theoretically  that 
a  significant  contribution  to  sputtering,  called  "single 
reflective  collisions",  arises  from  the  reflection  of  lower 
mass  atoms  from  underlying  heavier  atoms.  For  a  detailed 
discussion  of  momentum  reversal,  see  Harrison  and  Magnuson 
[11  ]. 

In  1962,  Wehner  et  al  [12]  published  a  compilation  of 
sputtering  yields  for  metals  and  semiconductors  in  the 
100-600  eV  range.  In  1966,  Wehner,  Anderson,  and  KenKnight 
[13]  published  the  results  of  a  study  that  covered  almost  a 
decade  of  research  in  low-energy  sputtering.  A  few  of  the 
numerous  topics  discussed  included  more  accurate  and 
improved  techniques  for  analyzing  the  sputtered  mass,  the 
temperature  dependence  of  sputtering  yields  of  variolous 
metals  bombarded  by  atomic  and  polyatomic  ions  in  the  2.8  to 
10  keV  range,  the  effect  of  oxygen  on  ion  ejection  patterns, 
and  The  spatial  distribution  of  ejected  atoms  in  metals  and 
semiconductors. 

A  few  low-energy  sputtering  experiments  [14,  15,  16,  17, 
13],  which  lend  themselves  to  direct  comparison  with  these 
simulations,  investigate  the  behavior  of  the  sputtering 
ratio  as  a  function  of  ion  energy,  lattice  orientation, 
target  material,  and  ion  angle  of  incidence. 
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Paralleling  the  experimental  investigation  of  the 
sputtering  yield  and  the  development  of  mass  spectrographic 
analysis  was  the  attempt  to  model  the  sputtering  mechanism 
through  the  use  of  computer  simulation.  In  1960,  Gibson, 
Goland,  Milgram,  and  Vineyard  [19],  the  first  investigators 
to  use  computer  simulation  as  a  means  of  studying  radiation 
damage,  set  forth  the  criteria  that  must  be  satisfied  before 
the  simulation  method  could  be  applied  to-  crystals. 
Reference  [19]  gives  a  very  detailed  description  of  the 
simulation  method.  Basically,  they  built  a  computer  model 
to  represent  metallic  copper  and  studied  radiation  damage 
events  at  low  and  moderate  ranges  of  energy  up  to  400  eV. 
In  their  model,  one  atom  in  the  stationary  lattice  was  given 
an  arbitrary  kinetic  energy  and  direction  of  motion,  as 
though  it  had  been  struck  by  a  bombarding  particle,  and  was 
allowed  to  interact  with  the  remaining  atoms  in  the  lattice. 
This  resulted  in  cascades  of  independent  binary  collisions. 
It  is  well  established  that  sputtering  of  solids  by  ion 
bombardment  is  the  result  of  successive  independent 
approximately  binary  collisions  somewhat  as  if  these 
collisions  take  place  with  or  between  separate  gas  atoms. 
Low-energy  ion  scattering  studies  [20]  have  shown  that  the 
approximately  binary  scattering  behavior  still  holds  at 
bombarding  energies  as  low  as  those  approaching  the 
sputtering  thresholds. 

In  1967,  Harrison,  Levy,  Johnson,  and  Effron  [21]  used  a 
computer  to  simulate  the  bombardment  of  f ace-centered-cubic 
(fee)  copper  single  crystals  by  an  argon  ion.  From  their 
studies,  the  sputtering  mechanisms  in  the  fee  crystal  were 
isolated  and  identified.  They  reported,  for  ion  energies 
less  than  10  keV,  the  sputtering  process  occurred 
predominately  within  three  atomic  layers  of  the  surface. 
The  computer  simulation  included  only  a  repulsive  force 
between    the   atoms.    The   equations   of   motion   for   the 
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many-body  system  were  computed  using  an  "average- force" 
algorithm  rather  than  a  "central-difference"  algorithm,  see 
Harrison,  Gay,  and  Effron  [22].  In  1972,  Harrison,  loore, 
and  Holcombe  [23]  in  continuing  efforts  to  develop  a  more 
realistic  computer  .model  incorporated  a  Cu-Cu  potential 
function  which  included  an  attractive  potential,  surface 
layer  relaxation,  and  surface  binding  energy  for  the  (100), 
(110)  ,  and  (111)  copper  planes.  They  reported  improved 
agreement  between  experimental  and  simulation  results,  but 
the  original  interpretation  of  the  low-energy  sputtering 
mechanisms  in  fee  copper  crystals  remained  unchanged. 

Computer  simulation  of  physical  processes  is  still  in  a 
very  infant  state  because  the  available  processing  speed  of 
existing  computers  severely  limits  the  sophistication  of  the 
models.  However,  past  investigations  were  quite  successful 
in  identifying  the  mechanisms  which  occur  when  energetic 
particles  interact  with  crystal  atoms.  The  ability  to 
observe  these  mechanisms  is  an  advantage  peculiar  to  the 
computer  simulation  because  each  crystal  atom  is  identified 
by  a  number  for  the  mathematical  calculation,  and  its 
complete  track  can  be  labeled  and  recorded.  The  tracks  of 
selected  atoms  can  be  examined  to  determine  the  dynamics  of 
the  atoms  which  exit  the  front  surface. 
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II.   RESEARCH  OBJECTIVES 


A.   ANGULAR  DISPERSION  OF  ISOTOPES 


Recently,  Wehner  [24]  reported  that  the  isotope  ratio  in 
metal  deposits  obtained  by  sputtering  a  flat  target  at  low 
ion  energy  (100e7)  under  normal  incidence  is  a  ^unction  of 
the  ejection  angle.  He  investigated  the  bombardment  of  a 
Cu  (63-65)  target  with  100  eV  Hg  ions  .  He  reported  that  the 
lighter  isotopes  are  enriched  in  a  direction  normal  to  the 
target  surface.  This  research  was  a  computer  simulation 
investigation  of  Wehner1 s  studies  of  angular  dispersion  of 
isotopes  in  a  monocrystalline  copper  lattice. 


B.   LIQUID  SPLASH  EFFECT 


In  preliminary  NPS  simulation  studies,  after  analysis  of 
data,  a  predominance  of  sputtered  atoms  occurred  at  certain 
distances  from  the  center  of  the  lattice  representative  area 
in  which  the  ion  impacts  occur.  Figures  1  and  2  [  25  ]  depict 
the  top  layer  of  a  copper  fee  crystal  respectively  bombarded 
by  neon  and  gold  ions  at  2  keV.  The  numbers  represent  the 
number  of  times  that  the  atom  was  sputtered  in  36  events 
(shots) .  The  dark  triangle  in  the  center  represents  the 
target  impact  area.  All  36  events  are  directed  into  this 
representative  area.  The  darkened  elipses  highlight  the 
atoms  most  often  sputtered  and  suggest  a  ring  outward  from 
the  center  (especially  for  the  Cu-Au  interaction) . 


21 


Macklin  and  Metaxas  [26]  experimentally  investigated  the 
splashing  of  liquid  drops  on  deep  and  shallow  liquid  layers. 
They  published  high  speed  photographs   which   showed   liquid 
splashes   that   might   be   analogous   to  surface  atoms  baing 
ejected  from  the  surface  and  a  few  eventually  sputtering. 

This  research  was  an  investigation  of  "ring  formation" 
or  "liquid  splash"  effect  of  sputtered  atoms.  Determination 
was  made  of  the  probability  that  a  congregate  of  atoms  at  a 
given  distance  from  the  impact  area  reference  will  sputter. 


C.   ENERGY  CONTAINMENT 


As  with  any  computer  simulation,  a  trade-off  must  be 
made  between  the  accuracy  and  detail  desired  and  the  size 
and  time  requirements  imposed  by  computer  availability. 
Because  of  this  constraint  a  limit  exists  on  the  size  of  the 
monocrystalline  lattice  used  in  the  simulation.  This 
research  was  an  investigation  to  determine  the  size  lattice 
required  to  contain  a  low  energy  event. 


D.   ION/ATOM  MASS  RATIO 


In  July  of  1976,  Eer  Nisse  [27],  reported  experimental 
confirmation  of  a  monotonic  dependence  of  Y  with  the  atomic 
number.  He  investigated  the  yield  of  a  gold  target  impacted 
by  40  keV  ions  up  to  the  atomic  number  of  gold.  This 
monotonic  dependence  is  predicted  by  a  sputtering  theory  of 
Sigmund  [28].  However,  computer  simulation  data  of  previous 
NPS  studies  [29]  indicate  that  the  monotonic   dependence   of 
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the  sputtering  ratio  on  the  ion  mass  ceases  to  be  valid  for 
ion/atom  mass  ratios  greater  than  one.  This  research  was  an 
investigation  of  the  sputtering  ratio  versus  ion  atomic 
number  for  a  wide  range  of  ion/atom  ratios. 
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III.   SIMULATION  MODEL  MODIFICATIONS 


The  basic  simulation  model  used  in  this  research  has 
been  developed  by  the  NPS  group  over  the  last  decade  [21, 
23,  31,  32].  A  brief  description  of  this  model  and  the 
alterations  made  to  it  for  this  research  follow.  For 
detailed  discussions  of  the  simulation  method,  the  reader  is 
directed  to  references  19,  20,  21,  and  23.  A  detailed 
description  of  the  mathematical  model  has  been  published 
[22]. 

In  the  development  of  the  computer  model,  the  lattice 
unit  (LU)  proves  to  be  a  convenient  unit  of  length,  and  is 
used  extensively.  One  LU=a0/2,  where  a„  is  the  cubic 
lattice  constant,  the  atomic  spacing  as  determined  by  x-ray 
crystallographic  studies.   For  copper,  a0 =3.615  angstroms. 

Certain  properties  of  the  program  and  results  are  best 
considered  in  the  framework  of  a  specific  run.  Each  run  is 
begun  by  constructing  a  target  microcrystallite  whose  sites 
represent  the  equilibrium  positions  of  copper  atoms  near  the 
surface  of  a  fee  crystal.  The  target  atom  -  target  atom 
interaction  is  represented  by  the  Born-Mayer  Gibson  type  II 
(Gibson-II)  repulsive  potential  function. 

Each  ion  of  the  beam  is  represented  by  a  single  neutral 
atom  of  the  desired  element  [30].  The  ion  is  placed  at  the 
point  on  its  trajectory  sach  that  its  velocity  vector 
determines  the  point  of  impact.  The  term  bullet  and  ion  are 
used  interchangeably  throughout  this  report  for  the  incident 
particle  to  avoid  confusing  it  with  target  atoms. 
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The  bullet-target  potential  function  was  always  either 
the  Gibson-II  or  that  determined  by  the  Kinetic  Secondary 
Electron  (KSE-B)  Potential  [3°]-  These  two  potential 
functions  seem  to  represent  a  reasonable  range  of  "hard"  and 
"soft"  potential  functions. 

As  reported  by  Harrison,  ?1oore,  and  Holcombe  [23],  the 
attractive  forces  between  the  atoms  in  the  crystal  are  so 
small,  when  compared  to  the  forces  which  occur  in  even 
moderately  energetic  collisions,  that  they  may  be  safely 
neglected  in  sputtering  simulations.  Also,  in  order  to  use 
the  smaller  reduced  target  areas,  the  ions  must  be  normally 
incident.  Even  though  the  computer  program  still  retained 
the  capability  to  utilize  attractive  potentials  and 
obliguely  incident  ions,  they  were  not  used  in  this 
investigation . 


A.   MICEOCRYSTALLITE  SIZE 


Initial  discussion  to  establish  research  objectives 
revealed  the  advantages  gained  from  using  a  symmetrical 
lattice  that  could  be  folded  and  rotated  so  -chat  atoms 
sputtered  by  ions  incident  at  all  possible  impact  points 
would  become  available  for  analysis.  A  lattice  with  each 
layer  a  square  and  having  an  odd  number  of  atoms  along  each 
boundary  was  needed.  The  lattice  generation  program  as 
previously  used  could  handle  this  change  with  no 
alterations. 

Figures   3,  4,  and  5  show  a  roughly  isometric  projection 

of  the  labeled  atoms  in  a  9x4x9  (100)  orientation,   a   9x5x9 

(110)   orientation,   and   a  9x4x9  (111)  orientation  of  a  fee 

microcrystallite  respectively.    This   type   of   projection, 

first  used  by  Harrison  and  Delaplain  [31],  with  the  distance 
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between  layers  in  the  Y  direction  exaggerated,  and  the  atoms 
represented  by  20°  ellipses  has  proved  to  be  a  valuable  tool 
for  analyzing  the  computer  programs  and  for  displaying  data. 
The  first  number  in  each  lattice  is  atom  number  two,  with 
atom  number  one  being  the  bullet. 

From  previous  studies,  it  was  clear  that  at  the  low 
energy  levels  considered  the  sputtering  mechanisms  were 
confined  to  the  surface  layers  [21].  For  this  reason,  the 
(100)  and  (111)  orientation  microcrystallites  were  only  four 
layers  deep.  The  (110)  orientation  lattice  generation 
program  constructs  a  lattice  which  has  significantly  fewer 
atoms  and  its  layers  are  closer  together.  In  an  attempt  to 
increase  the  number  of  atoms  in  the  lattice  while  keeping 
the  required  computer  time  below  that  of  the  (100)  and  (111) 
orientation  programs,  the  (110)  orientation  microcr ystallite 
was  set  at  six  layers  deep. 


B.   CRYSTAL  EOONDARIES 


The  previous  simulation  program  contained  a 
microcrystallite  with  established  boundaries  that  were  not 
consistent  with  the  development  of  the  potential  functions. 
Since  the  potential  function  is  truncated  at  the 
nearest-neighbor  separation  (ROE) ,  all  boundaries  (except 
the  top  surface)  were  set  at  ROE.  This  is  justified  by  the 
fact  that  an  atom  outside  these  boundaries  is  farther  from 
the  crystal  than  the  potential  truncation  radius.  The  upper 
surface  boundary  was  set  by  starting  a  small  distance  from 
the  lattice  surface,  and  increasing  the  distance  in 
succeeding  runs  until  the  total  yield  was  independent  of  the 
distance.  For  computational  efficiency  this  distance  should 
be  set  as  small  as  possible. 
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Figure  6  shows  a  drawing  of  the  microcrystallite  in 
solid  lines  and  the  plane  boundaries  in  dashed  lines.  The 
boundaries  which  are  defined  in  the  Glossary  (Appendix  C) 
are  shown  as  planes  at  the  distance  from  the 
microcrystallite  surface  cited  above. 


C.   ISOTOPE  GENERATION 


In  order  to  study  the  angular  dispersion  of  sputtered 
isotopes,  a  method  for  generating  a  microcrystallite  with 
the  correct  isotope  ratio  and  masses  was  needed.  The  method 
used  assigned  an  additional  variable  to  each  atom  that 
identified  the  atom  "type".  The  atoms  were  classified 
according  tc  two  types;  type  1  -  least  abundant,  heavier 
isotope  and  type  2  -  most  abundant,  lighter  isotope. 
Appendix  A  is  a  detailed  discussion  of  how  the  correct 
isotope  ratio  was  obtained  and  each  atom  properly  identified 
with  the  correct  mass. 


D.   TARGET  IMPACT  AREAS 


A  representative  area  must  contain  all  possible  impact 
points  upon  a  particular  surface,  that  is  any  point  in  the 
surface  must  be  projectable  into  a  point  of  the 
representative  area.  Figure  7  shows  the  intrinsic  geometry 
of  the  three  fee  orientations  and  the  complete  sets  of 
impact  points  for  each  representative  area.  Figure  8  shows 
the  reduced  sets  of  impact  points  from  which  the  complete 
sets  are  generated  by  reflection  and  rotation.  Figures  9, 
10,  and  11  show  an  expanded  view  of  the  target  impact  areas 
and  the  actual  impact  points  used.  The  individual  impact 
points  are  those  used  in   previous   simulations   for   normal 
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incidence.  This  simplification,  which  greatly  reduces  the 
computer  running  time,  is  only  possible  when  the  ion  is 
normally  incident. 


E.   LATTICE  FOLDING 


In  order  to  investigate  the  ring  formation  or  "liquid 
splash"  effect,  it  was  necessary  to  develop  a  means  of 
folding  the  results  obtained  from  the  reduced  target  area. 
This  requirement  was  satisfied  by  writing  three  subroutines; 
KFOLDA,  KFOLDB,  and  KFOLDC.  Appendix  B  gives  a  detailed 
description  of  the  "folding"  logic  and  the  geometry  of  each 
orientation.  Figure  12  shows  how  the  (100)  and  (110) 
orientation  reduced  representative  areas  were  folded  to 
obtain  the  complete  representative  areas.  Figure  13  shows 
how  the  (111)  orientation  reduced  representative  area  was 
folded  and  summed  to  obtain  the  complete  representative 
area. 


F.   ATOM  SPUTTERING  PROBABILITY 


A  means  of  taking  the  folded  data  and  calculating  a 
conditional  probabilty  that  an  atom  would  sputter  from  a 
particular  ring  was  needed.  This  was  accomplished  by 
writing  the  subroutine  KPROB.  In  this  subroutine,  the 
distance  from  each  atom  to  the  reference  (0,0)  point  in  tne 
reduced  representative  impact  area  was  calculated.  then  all 
atoms  at  each  distance  were  checked  and  the  atom  sputtering 
count  was  summed.  This  sum  was  divided  by  the  total  number 
of   sputs   from   all   the   shots.    This   quotient   was   the 
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conditional  probabilty  that  a  given  ring  at  a  known  distance 
would  sputter  an  atom  given  the  atoms  that  had  already 
sputtered . 

KPR03  performs  the  above  calculations.  However,  a  very 
interesting  result  is  obtained  by  dividing  the  conditional 
probabilty  obtained  above  by  the  number  of  atoms  in  a 
particular  ring.  This  result  was  a  conditional  probability 
that  an  atom  at  a  given  distance  from  the  impact  point  would 
sputter. 


G.   EJECTION  ANGLES 


The  ejection  angle  is  defined  as  the  angle  between  the 
atom's  velocity  vector  and  the  surface  normal.  A  tally  was 
kept  of  each  atom's  ejection  angle,  grouped  in  five  degree 
intervals.  This  tally  was  recorded  separately  for  each 
isotope.  To  provide  a  guick  analysis  of  the  data,  a  program 
(HISTO)  was  written  to  print  out  a  histogram  of  the  number 
of  sputtered  atoms  ejected  in  each  interval  for  each 
isotope. 


H.   VAHIABLE-SHOT  SELECTION 


The  previous  program  was  only  flexible  enough  to  allow 
either  a  single-shot  run  or  all  36  shots.  As  the  size  of 
the  lattice  and  the  ion's  energy  were  increased,  the 
computer  run  time  for  a  set  of  36  shots  became  unmanageable 
(12  -  18  hours) .  At  this  point,  a  variable  size  shot 
"package"  was  developed.  The  features  included  the  ability 
to  specify  the  number  of  shots  desired  and  the  starting 
impact   point   of   the  "package".   The  final  summary  for  the 
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particular  package  run  was  printed  out.  The  packages  had  to 
be  combined  by  hand  calculations  in  order  to  get  the  total 
information  for  a  particular  run.  Even  so,  this  procedure 
allowed  runs  to  be  made  with  larger  lattices  and  higher 
energies  (i.e.  the  containment  studies)  that  previously 
would  have  been  prohibitively  long. 


I.   RANDOMLY  SELECTED  IMPACT  POINTS 


To  obtain  more  angular  dispersion  data  for  the 
low-energy,  (100)  orientation  events,  a  modification  was 
made  to  the  simulation  model.  This  modification  replaced 
the  systematic  method  of  incrementing  the  impact  point  over 
a  grid  of  points  in  the  representative  target  area  with  a 
random-shot  capability.  Again,  the  IMSL  GGU3  routine  was 
used  to  generate  uniformly  distributed,  psuedo-random 
numbers.  The  first  random  number  was  labeled  the  x 
coordinate  in  the  impact  area.  The  second  random  number  was 
labeled  the  z  coordinate  and  tested  to  determine  if  it  was 
inside  the  impact  area  z  boundaries  determined  by  the  value 
of  the  x  coordinate.  If  this  test  was  passed  sucessfully, 
these  values  for  the  x  and  z  coordinates  were  used. 
However,  if  the  test  did  not  pass  (the  x  and  z  coordinates 
were  outside  the  impact  area),  the  two  original  numbers  were 
discarded,  two  more  were  selected,  and  the  process  was 
repeated  again. 

This   random   impact  point  did  not  effect  the  variable  - 

shot  selection  feature.  But,  enough  random  numbers   had   to 

be   generated   to   allow  for  those  points  that  did  not  fall 
into  the  impact  area. 
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The  results  of  randomly  bombarding  the  impact  area  will 
be  discussed  later.  This  feature  was  incorporated  in 
response  to  other  investigators  suggestion  that  the  grid  of 
impact  points  be  compared  with  a  randomly  bombarded  impact 
area. 
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IV.   RESULTS 


The  three  lattice  orientation  models  for  the  copper 
target  were  utilized  in  computer  simulation  runs  in  which 
100-eV,  200-eV,  500-eVf  and  1-keV  normally  incident  ions 
bombarded  the  target.  The  ion  masses  varied  from  '4  amu 
(He+)  to  over  200  amu  (Hg+)  .  The  lattice  was  constructed 
with  copper  atoms  of  either  the  natural  mass  (63.54  amu)  or 
the  proper  ratio  of  isotopes  (63.93  amu  and  64.93  amu)  .  The 
lattice  size  varied  from  a  9x4x9,  (100)  orientation  of  162 
atoms  to  a  17x4x17,  (100)  orientation  of  578  atoms. 

Computer  runs  were  made  for  various  permutations  of  ion 
energy,  ion  mass,  lattice  orientation,  lattice  size,  and 
potential  function.  The  summary  of  each  run  contained  the 
number  of  atoms  sputtered  from  each  impact  point,  the  number 
of  times  each  atom  sputtered,  a  histogram  showing  the 
angular  dispersion  of  the  sputtered  atoms  for  each  isotope, 
the  number  of  spurs  for  each  atom  after  the  folding  process, 
the  distance  from  each  ring  to  the  impact  area  reference 
point,  the  number  of  atoms  and  the  sum  of  their  sputs  in 
each  ring,  and  the  probability  that  each  ring  contained  a 
sputtered  atom. 

In  the  sputtering  summaries  the  numbers  in  the  ellipses 
indicate  the  number  of  times  that  particular  atom  spattered 
in  the  simulation.  The  dark  areas,  a  triangle  for  the  (100) 
and  (111)  orientations  and  a  rectangle  for  the  (110) 
orientation,  indicate  the  location  of  the  target  impact 
area.  The  atom  that  is  struck  by  the  ion  is  called  the 
"primary  knock-on  atom"  (hereafter  PKA)  .  Since  different 
seeds   were   used   to   generate   the  random  isotope  sites,  a 
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letter  L  or  H  denotes  whether  the  sputtered  atom  was  a 
light  or  heavy  isotope.  Wherever  possible,  the  results  for 
the  Gibson-II  and  KSE-B  potential  functions  were  combined. 
In  the  figures  displaying  the  sputtering  summary  and  the 
mechanism  tracing,  some  of  the  lower  layers  were  omitted. 
However,  all  atoms  were  included  in  all  calculations.  In 
the  100-eV  results  all  the  sputtered  atoms  are  shown. 
Wherever  a  direction  in  the  lattice  is  specified,  it  is 
always  with  respect  to  that  particular  lattice  axes.  The 
top  surface  is  always  an  XZ  plane. 


A.   9X9  3ICR0CRYSTALLITE 


Initially,  all  simulations  used  a  9x9  lattice  to 
determine  the  sputtering  ratio  as  a  function  of  the  ion  mass 
at  100  eV  and  1  keV.  The  model  was  modified  to  include  the 
isotope  generator,  the  symmetrical  folding  of  the  summary, 
the  angular  dispersion  of  sputtered  isotopes,  and  the  "ring 
sputtering"  probabilities. 

1 •   110 0L  Orientation 

The  (100)  orientation  microcrystallite  was 
constructed  in  a  lattice  consisting  of  162  atoms.  The 
lattice  was  9  atoms  in  the  X  and  Z  directions  and  4  atoms 
deep  in  the  Y  direction,  hence  9x4x9.  Figure  3  shows  an 
isometric  projection  of  the  lattice  and  its  atom  numbering 
scheme.  Figure  9  shows  an  expanded  view  of  the  grid  of  36 
impact  points  used. 
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a.   100-eV  Ions 

Of  primary  importance  was  the  investigation  of 
the  angular  dispersion  of  isotopes  reported  by  Wehner  [24  1. 
His  experiments  were  conducted  with  100-eV  mercury  ions 
bombarding  a  copper  target.  These  conditions  were  ideal  to 
simulate  with  modifications  to  the  existing  model. 

The  sputtering  simulations  for  a  normally 
incident,  100-eV  ion  were  run  using  He+,  Ne+,  Ar+,  Cu+, 
Kr+,  Xe  +  ,  Au+,  and  Hg+  ions.  Figure  14  shows  the 
sputtering  ratio  as  a  function  of  the  ion  mass  for  both 
potential  functions.  As  mentioned  in  the  objectives,  a 
decrease  in  the  yield  was  evident  for  the  heavier  ions.  The 
KSE-B  curve  shows  a  more  natural,  smooth  behavior  than  that 
of  Gibson-II. 

(1)  Helium  Ions.  The  results  when  100-eV 
helium  icns  impact  the  target  area  are  shown  in  figures  15, 
16,  and  17.   Figure  15  shows  the  sputtering  summaries   using 

(a)  Gibson-II,  resulting  in  a  sputtering  ratio  of  0.444,  and 

(b)  KSE-B,  resulting  in  a  sputtering  ratio  of  0.111.  The 
predominance  of  the  sputtered  atoms  in  the  <100>  and  <  1 0 1 > 
directions  is  quite  evident.  Figure  16  shows  the 
conditional  probability  that  a  ring  of  atoms,  as  a  function 
of  the  ring  distance  from  the  impact  area  reference  point, 
will  contain  an  atom  that  sputtered  for  (a)  Gibson-II  and 
(b)   KSE-B.    The   term   "ring    sputtering"    probability, 

henceforth  ring  probability,  was  used  as  a  label  -co  mean  the 
above  conditional  probability.  The  relative  amplitudes  of 
the  ring  probabilties  could  be  misleading  because  the 
calulaticn  does  not  take  into  effect  the  number  of  atoms  in 
each  ring.  Again,  the  predominance  of  rings  that  contain 
atoms  in  the  <100>  and  <101>   directions   is   evident.    The 
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peak  of  the  KSE-3  curve  was  closer  to  the  center  than  the 
Gibson-II  peaks,  which  were  spread  out  all  the  way  to  the 
lattice  edges.  Figure  17  shows  the  angular  dispersion 
results  for  (a)  Gibson-II  and  (b)  KSE-3.  The  particular 
seed  used  in  the  psuedo-r andom  number  generator  to  locate 
the  isotope  random  sites  was  such  that  only  light  isotopes 
sputtered  in  this  case. 

(2)  Neon  Ions .  The  results  when  100-eV  neon 
ions  impact  the  target  area  are  shown  in  figures  18,  19,  and 
20.  Figure  19  shows  that  for  the  light  mass  of  the  neon 
ion,  atoms  on  the  lattice  boundary  are  sputtered  using  (a) 
Gibson-II  and  (b)  KSE-3.  Gibson-II  resulted  in  a  sputtering 
ratio  of  0.889,  while  the  sputtering  ratio  for  KSE-B  was 
0.806.  Figure  19  shows  that  the  Ne  +  sputtered  atoms  were 
further  away  from  the  lattice  center  than  the  He+,  for  both 
(a)  Gibson-II  and  (b)  KSE-B.  This  was  confirmed  by  the 
increased  number  of  next-nearest  neighbor  atoms  (in  the 
<101>  direction)  sputtering.  Figure  20  indicates  that  some 
of  the  heavier  isotopes  sputtered  for  (a)  Gibson-II  and  (b) 
KSE-B.  There  were  not  enough  sputtered  atoms  to  notice  any 
trend  in  the  isotope  angular  dispersion. 

(3)  Argon  Ions^  The  argon  simulations  w°re 
used  to  investigate  the  "reflective  collision"  theory. 
Argon  ions  of  100-eV  energy  bombarded  the  target  area  using 
both  of  the  potential  functions.  In  the  random  isotope  site 
locations,  the  PKA  (atom  22)  was  a  heavy  isotope,  while 
three  (52,  58,  and  67)  of  the  four  atoms  directly  underneath 
were  light  isotopes.  Four  simulations  were  run  using  each 
of   the   potential   functions;  (1)  the  random  isotope  sites, 

(2)  forcing  the  PKA  to  a  light  isotope,  (3)  forcing  the  PKA 
to  a  heavy  isotope  and  all  four  atoms  underneath  to  a  light 
isotope,  and  (4)  forcing  the  PKA  to  a  light  isotope  and  all 
four  atoms  underneath  to  a  heavy  isotope.  The  remainder  of 
the  lattice  isotopes  were  unchanged. 
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Figure  21  shows  the  sputtering  summaries 
for  all  four  combinations  using  (a)  Gibson-II  and  (b)  KSE-B. 
In  all  four  cases,  identical  atoms  sputtered  independent  of 
isotope  locations.  Gibson-II  resulted  in  sputtering  ratios 
of  1.056,  while  KSE-B  resulted  in  sputtering  ratios  of 
0.944.  Figure  22  illustrates  the  ring  probabilty  in  all 
four  cases  for  (a)  Gibson-II  and  (b)  KSE-B.  Figure  23 
illustrates  the  angular  dispersion  of  ejected  atoms,  using 
Gibson-II,  for  (a)  case  (1),  (b)  case  (2),  (c)  case  (3),  and 
(d)  case  (4) .  No  light  isotopes  were  sputtered.  Figure  24 
illustrates  the  same  results  for  KSE-B. 

Since  at  least  four  atoms  were  ejected  into 
a  different  five  degree  window  for  Gibson-II,  the  sputtering 
mechanisms  were  checked.  The  most  active  event  occurred 
when  the  impact  area  was  struck  by  the  ion  at  impact  point 
4060  (impact  displaced  0.40  LU  in  the  X  direction  and  0.60 
LO  in  the  Z  direction  from  the  center  of  the  PKA,  hence 
4060)  .  Four  atoms  were  sputtered  at  this  point.  A 
single-impact  program  was  run  at  this  point  with  the  output 
listing  the  positions,  velocities,  kinetic  energy,  and 
potential  energy  of  each  atom  with  a  total  energy  greater 
than  0.1  eV.  The  listing  was  printed  once  every  five 
timesteps.  A  timestep  is  the  time  interval  between  dynamic 
calculations.  It  is  variable  and  determined  as  the  time 
reguired  for  the  fastest  moving  atom  to  move  0.1  LU  from  its 
previous  position.  In  real  time  the  timestep  is  on  the 
order  of  10~15  seconds.  Figure  25  illustrates  atom 
trajectories  in  the  sputtering  of  atoms  17,  32,  57,  and  77. 
In  cases  (3)  and  (4) ,  the  sputtering  mechanisms  were  the 
same  displaying  only  a  slight  difference  in  the  ejection 
angles  and  the  velocities  of  the  sputtered  atoms. 
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The  argon  ion  collides  with  atom  22  (PKA) 
and  atom  27  by  timestep  5,  pushing  them  down  and  under  atoms 
17  and  32  respectively.  The  ion  continued  downward  and 
collided  with  atom  67  by  timestep  15.  The  ion  reflected 
back  up  through  the  hole  7acated  by  atoms  22  and  27.  By 
timestep  15,  atoms  17  and  32  were  moving  upward.  Atom  67 
moved  down  colliding  with  atoms  107  and  112  by  timestep  20. 
Atoms  107  and  112  squeezed  [21]  under  atoms  102  and  117 
respectively.  Atoms  102  and  117  were  moving  up  by  timestep 
30,  and  struck  atoms  57  and  77  from  below  by  timestep  40. 
By  this  time,  atoms  17  and  32  had  moved  far  enough  above  the 
lattice  with  sufficient  upward  velocity  to  be  classified  as 
"sputs".  Atom  17  was  ejected  at  an  angle  of  33.56°  to  the 
surface  normal  with  an  upward  velocity  of  2924  m/sec  for 
case  (3) ,  and  at  an  angle  of  33.36°  with  a  velocity  of  2974 
m/sec  for  case  (4)  .  Atom  32  was  ejected  at  an  angle  of 
33.65°  with  an  upward  velocity  of  2926  m/sec  for  case  (3), 
and  at  an  angle  of  33.56°  with  a  velocity  of  2918  m/sec  for 
case  (4) .  Atoms  57  and  77  continued  upward  through  the 
holes  left  by  the  sputtered  atoms  17  and  32.  As  atom  57 
passed  through  this  hole  its  trajectory  was  bent  closer  to 
the  surface  normal  by  atoms  12,  16,  and  21.  Likewise,  atom 
77  was  reflected  upward  by  atoms  28,  33,  and  37.  By 
timestep  75,  atoms  57  and  77  were  far  enough  above  the 
lattice  surface  with  sufficient  upward  velocity  to  be 
classified  as  "sputs".  Atom  57  was  ejected  at  an  angle  of 
18.12°  with  an  upward  velocity  of  2967  m/sec  for  case  (3), 
and  at  an  angle  of  17.75°  with  a  velocity  of  2882  m/sec  for 
case  (4) .  Atom  77  sputtered  at  an  angle  of  25.04°  with  an 
upward  velocity  of  2799  m/sec  for  case  (3) ,  and  at  an  angle 
of  24.75°  with  a  velocity  of  2859  m/sec  for  case  (4). 
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It   is   noted  that  when  the  PKA  was  a  light 

isotope,  all  of  the  sputtered  atoms  were  ejected   closer  to 

the   surface   normal.   Also  noted  is  the  remarkable  symmetry 

of  ejection  angles  and  velocities  of  the  sputtered  atoms. 

(4)  Co£joer  Ions_j_  The  results  when  100-eV 
copper  ions  impact  the  target  area  are  shown  in  figures  26, 
27,   and   28.    Figure  26  shows  the  sputtering  summaries  for 

(a)  Gibson-II,  resulting  in  a  sputtering   ratio   0.944,   and 

(b)  KSE-B,  resulting  in  a  sputtering  ratio  of  1.111. 
Gibson-II  simulations  with  Ar+  and  Cu+  were  the  only  runs 
that  had  atoms  sputter  from  the  second  layer.  Also,  these 
two  ions  resulted  in  the  highest  sputtering  ratios  for  both 
potential  functions.  These  results  further  indicate  a 
maximum  sputtering  ratio  occurs  near  an  ion/atom  mass  ratio 
of  unity.  Figure  27  illustrates  a  continued  spread  outward 
of  the  ring  probabilty  for  (a)  Gibson-II  and  (b)  KSE-3. 
Figure  28  illustrates  the  angular  dispersion  results  for  (a) 
Gibson-II  and  (b)  KSE-B. 

(5)  E£I£ton  Ions_.  The  results  when  100-eV 
krypton  ions  impact  the  target  area  are  shown  in  figures  29, 
30,   and   31.    Figure  29  shows  the  sputtering  summaries  for 

(a)  Gibson-II,  resulting  in  a  sputtering  ratio  of  0.694,  and 

(b)  KSE-B,  resulting  in  a  sputtering  ratio  of  1.139.  Here, 
a  maximum  in  the  sputtering  ratio  was  reached  for  KSE-B, 
while  the  sputtering  ratio  for  Gibson-II  dropped 
drastically.  Figure  30  illustrates  a  continued  outward 
spread  of  the  ring  probabilities  for  (a)  Gibson-II  and  (b) 
KSE-B.  It  is  noticed  that  the  outward  spread  still  occured, 
even  though  the  Gibson-II  yield  decreased.  Figure  31 
illustrates  the  angular  dispersion  results  for  (a)  Gibson-II 
and  (b)  KSE-B. 
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(6)  Xenon  Igns^   The  results  when  100-eV   xenon 
ions  impact  the  target  area  are  shown  in  figures  32,  33f  and 
34.   Figure   32   shows   the   sputtering   summaries   for   (a) 
Gibson-II,  resulting  in  a  sputtering  ratio  of  0.556,  and  (b) 
KSE-B,  resulting  in  a  sputtering  ratio  of  1.028.   Figure   33 
illustrates  the  ring  probabilities  for  (a)  Gibson-II  and  (b) 
KSE-B.   The  outward  spread  of  the  ring  probability  peaks   is 
still   seen   for   KSE-B,   while  Gibson-II  remained  constant. 
Figure  34  illustrates  the  angular  dispersion  results  for  (a) 
Gibson-II  and  (b)  KSE-B. 

(7)  Gold  Ions^  The   results   when   100-eV   gold 
ions  impact  the  target  area  are  shown  in  figures  35,  36,  and 
37.   Figure   35   shows   the   sputtering   summaries   for   (a) 
Gibson-II,  resulting  in  a  sputtering  ratio  of  0.528,  and  (b) 
KSE-B,  resulting  in  a  sputtering  ratio  of  0.778.   Figure   36 
illustrates  the  ring  probabilities  for  (a)  Gibson-II  and  (b) 
KSE-B.   Gibson-II  has  settled  out  with   approximately   egual 
probability   of  sputtering  from  atoms  in  the  first  two  rings 
in  the  <  1 0 1 >  direction.   In  the  case  of  KSE-B,  a  shift   from 
the  first  and  second  rings  in  the  <100>  direction  toward  the 
second  ring  in  the  <10  1>  direction  is   noticed.    Figure   37 
illustrates  the  angular  dispersion  results  for  (a)  Gibson-II 
and  (b)  KSE-B. 

(8)  Mercury  Ions .  These  Hg+  bombardment 
results  are  compiled  from  the  sets  of  simulation  runs  used 
to  investigate  the  "reflective  collision"  theory.  Figure  38 
shows  the  sputtering  summaries  for  all  four  cases,  in  which 
the  masses  of  the  PKA  and  the  four  atoms  underneath  are 
varied.  The  results  for  (a)  Gibson-II  were  identical  in  all 
four  cases,  resulting  in  a  sputtering  ratio  of  0.528.  The 
sputtering  ratio  was  0.743  using  KSE-B  in  a  lattice 
consisting  of  atoms  with  only  the  natural  mass.  The  results 
for  (b)  KSE-3  depended  on  the  order  of  the  heavy  over   light 
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interaction.  When  the  PKA  was  a  heavy  isotope,  shown  in 
figure  38(b),  the  sputtering  ratio  was  0.7229  However,  when 
the  PKA  was  a  light  isotope,  atom  30  sputtered  an  additional 
time  (denoted  by  A  )  and  the  sputtering  ratio  became  0.750. 
Figure  39  illustrates  the  ring  probabilties  for  (a) 
Gibson-II,  which  remained  relatively  unchanged,  and  (b) 
KSE-B,  which  still  displayed  a  shift  to  the  second  ring  in 
the  101  direction.  Figure  40  illustrates  the  angular 
dispersion  results  for  Gibson-II.  Figure  41  shows  the 
angular  dispersion  results  for  KSE-B.  Again,  atom  30 
sputtered  an  additional  time.  The  additional  ejection  angle 
is  denoted  by  A . 

b.   1-keV  Results 

The  1-keV  events  were  obviously  not  contained  in 
the  9x4x9  microcrystallite .  However,  the  simulations  were 
run  to  compare  the  sputtering  ratios  with  experimental 
results.  Also,  the  dependence  of  the  sputtering  on  ion  mass 
was  investigated.  Useful  information  can  be  obtained  about 
the  ring  probabilties  and  isotope  angular  dipersion  at  this 
higher  energy.  The  162  atom,  9x4x9  microcrystallite,  with 
the  correct  ratio  of  copper  isotopes  at  random  sites,  was 
used.  Even  though  all  atoms  were  used  in  the  calculations 
and  the  printout  provided  information  on  all  atoms,  only  the 
top  layer  results  are  shown  in  the  figures.  Even  in  the 
worse  case  (Cu/Ar+,  Gibson-II),  more  than  86  per  cent  of 
the  sputtered  atoms  were  from  the  top  layer. 


The 


sputtering        simulations        for        normally 


incident,    1-keV    ions    were    run      using      He+,       Ne+,       Ar+, 


Cu+, 


Kr+,  Xe+,  Au  +  ,  and  Hg+  ions.  Figure  42  shows  the  sputtering 
ratio  as  a  function  of  the  ion  mass  for  both  potential 
functions.     As    expected,    the   sputtering   ratios   have 
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increased.  However,  they  are  no  longer  smooth  and 
continuous,  but  exhibit  a  randomness,  as  if  each  was 
independent  of  the  others. 

All  the  1-keV  simulations  used  the  same  random 
isotope  sites.  The  letters  L  and  H,  denoting  light  or  heavy 
isotopes,  appear  only  in  the  Cu/Ne+,  KSE-B  summary  (more 
different  atoms  sputtered  -  but  not  the  highest  sputtering 
ratio) . 

0)  Helium  Ions^  The  results  when  1-keV  helium 
ions  impact  the  target  area  are  shown  in  figures  43,  44,  and 
45.  Figure  43  shows  the  sputtering  summaries  for  (a) 
Gibson-II,  resulting  in  a  sputtering  ratio  of  1.914,  and  (b) 
KSE-B,  resulting  in  a  sputtering  ratio  of  2.833.  In  both 
cases,  boundary  atoms  sputtered.  Figure  44  illustrates  the 
ring  probabilities  for  (a)  Gibson-II  and  (b)  KSE-B.  Both 
display  a  large  action-reaction  effect,  as  if  a  bomb  were 
dropped  at  the  center  and  everything  out  from  the  center  was 
blasted  upward.  Figure  45  illustrates  the  angular 
dispersion  results  for  (a)  Gibson-II  and  (b)  KSE-B. 

(2)  Neon  Ions^  The  results  when  1-JceV  neon  ions 
impact  the  target  area  are  shown  in  figures  46,  47,  and  48. 
Figure  46  shows  the  sputtering  summaries  for  (a)  Gibson-II, 
resulting  in  a  sputtering  ratio  of  3.686,  and  (b)  KSE-B, 
resulting  in  a  sputtering  ratio  of  4.500.  Host  of  the 
surface  atoms  were  sputtered,  indicating  that  the  energy  ( 
or  momentum  )  spread  out  over  the  entire  lattice  surface. 
Figure  47  illustrates  the  ring  pr ooabilities  for  (a) 
Gibson-II  and  (b)  KSE-B.  The  KSE-3  simulation  indicates  a 
spreading  outward  of  the  sputtered  atoms.  The  Gibson-II 
simulation  not  only  exhibits  increased  sputtering  of  atoms 
further  from  the  center,  but  it  also  clearly  shows  a  shift 
in  the  number  of  sputtered  atoms  from  those  at  2.83  LU  in 
the  <101>  direction   to   those   at   3.16   LU   in   the   <103> 
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direction.  As  mentioned  earlier,  this  could  be  artificial 
because  the  ring  at  3.16  LO*  is  the  first  ring  to  contain 
eight  atoms.  Figure  48  illustrates  the  angular  dispersion 
results    for    (a)    Gibson-II    and    (b)     KSS-B. 

(3)         Argon      Ions^      The    results    when    1-keV   argon 
ions   impact    the   target   area    are    shown    in    figures    49,    50,    and 
51 .         Figure      49      shows      the      sputtering      summaries      for    (a) 
Gibson-II,    resulting    in   a    sputtering    ratio    of   3.806,    and    (b) 
KSE-3,      resulting    in    a   sputtering   ratio   of    4.611.       Figure   50 
illustrates    the   ring    probabilities    for    (a)     Gibson-II    and    (b) 
KSE-B.         The      Gibson-II      simulation   continued   to    increase    in 
the   outward    direction    at    the    expense    of   the    atoms   at    3.16    LU 
in      the      <103>      direction.         However,       the      KSE-B    simulation 
reversed    its    trend    and  the    ring    probability    for    the    close-in 
rings       (1.41    LU    and   2.00    LU)     increased    at    the   expense    of    the 
ring   at    2.83   LU    and   the   rings  further   out    from      the      center. 
Figure    5  1    illustrates    the    angular    dispersion    results    for    (a) 
Gibson-II   and    KSE-3. 

CO       Q.opp_er    Ions^     The   results    when    1-keV    copper 
ions    impact    the   target   area    are    shown    in   figures    52,    53,    and 
54.         Figure      52      shows      the      sputtering      summaries      for    (a) 
Gibson-II,    resulting    in   a   sputtering    ratio    of   3.028,    and    (b) 
KSE-3,      resulting    in   a    sputtering    ratio   of    3.972.       These    low 
sputtering      ratios      not      only      disagree      with      a         monatonic 
increase    in    the    yield    with    mass,    but    are    drastically    low    for 
the    maximum    to   occur   at   a      ion/atom      mass      ratio      of      unity. 
Figure        53      illustrates      the      ring      probabilities      for      (a) 
Gibson-II    and    (b)     KSE-B.       The   Gibson-II    simulation    showed   an 
ubexpected      change.       The      ring    probability    increased    for    the 
rings    at    2.00      LU       (<100>      direction)       and      3.16      LU       (<103> 
direction) .         This      trend      was      also      displayed    in    the   KSE-B 
simulation,      where      it      starts      to      occur      in      the  Cu/Ar+ 

simulation.  Figure  54  illustrates  the  angular  dispersion 
results    for     (a)    Gibson-II    and    (b)     KSE-B. 
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(5)  Krypton  Ions_.  The  results  when  1-keV 
krypton  ions  impact  the  target  area  are  shown  in  figures  55, 
56,   and   57.    Figure  55  shows  the  sputtering  summaries  for 

(a)  Gibson-II,  resulting  in  a  sputtering  ratio  of  2.833,  and 

(b)  KSE-B,  resulting  in  a  sputtering  ratio  of  3.917.  Note 
that  the  sputtering  ratio,  for  both  potential  functions, 
continues  to  decrease  with  increased  ion  mass.  Figure  56 
illustrates  the  ring  probabilities  for  (a)  Gibson- II  and  (b) 
KSE-B.  Both  potential  function  simulations  continue  to 
exhibit  the  trend  to  establish  a  peak  at  2.00  LU  in  the 
<100>  direction  and  at  3.16  LU  in  the  <103>  direction. 
However,  the  KSE-B  simulation  shows  significantly  larger 
ring  probability  at  3.16  LU.  Figure  57  illustrates  the 
angular  dispersion  results  for  (a)  Gibson-II  and  (b)   KSE-B. 

(6)  Xe_H2n_   Ions^   The  results  when  1-keV  xenon 
ions  impact  the  target  area  are  shown  in  figures  58,  59,  and 
60.    Figure   58   shows   the   sputtering   summaries   for  (a) 
Gibson-II,  resulting  in  a  sputtering  ratio  of  2.333,  and  (b) 
KSE-B,   resulting   in   a   sputtering   ratio   of  3.472.   3oth 
potential  functions  continued  to   show   a   decrease   in   the 
sputtering   ratio   with   an  increase  in  the  mass  of  the  ion. 
Figure   59   illustrates   the   ring   probabilties    for    (a) 
Gibson-II  and  (b)  KSE-B.   All  the  previous  trends,  except  an 
increased   ring   probability   at   3.16   LU   were   shattered. 
Figure  60  illustrates  the  angular  dispersion  results  for  (a) 
Gibson-II  and  (b)  KSE-B. 

(7)  Gold  Ions^  The  results  when  1-keV  gold 
ions  impact  the  target  area  are  shown  in  figures  61,  62,  and 
63.  Figure  61  shows  the  sputtering  summaries  for  (a) 
Gibson-II,  resulting  in  a  sputtering  ratio  of  2.778,  and  (b) 
KSE-3,  resulting  in  a  sputtering  ratio  of  4.167.  A  definite 
increase  in  the  sputtering  ratio  occured.  Figure  62 
illustrates  the  ring  probabilities  for  (a)  Gibson-II  and  (b) 
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KSE-B.  No  new  trends  were  noted,  except  that  the  ring  at 
3.16  Lu"  continues  to  dominate  and  the  other  ring 
probabilities  fluctuate.  Figure  63  illustrates  the  angular 
dispersion  results  for  (a)  Gibson-II  and  (b)  KSE-3. 

*  (8)    Mercury   Ions.    The   results   when   1-keV 

mercury  ions  impact  the  target  area  are  shown  in  figures  64, 
65,  and  66.   Figure  64  shows  the   sputtering   summaries   for 

(a)  Gibson-II,  resulting  in  a  sputtering  ratio  of  2.778,  and 

(b)  KSE-B,  resulting  in  a  sputtering  ratio  of  4.278.  These 
simulations  resulted  in  no  change  in  the  sputtering  ratio 
for  Gibson-II,  and  only  a  slight  increase  for  KSE-B.  Figure 
65  illustrates  the  ring  probabilities  for  (a)  Gibson-II  and 
(b)  KSE-B.  Only  insignificant  changes  occured  compared  to 
the  gold  icn  results.  Figure  66  illustrates  the  angular 
dispersion  results  for  (a)  Gibson-II  and  (b)  KSE-B. 

c.   2-keV  Results 

The  2-keV  ion  simulations  were  the  earliest  runs 
made.  No  isotope  information  was  obtained.  As  a 
consequence,  the  only  results  are  the  sputtering  ratio 
versus  ion  mass  curve,  figure  67,  for  both  potential 
functions.  They  exhibit  an  even  more  erratic  behavior  than 
did  the  1-keV  simulations.  The  sputtering  summaries 
frequently  show  atoms  sputtered  from  the  third  layer,  and 
rarely  an  atom  sputtered  from  the  fourth  layer.  At  certain 
impact  points  as  many  as  25  atoms  sputtered.  The  small, 
9x4x9  microcrystallite  was  effectively  blasted  apart  by  the 
2-keV  ions. 
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2-  Ill CL  Orientation  100  eV 

The  effort  devoted  to  (110)  orientation  simulation 
runs  was  very  limited  due  to  the  higher  priority  assigned  to 
the  (100)  orientation  runs.  As  a  result,  only  the  argon  and 
mercury  ion  simulations  at  100  eV  were  run.  Still  needing 
investigation  are  the  sputtering  ratio  as  a  function  of  ion 
mass  and  sputtering  mechanism  traces. 

The  (110)  orientation  microcrystalli te  was 
constructed  in  a  lattice  consisting  of  123  copper  isotopes. 
In  the  random  isotope  site  locations,  the  PKA  (14)  was  a 
light  isotope,  as  were  three  (33,  36,  and  37)  of  the  four 
atoms  underneath  the  PKA.  Figure  4  shows  an  isometric 
projection,  of  the  lattice  and  its  numbering  scheme.  Figure 
10  shows  an  expanded  view  of  the  grid  of  36  impact  points 
used.  A  9x6x9  lattice  was  used  with  no  increase  in  computer 
time  required  to  complete  a  run. 

Two  sets  of  simulations  were  run  to  investigate  the 
"reflective  collision"  theory.  Four  runs  were  made  in  each 
set  of  simulations  for  each  potential  function.  The  impact 
area  was  bombarded  by  100-eV,  argon  and  mercury  ions.  The 
four  runs  consisted  of;  (1)  the  random  isotope  sites  (PKA 
was  a  light  isotope) ,  (2)  forcing  the  PKA  to  a  heavy 
isotope,  (3)  forcing  the  PKA  to  a  light  isotope  and  all  four 
atoms  underneath  to  a  heavy  isotope,  and  (4)  forcing  tha  PKA 
to  a  heavy  isotope  and  all  four  atoms  underneath  to  a  light 
isotope.  The  remainder  of  the  lattice  isotopes  were 
unchanged. 
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a.   Argon  Ions 

The  results  when  100-eV  argon  ions  impact  the 
target  area  are  shown  in  figures  68-75.  Figure  63  shows  the 
sputtering  summaries  using  Gibson-II  for  all  four  cases. 
Case  (1)  and  case  (3)  ,  shown  in  figure  63(a)  ,  resulted  in  a 
sputtering  ratio  of  0.839.  Case  (2)  and  case  (4),  shown  in 
figure  68(b),  resulted  in  a  sputtering  ratio  of  0.361.  A 
predominance  of  sputtered  atoms  in  the  <100>  and  < 10 1> 
directions  is  evident.  Figure  69  illustrates  the  ring 
probabilities  for  (a)  case  (1)  and  case  (3),  and  (b)  case 
(2)  and  case  (4)  using  Gibson-II.  Notice  that  an  atom  at 
the  boundary  (  in  the  <101>  direction  )  sputtered.  Figure 
70  illustrates  the  angular  dispersion  results  for  (a)  case 
(1),  and  (b)  case  (3)  using  Gibson-II.  Figure  71  shows  the 
same  information  for  (a)  case  (2)  ,  and  (b)  case  (4)  . 

Figure  72  shows  the  sputtering  summaries  using 
KSE-B  for  all  four  cases.  Case  (3)  resulted  in  a  sputtering 
ratio  of  0.750,  while  case  (1)  sputtered  atom  26  one 
additional  time  (denoted  by  A )  and  had  a  sputtering  ratio 
of  0.778.  Case  (1)  and  case  (3)  are  shown  in  figure  72(a). 
Case  (2)  resulted  in  a  sputtering  ratio  of  0.778,  while  case 
(4)  sputtered  atom  8  one  additional  time  and  had  a 
sputtering  ratio  of  0.805.  Case  (2)  and  case  (4)  are  shown 
in  figure  72(b) .  Again,  it  is  evident  that  only  atoms  in 
the  <100>  and  <  1 1 0  >  directions  sputtered.  Figure  ~?3 
illus-crates  the  ring  probabilities,  using  KSF-B,  for  (a) 
case  (1)  and  case  (3)  ,  and  (b)  case  (2)  and  case  (4)  .  The 
close-in  trend,  similar  to  that  of  Gibson-II,  is  evident,  as 
are  the  sputtered  atoms  at  the  boundary  in  the  <101> 
direction.    Figure   74   illustrates   the  angular  dispersion 
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results,  using  KSE-B,  for  (a)  case  (1)  ,  and  (b)   case   (3)  . 
Figure   75   shows  the  same  information  for  (a)  case  (2),  and 
(b)  case  (4)  . 

b.   Mercury  Ions 

The  results  when  100-eV  mercury  ions  impact  the 
target  area  are  shown  in  figure  76,  77,  and  78.  Figure 
76(a)  shows  the  sputtering  summaries  for  all  four  cases, 
using  Gibson-II.  Case  (1),  case  (2),  and  case  ('4)  resulted 
in  a  sputtering  ratio  of  0.278.  Case  (3)  sputtered  atom  13 
one  additional  time  (denoted  by  A )  and  had  a  sputtering 
ratio  of  3.05.  Figure  76(b)  shows  the  sputtering  summaries 
for  all  four  cases,  using  KSE-B,  which  resulted  in 
sputtering  ratios  of  0.222.  Figure  77  illustrates  the  ring 
probabilities  for  (a)  Gibson-II  and  (b)  KSE-B.  Figure  78 
illustrates  the  angular  dispersion  rasults  for  (a)  Gibson-II 
and  (b)  KSE-B. 

No  significant  trends  were  noted  in  the  (110) 
orientation.  However,  the  ease  with  which  atoms  sputtered 
in  the  <101>  direction  is  evident. 

3  .    J1U1  0£L§nta tion  100e V 

Like  the  (110)  orientation,  the  (111)  orientation 
did  not  receive  extensive  attention.  Both  the  sputtering 
ratio  as  a  function  of  the  ion  mass,  and  sputtering 
mechanism  traces  need  further  investigation. 

The  (111)  orientation  microcrystallite  is 
constructed  in  a  9x4x9  lattice  consisting  of  163  atoms.  In 
the  random  isotope  site  locations,  the  PKA  (22)  was  a  light 
isotope,   as   were   three  (58,  63,  and  67)  of  the  four  atoms 
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underneath  the  PKA.  Figure  5  shows  an  isometric  projection 
of  the  lattice  and  its  numbering  scheme.  Figure  11  shows  an 
expanded  view  of  the  grid  of  21  impact  points  used. 

Two  sets  of  simulations  were  run  to  investigate  the 
"reflective  collision"  theory.  Four  runs  were  made  in  each 
set  of  simulations  for  each  potential  function.  The  target 
was  bombarded  by  100-eV  argon  and  mercury  ions.  The  four 
runs  consisted  of;  (1)  the  random  isotope  sites  (PKA  was  a 
light  isotope),  (2)  forcing  the  PKA  to  a  heavy  isotope,  (3) 
forcing  the  PKA  to  a  light  isotope  and  all  four  atoms 
underneath  to  a  heavy  isotope,  and  (4)  forcing  the  PKA  to  a 
heavy  isotope  and  all  four  atoms  underneath  to  a  light 
isotope.  The  remainder  of  the  lattice  isotopes  rfere 
unchanged. 

a.   Argon  Ions 

The  results  when  100-eV  argon  ions  impact  the 
target  area  are  shown  in  figures  79-83.  Figure  79  shows  the 
sputtering  summaries  for  (a)  Gibson-II,  resulting  in 
sputtering  ratios  of  0.857,  and  (b)  KSE-B,  resulting  in 
sputtering  ratios  of  0.762.  For  KSE-B,  case  (2) ,  atom  16 
sputtered  one  additional  time  (denoted  by  A),  resulting  in 
a  sputtering  ratio  of  0.310.  Figure  80  illustrates  the  ring 
probabilities  for  (a)  Gibson-II  and  (b)  KSE-B.  It  was 
unexpected  that,  for  KSE-B,  all  of  the  sputtered  atoms  are 
in  the  <101>  direction  in  the  ring  at  2.83  LU.  Figure  31 
illustrates  the  Gibson-II  angular  dispersion  results  for  (a) 
case  (1)  and  case  (3)  ,  and  (b)  case  (2)  and  case  (4)  . 
Figure  82  illustrates  the  KSE-B  angular  dispersion  results 
for  (a)  case  (1)  and  (b)  case  (3)  .  Figure  83  shows  the  same 
results  for  (a)  case  (2)  ,  and  (b)  case  (4)  . 
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b.   Mercury  Ions 

The  results  when  100-eV  niercury  ions  impact  the 
target  area  are  shown  in  figures  84-89.  Figure  84  (a)  shows 
the  sputtering  summaries,  using  Gibson-II,  for  case  (1), 
case  (2) ,  and  case  (4) ,  which  resulted  in  sputtering  ratios 
of  0.476.  Case  (3)  had  a  sputtering  ratio  of  0.381  and  is 
shown  in  figure  84(b)  .  Figure  85  shows  the  sputtering 
summaries  using  the  KSE-B.  Case  (3)  resulted  in  a 
sputtering  ratio  of  0.762,  while  case  (1)  sputtered  two 
additional  atoms  (12  and  14  -  denoted  by  A)  and  has  a 
sputtering  ratio  of  0.857.  Case  (1)  and  case  (3)  are  shown 
in  figure  85(a)  .  Case  (2)  and  case  (4)  ,  shown  in  figure 
85(b),  have  sputtering  ratios  of  0.952.  Figure  86 
illustrates  the  ring  probabilities  for  (a)  Gibson-II  and  (b) 
KSE-B.  Again,  it  is  evident  that  the  only  sputtered  atoms 
are  in  the  ring  at  2.83  LU,  in  the  < 1 0 1 >  direction.  Figure 
87  illustrates  the  Gibson-II  angular  dispersion  results  for 
(a)  case  (1)  ,  case  (2)  ,  and  case  (4)  ,  and  (b)  case  (3)  . 
Figure  88  illustrates  the  KSE-B  angular  dispersion  results 
for  (a)  case  (1)  and  (b)  case  (3)  .  Figure  89  shows  the  same 
information  for  (a)  case  (2)  ,  and  (b)  case  (4)  . 


B.   LARGER  LATTICES 


The  erratic  behavior  of  the  sputtering  ratio  as  a 
function  of  the  ion  mass,  indicates  that  the 
microcrystallite  is  not  containing  the  event.  Originally, 
"containment"  was  defined  as  "no  boundary  atoms  sputtering." 
Based  on  this  definition,  none  of  the  100-eV  events  were 
contained  in  the  (100)  orientation  or  (110)  orientation  9x9 
microcrystallites .  However,  in  the  (111)  orientation, 
100-eV  events  were  contained.   The  sputtering  summaries   for 
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mercury  ions  bombarding  a  13x6x13,  (110)  orientation 
microcrystallite  containing  255  atoms  are  shown  in  figure 
90.  Figure  90(a)  is  for  Gibson-II,  while  figure  90(b)  is 
for  KSE-B.  Comparison  of  figure  90  with  figure  76, 
confirmed  that  the  9x5x9  microcrystallite  did  not  contain 
the  100-eV,  mercury  event.  Based  on  the  original  definition 
of  containment,  the  13x6x13  microcrystallite  did  not  contain 
the  Gibson-II  simulation  run. 

In  the  (111)  orientation,  the  sputtering  summaries  for 
mercury  ions  bombarding  a  13x4x13  microcrystallite 
containing  339  atoms  are  shown  in  figure  91.  Figure  91(a) 
is  for  Gibson-II,  while  figure  91(b)  is  for  KSE-3. 
Comparison  of  figure  91  with  figures  84  and  85,  indicate 
that  even  though  the  9x4x9  lattice  had  no  boundary  atoms 
sputtered,  the  sputtering  ratio  increases  when  a  13x4x13 
lattice  is  used. 

1 •   100-eV  Containment 

Based  on  the  above  information,  the  definition  for 
"containment"  was  changed  to:  "an  event  is  contained  when  an 
increase  in  the  microcrystallite  size  does  not  result  in  any 
further  increase  in  the  yield."  Earlier  simulations 
resulted  in  larger  yields  for  KSE-B.  For  this  reason,  KSE-B 
was  used  to  investigate  the  size  microcrystallite  needed  to 
contain  a  100-eV  event. 

The  sputtering  ratio  as  a  function  of  the  ion  mass 
was  examined  for  the  (100)  orientation.  Figure  92  shows  the 
results  for  9x4x9,  11x4x11,  and  11x6x11  microcrystallites . 
Additionally,  simulations  were  run  for  copper  (largest 
sputtering  ratio)  and  mercury  ions  for  13x4x13,  13x6x13,  and 
17x4x17  microcrystallites.  Based  on  the  new  definition  of 
containment,  100-eV  helium  and  xenon  events  are  contained  in 
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a  9x4x9  microcrystallite  because  an  11x4x11  microcrystallite 
produces  the  same  yield.  However,  the  argon,  copper,  and 
krypton  events  were  not  contained.  Further  increase  of  the 
microcrystallite  size  demonstrated  that  the  100-eV,  copper 
event  is  contained  by  a  13x4x13  microcrystallite.  Figure  92 
also  illustrates  that  the  maximum  sputtering  ratio  occurs 
when  the  ion/atom  mass  ratio  is  approximately  unity.  Figure 
93  shows  the  same  information  in  a  plot  of  sputtering  ratio 
versus  lattice  size  for  ill  eight  ions. 

2  •  llx.l_x.ll   microcrystallite 

Asuming  that  all  100-eV  events  are  contained  by  a 
17x4x17  microcrystallite,  the  effect  of  increasing  the  ion 
energy  was  investigated.  The  reason  for  this  investigation 
was  to  determine  how  energetic  an  event  could  be  contained 
by  a  17x4x17  microcrystallite.  And  to  compare  higher  energy 
events  with  experimental  results.  The  reason  for  not 
increasing  the  lattice  size  without  limit  was  obviously  to 
keep  the  computer  run  time  within  reason.  Again,  only  KSE-B 
with  mercury  and  copper  ions  were  investigated. 

a.   Copper  Ions 

The  results  when  100-eV  copper  ions  impact  the 
target  area  are  shown  in  figures  94,  95,  and  96.  Figure  94 
shows  the  sputtering  summary  with  a  sputtering  ratio  of 
1.389.  Containment  for  the  100-eV  copper  event  had  already 
been  established.  Figure  95  illustrates  the  "atom 
sputtering  probability",  henceforth  called  the 
"probability".  The  "probability"  is  the  "ring  sputtering 
probability"  divided  by  the  number  of  atoms  in  each  ring. 
The  dimensions  for  "probability"  is  sputs  per  atom  divided 
by   the  total  sputs  at  a  fixed  distance  from  the  impact  area 
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reference  point.  Figure  95  indicates  an  enhancement  of 
sputtered  atoms  in  the  <100>  and  < 1 0 1 >  directions.  Figure  96 
illustrates  the  angular  dispersion  results. 

Results  for  the  200-eV  copper  ion  simulations 
are  shown  in  figures  97,  98,  and  99.  Figure  97,  the 
sputtering  summary  with  a  sputtering  ratio  of  2.839,  shows 
that  boundary  atoms  are  sputtered,  and  that  the  event 
definitely  is  not  contained.  Figure  98  illustrates  the 
"probability"  results  and  indicates  that  an  atom  farther  out 
had  increased  chances  of  sputtering.  Figure  99  illustrates 
the  angular  dispersion  results.  A  bias  toward  light 
isotopes  closer  to  the  normal  is  evident. 

The  results  for  the  500-eV  copper  ion 
simulations,  with  a  sputtering  ratio  of  4.028,  are  shown  in 
figures  100,  101,  and  102.  Again,  boundary  atoms  are 
sputtered  and  a  slight  shift  outward  of  the  "probability"  is 
noticed.  Figure  102  illustrates  the  angular  dispersion 
results  which  showed  that  more  light  isotopes  are  ejected 
more  nearly  normal  to  the  lattice  surface. 

the  1-keV,  copper  ion  results  with  a  sputtering 
ratio  of  7.083,  are  shown  in  figures  103,  104,  and  105. 
Figure  103,  the  sputtering  summary,  shows  a  very  large 
number  of  the  surface  atoms  sputtered.  Boundary  atoms  in 
every  direction  are  sputtered.  This  increased  spreading 
outward  was  confirmed  in  figure  104,  which  shows  the 
"probabilities"  relatively  level.  Figure  105  illustrates 
the  angular  dispersion  results. 

Obviously,  the  200-eV,  500-eV,  and  1-keV  copper 
events  are  not  contained,  even  in  a  17x4x17 
microcrystallite. 
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b.   Mercury  Ions 

The  results  for  the  100-eV  mercury  ion 
simulations  are  shown  in  figures  1 06 ,  107,  and  108.  Again, 
the  11x4x11  lattice  has  already  been  determined  to  contain 
the  100-eV,  mercury  event.  Figure  106  shows  the  sputtering 
summary  with  a  sputtering  ratio  of  0.833.  Figure  107 
illustrates  the  "probability",  and  like  the  copper 
simulation,  was  predominant  in  the  < 1 0 1 >  direction.  Figure 
108  illustrates  the  angular  dispersion  results. 

The  200-eV  mercury  ion  results,  with  a 
sputtering  ratio  of  1.722,  are  shown  in  figures  109,  110, 
and  111.  Figure  109,  the  sputtering  summary,  shows  an 
outward  spreading  of  the  sputtered  atoms.  However,  no 
boundary  atoms  were  sputtered  and  the  200-eV,  mercury  event 
is  tentatively  classified  as  "contained".  Figure  110 
confirms  the  outward  shift  in  "probability".  Figure  111 
illustrates  the  angular  dispersion  results,  which  indicate  a 
greater  number  of  heavy  isotopes  sputtering  than  the  isotope 
ratio  would  justify. 

The  500-eV  mercury  ion  results,  with  a 
sputtering  ratio  of  3.667,  are  shown  in  figures  112,  113, 
and  114.  Figure  112,  the  sputtering  summary,  indicates  that 
boundary  atoms  are  sputtered.  Figure  113  illustrates  the 
continued  outward  spreading  of  sputtered  atoms.  Figure  114 
illustrates  the  angular  dispersion  results,  and  the  number 
of  heavy  to  light  isotopes  closer  to  the  correct  ratio. 

The  1-keV  mercury  ion  results,  with  a  sputtering 
ratio  of  6.139,  are  shown  in  figures  115,  115,  and  117. 
Figure  115  shows  that  boundary  atoms  are  sputtered.    Figure 

116  shows  the  "probability"  still  spreading  outward.   Figure 

117  illustrates  the  angular  dispersion  results. 
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The  "probability"  figures  for  both  the  copper 
and  mercury  ions  are  remarkably  similar  and  display  the  same 
trends. 

c.   Sputtering  Mechanism  Traces 

The  containment  studies  indicated  that,  as  the 
lattice  size  was  increased,  the  sputtering  ratio  increased. 
Investigation  into  the  similarities  of  the  17x4x17  and  9x4x9 
microcrystallites,  indicated  that  atom  108  (17x4x17) 
sputtered  twice  in  36  events.  This  corresponds  to  atom  40 
in  the  9x4x9  microcrystallite,  which  did  not  sputter  in  36 
events.  Sputtering  mechanism  trace  runs  were  performed  to 
determine  the  sputtering  mechanism  involved  in  the 
sputtering  of  atom  108.  Atom  108  sputtered  when  ions 
impacted  at  points  50  (  impact  point  displaced  0.50  LU  in 
the  Z  direction  only)  and  60  (  impact  point  displaced  0.6  LU 
in  the  Z  direction  only)  . 

Figure  118  illustrates  the  atom  trajectories  in 
the  sputtering  of  atom  108  from  impact  point  50  (  the 
mechanism  is  the  same  for  impact  point  60) .  The  copper  ion 
struck  atoms  74  (PKA)  and  91  forcing  them  down  by  timestep 
15.  As  atom  91  was  moving  downward,  atoms  99,  100,  and  108 
prevented  atom  91  flying  off  in  the  positive  Z  direction, 
and  channeled  it  downward.  3y  timestep  25,  atom  91  had 
collided  with  atom  224,  reversing  its  direction.  At 
timestep  30,  atom  91  was  forcing  atom  108  upward  from 
underneath.  Atom  91  transfered  almost  all  of  its  energy 
(momentum)  to  atom  108  and  had  less  than  0.1  eV  by  timestep 
45.  As  atom  108  was  going  upward,  it  also  had  a  velocity  in 
the  positive  Z  direction  of  greater  than  6000  m/sec.  Atoms 
116,  117,  and  125  bent  atom  108  upward  quite  radically.  The 
velocity  in  the  Z  direction  had  been  reduced  to  less  than 
700  m/sec  by  timestep  50.   By  timestep  60,  atom  108  was   far 
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enough  above  the  lattice  surface,  with  sufficient  velocity 
to  be  classified  as  a  sput.  Atom  108  was  ejected  at  an 
angle  of  19.23°  to  the  surface  normal  with  an  upward 
velocity  of  3206  m/sec. 

Quite  evident  is  the  role  atoms  116,  117,  and 
125  have  in  channeling  atom  108  upward.  In  the  9x4x9 
lattice,  atom  40  (  a  boundary  atom)  corresponds  to  atom  108. 
And  nothing  equivalent  exists  for  atoms  116,  117,  and  125. 
Investigation  was  then  made  to  determine  whether  an  atom 
corresponding  to  atom  103  could  be  sputtered  in  the  9x4x9 
lattice  by  using  the  same  impact  point  on  a  different  PKA. 
In  this  simulation,  atom  31,  corresponding  to  atom  109  in 
the  larger  lattice,  was  expected  to  sputter.  Figure  119 
illustrates  the  atom  trajectories  when  the  copper  ion  struck 
the  new  PKA  (atom  13)  at  impact  point  50.  Atom  31  was 
ejected  at  an  angle  of  20.17°  with  an  upward  velocity  of 
3127  m/sec.  The  mechanism  is  the  same,  and  the  importance  of 
boundary  atoms  is  very  evident. 

Additional  runs  were  made  on  the  9x4x9  lat-tice 
bombarding  atom  13  at  impact  points  50  and  60.  At  both 
points,  comparison  is  made  between  runs  with  (1)  only  atoms 
of  the  natural  mass,  ,  and  natural  mass,  and  (2)  forcing 
atom  22  (91)  to  a  heavy  isotope  and  atom  67  (224)  to  a  light 
isotope.  In  all  four  cases,  the  sputtering  mechanism  is  the 
same.  The  differences  in  the  ejection  angles  and  upward 
velocities,  for  each  impact  point,  are  less  than  0.4°  and  75 
m/sec  respectively. 

3 •   Ion/Atom  Mass  Ratio 

Sigmund  [28]  proposed  a  theory  that  the  maximum 
spattering  yield  increases  as  a  function  of  energy. 
However,   results   of  this  investigation  strongly  indicate  a 
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dependence  on  the  ratio  of  the  ion  mass  to  the  targe-c  atom 
mass.  A  rearrangement  of  Sigmund's  theoretical  curves 
(figure  16b,  ref.  28)  to  show  the  sputtering  ratio  as  a 
function  of  the  ion/atom  mass  ratio  results  in  figure  120. 
The  curves  shown  by  Sigmund  are  predicted  for  200-eV, 
600-eV#  1-keV,  and  5-keV  ions  bombarding  an  aluminum  target. 
Figure  120  illustrates  that  Sigmund's  theory  would  predict 
the  sputtering  ratio  to  increase  as  a  function  of  the 
ion/atom  mass  ratio.  Figure  121  shows  experimental  results 
obtained  by  Wehner  and  is  presented  by  Sigmund  (  figure  12, 
ref.  28  )  to  support  his  theory.  This  figure  has  been 
replotted  to  show  the  sputtering  ratio  as  a  function  of  the 
ion/atom  mass  ratio.  The  results  were  obtained  by 
bombarding  a  copper  target  with  100-eV,  200-eV,  and  600-eV 
ions.  Experimental  results  plotted  in  this  manner  show  a 
decrease  in  the  sputtering  ratio  with  an  ion/atom  mass  ratio 
greater  than  unity,  in  disagreement  with  the  theory  proposed 
by  Sigmund.  Results  obtained  in  the  present  investigation 
completely  agree  with  the  shape  of  experimental  data  shown 
in  figure  121. 


C.   RANDOM  IMPACT  POINTS 


In  response  to  other  investigators'  suggestion  that  the 
target  impact  area  be  bombarded  randomly  rather  than  the 
systematic  grid  of  impact  points,  the  computer  program  was 
slightly  modified.  Since  additional  data  were  desired  for 
mercury  events,  Hg+  was  used  as  the  incident  ion.  The 
KSE-B  potential  function  was  used.  Earlier  results  revealed 
that  the  100-eV,  mercury  ion  event  was  contained  in  an 
11x4x11  raicrocrystallite. 
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Figure  122  illustrates  the  (100)  orientation  impact  area 
that  was  struck  at  the  144  randomly  selected  impact  points. 
A  few  gaps  are  evident,  but  the  randomness  is  very  apparent'. 
Figures  123-128  show  a  comparison  of  the  results  from  a 
randomly  bombarded  impact  area  and  the  results  from  the  grid 
of  impact  points.  Since  the  number  of  events  were  not  the 
same,  the  sputtering  summaries  indicate  a  sputtering  ratio 
for  each  atom.  Figure  123  shows  the  normalized  sputtering 
summary  for  the  randomly  selected  impact  points.  Figure  124 
shows  the  same  information  using  the  grid  of  impact  points. 
The  four  groups  of  36  randomly  selected  impact  points  had  a 
mean  sputtering  ratio  of  0.812  with  a  variance  of  0.0048. 
Additional  runs  (  not  shown  ),  combined  with  the  four 
previous  groups  to  give  a  total  of  21  groups  of  36  randomly 
selected  impact  points  (  756  shots  ) ,  resulted  in  a  mean 
sputtering  ratio  of  0.782  with  a  variance  of  0.016.  This 
compares  very  favorably  with  a  sputtering  ratio  of  0.778 
using  the  grid  of  impact  points.  Figures  125  and  126  show 
the  "probability"  results  for  the  randomly  selected  impact 
points  and  the  grid  of  impact  points  respectively.  It  is 
quite  evident  that  they  are  almost  identical.  Figure  127 
illustrates  the  angular  dispersion  results  for  randomly 
selected  impact  points.  A  definite  shift  of  the  lighter 
isotopes  toward  the  surface  normal  is  observed.  Figure  128 
illustrates  the  same  information  for  the  grid  of  impact 
points.  The  absence  of  the  shift  of  the  light  isotopes 
toward  the  normal  is  noted. 
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V.   CONCLUSIONS  AND  RECOMMENDATIONS 


A.   SIMULATION  MODEL 


The  more  sophisticated  a  simulation  model  becomes,  the 
more  uncertain  are  its  results.  Some  of  these  uncertainties 
are  now  beginning  to  appear  in  this  model. 

The  lattice  generators  can  now  construct  randomly 
positioned  isotopes.  Symmetrical  unfolding  of  the  event 
results  is  the  reverse  of  the  folding  required  to  reduce  the 
representative  impact  area.  The  model  can  now  display  the 
ejection  angle  of  the  sputtered  isotopes  in  a  histogram  for 
visual  comparison.  It  now  determines  the  distance  from  the 
surface  layer  center  to  each  ring  of  sputtered  atoms,  the 
number  of  atoms  and  the  number  of  times  each  atom  sputters 
in  each  ring,  and  the  probability  that  that  each  ring 
contained  an  atom  that  was  sputtered.  Other  modifications 
were  made  in  the  model  to  approximate  more  nearly  a 
meaningful  physical  decription.  These  include  changes  in 
the  distances  from  the  lattice  surface  at  which  calculations 
are  truncated,  the  depth  of  lattice  needed  to  prevent  any 
change  in  the  yield,  and  a  systematic  determination  of  the 
distance  above  the  surface  at  which  an  atom  is  declared  to 
have  sputtered. 
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B.   CONTAINMENT 


The  containment  of  an  event  is  paramount  in  a  complete 
description  of  the  sputtering  process.  However,  it  is  less 
significant  when  studying  mechanisms  and  trends. 

3ased  on  the  definition  of  containment  (  no  increase  in 
yield  with  an  increase  in  the  microcrystallite  size  ) ,  the 
100-eV  results  showed  that  a  9x4x9,  (100)  orientation 
microcrystallite  contained  only  the  helium  and  xenon  events. 
The  100-eV,  mercury  event  is  contained  in  an  1 1xUx11 
lattice.  To  contain  a  100-eV  copper  event,  which  had  the 
largest  yield,  a  13x4x13  lattice  was  needed.  Sputtering 
mechanism  traces  revealed  the  importance  of  boundary  atoms 
in  deflecting  and  channeling  atoms  in  the  sputtering 
process.  By  choosing  the  PKA  to  be  an  atom  other  than  the 
one  at  the  surface  layer  center,  the  effect  of  an  increased 
lattice  size  could  be  duplicated. 

Investigation  with  a  17x4x17  microcrystallite  revealed 
that  Gibson-II  potential  function  simulations  are  more 
difficult  to  contain  than  when  using  the  KSE-B  potential 
function.  Only  the  100-eV,  copper  event  was  contained  using 
Gibson-II,  while  using  KSE-B  the  100-eV  and  200-eV,  copper 
events  were  contained.  It  is  quite  obvious  that  to  contain 
any  higher-energy  events  would  require  a  microcrystallite  so 
large  that  the  computer  run  time  would  be  prohibitive. 
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C.   ANGULAR  DISPERSION 


Because  the  low-energy  events  resulted  in  so  few 
sputtered  atoms,  the  angular  dispersion  results  for  any 
particular  ion  are  inconclusive.  A  trend  of  light  isotope 
enrichment  is  observed  when  the  results  from  numerous  runs 
are  combined.  For  example,  the  histograms  of  the  four  sets 
of  36  randomly  selected  impact  points  show  a  definite  shift 
of  the  lighter  isotopes  to  smaller  ejection  angles.  In  his 
report  on  the  enrichment  of  lighter  isotopes  in  the  normal 
direction,  Wehner  [24]  noted  only  a  0.6%  enrichment  in  Cu 
(63) .  The  number  of  sputtered  isotopes  was  not  sufficient 
to  detect  this  small  difference. 

The  "reflective  collision"  theory,  which  Wehner  suggests 
to  be  the  cause,  was  observed  to  have  very  little  effect  on 
the  sputtered  atoms.  Investigation  revealed  that  in  most  of 
the  simulation  runs,  the  same  atoms  sputtered  independent  of 
its  isotope  type  or  the  various  positions  of  the  lattice 
isotopes.  When  the  PKA  was  a  light  isotope,  the  atoms 
that  did  sputter  were  ejected  more  nearly  normal  to  the 
surface.  Varying  the  isotope  type  of  the  PKA,  atoms  that 
were  known  to  sputter,  and  some  of  the  target  atoms  that 
they  collided  with  resulted  in  differences  in  the  ejection 
angle  normal  of  less  than  0.4°. 


D.   RING  FORMATION 


Results  definitely  show  that  the  sputtering  mechanisms 
do  produce  a  formation  of  rings  of  sputtered  atoms  around 
the  impact  area.   The  results  indicate  that  the  transfer   of 
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momentum  along  the  <100>  and  < 1 0 1 >  directions  in  the  surface 
layer  is  relatively  easy.  The  predominance  of  sputtering  of 
next-nearest  neighbors,  in  the  < 1 0 1 >  direction  is  evident. 
For  example,  at  100-eV,  all  simulations  using  an  ion  of  the 
mass  of  argon  or  greater  sputtered  the  same  atoms,  but  with 
different  freguencies.  The  particular  atoms  that  sputter  do 
not  appear  to  be  dependent  on  the  mass  of  the  incident  ion. 
However,  the  number  of  times  the  atom  sputters  does  appear 
to  be  dependent  on  the  ion  mass.  Ring  and  atom  sputtering 
probabilities  clearly  show  the  sputtered  atoms  from  certain 
site  directions  and  the  distances  of  the  sites  from  the 
center  of  the  PKA. 

A  trend  was  noted:  as  the  energy  and  mass  of  the  ion  is 
increased,  the  sputtered  atoms  come  from  sites  which  appear 
to  spread  outward  from  the  impact  area.  These  higher-energy 
events  give  the  appearance  of  a  "liguid  splash",  with  the 
size  of  the  crater  dependent  on  the  ion  energy  and  mass. 


E.   ION/ATOM  MASS  RATIO 


The  results  of  this  investigation  do  not  support  a 
monatonic  increase  in  the  sputtering  ratio  with  increased 
ion  mass  as  implied  in  figures  of  Sigmund[28]  and  Eer  Nisse 
[27].  Ml  the  plots  showing  the  sputtering  ratio  as  a 
function  of  the  ion  mass,  clearly  indicate  a  maximum  near  an 
ion/atom  mass  ratio  of  one.  At  ion/atom  mass  ratios  greater 
than  one,  the  sputtering  ratio  appears  to  become  asymtoptic 
to  a  lower  value.  This  may  be  an  artifact  of  the  computer 
model.  These  plots  also  highlight  the  containment  problem. 
The  behavior  of  the  sputtering  ratio,  for  larger  ion/atom 
mass  ratios,  is  erratic  because  the  microcrystallite  does 
not  contain  the  event. 
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The  general  shape  of  the  sputtering  ratio  carve,  when 
plotted  as  a  function  of  the  ion/atom  mass  ratio,  is  in 
agreement  with  the  experimental  data  of  Wehner  and  his 
co-workers.  Results  conclusively  show  that  the  ion/atom 
mass  ratio  is  an  important  factor  in  the  determination  of 
the  sputtering  ratio. 


F.   RANDOMLY  SELECTED  IMPACT  POINTS 


The  sputtering  ratio  obtained  using  the  grid  of  impact 
points  is  extremely  close  to  that  obtained  when  randomly 
selected  impact  points  were  used. The  similarity  between  the 
sputtering  summaries  and  the  atom  sputtering  probabilities 
are  guite  evident.  The  guantity  of  data  obtained  with  the 
randomly  selected  impact,  points  was  reguired  to  confirm  that 
light  isotope  enrichment  was  present. 

Based  on  the  comparison  of  results  obtained  from  the 
randomly  selected  impact  points  and  from  the  grid  of  impact 
points,  the  systematic  method  of  "walking"  through  the  grid 
seems  completely  justifisd. 


G.   FUTURE  RESEARCH 


Even   the  more  than  1000  hours  of  computer  time  used  in 

running  these  simulations  did  not   begin   to   encompass   the 

permutations  and  combinations  of  parameters  that  need  to  be 
investigated. 
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The  (110)  and  (111)  orientations  need  to  be  studied  in 
depth.  The  few  sputtering  mechanisms  observed  were  the  same 
as  those  for  the  (100)  orientation.  However,  the  sputtering 
ratio  as  a  function  of  the  ion  mass,  containment 
requirements,  ring  formation,  and  angular  dispersion  of 
isotopes  need  further  investigation.  It  is  highly 
recommended  that  the  effect  of  momentum  transfer  along 
certain  directions  be  investigated. 


63 


APPENDIX  A 


ISOTOPE  GENERATION 


All  atoms  in  the  lattice  vers  one  of  two  types;  type  1  - 
the  heavier,  least  abundant  isotope,  or  type  2  -  the 
lighter,  most  abundant  isotope.  The  average  atomic  weight 
of  copper  is  63.54  amu.  The  isotope  weights  are  62.93  amu 
and  64.93  amu  with  natural  abundances  of  69.1-?  and  30.9% 
respect  ively  . 

The  lattice  was  generated  with  either  the  isotope  masses 
and  their  corresponding  abundance  or  the  weighted  average  of 
the  isotope  masses  with  an  abundance  of  100%.  Initially, 
all  atoms  were  forced  to  the  type  2  atom  in  the  lattice 
generator.  If  the  isotope  lattice  was  desired,  the 
International  Mathematical  and  Statistical  Library  (I5ML) 
routine  GGUB  was  called.  This  routine  generates  a  uniformly 
distributed  set  of  psuedo  -  random  numbers  between  zero  and 
one.  Each  of  these  random  numbers  was  assigned  to  a 
particular  atom  and  was  tested  to  see  if  it  was  less  than  or 
equal  to  0.309.  If  the  number  was  less  than  this  test  value 
(the  abundance  of  the  least  abundant  isotope) ,  the  atom  was 
changed  to  a  type  1  atom.  A  seed  for  the  random  number 
generator,  for  the  particular  size  lattice  to  be  generated, 
was  chosen  by  trial  and  error  which  would  produce  the 
correct  number  of  atoms  above  and  below  the  test  value  to 
assure  the  correct  isotope  ratio. 
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A  feature  was  designed  to  allow  any  number  of  the  atom 
masses  (types)  to  be  selectively  changed  after  the  lattice 
had  been  generated.  This  '  feature  was  incorporated  to 
investigate  the  changes  and  effects  of  picking  heavy/light, 
heavy/heavy,  or  light/light  collision  atoms.  Phis  would 
facilitate  the  dynamic  investigation  of  the  "reflective 
collision"  theory  proposed  by  Winters  and  Sigmund  [10]. 

The  dynamics  of  the  simulation  remained  unchanged  except 
that  all  calculations  that  were  dependent  on  target  atom 
mass  were  modified  to  use  the  correct  mass. 

The  simulation  was  nDw  capable  of  generating  a  lattice 
with  the  correct  isotope  ratio  and  masses  with  randomly 
selected  sites. 
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APPENDIX  B 


LATTICE  FOLDING 


As  mentioned  in  the  simulation  model  development,  a 
lattice  that  was  symmetrical  about  the  folding  axes  was 
required.  Once  this  was  achieved  by  the  lattice  generators, 
it  was  necessary  to  fold  the  lattice  in  a  manner  such  that, 
after  the  folds,  the  reduced  representative  impact  araa  was 
identical  to  the  complete  representative  impact  area.  After 
these  folds'  were  performed  correctly  and  in  the  right 
sequence,  the  results  were  identical  to  that  obtained  if  the 
complete  representative  impact  area  had  been  bombarded 
systematical ly. 

After  each  fold,  the  sputtering  count  for  all  atoms  that 
were  symmetric  about  the  folding  axis  were  summed.  The 
summing  process  started  with  the  first  atom  in  a  layer  to 
the  highest  numbered  atom  in  that  layer  After  this  process, 
the  sputtering  count  for  each  atom  was  twice  what  the  sum 
was  since  the  sputs  were  summed  each  time  a  symmetric  atom 
was  encountered.  This  doubling  effect  of  each  fold  was 
corrected  by  scaling  the  sputtering  count  of  each  atom  after 
all  the  folds  in  a  layer  were  completed.  All  necessary 
folds  were  then  performed  on  each  successive  layer. 
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A.    (100)  ORIENTATION  FOLDING 


Figure  12  shows  the  (100)  orientation  reduced 
representative  impact  area  folding  geometry.  The  reduced  36 
impact  point  area,  when  folded,  represents  the  complete  298 
impact  point  area. 

Since  the  (100)  orientation  impact  area  is  a  45°-45° 
triangle,  three  folds  were  required  to  complete  the  folding 
process.  The  first  fold  was  along  the  diangon  The  first 
fold  was  along  the  diagonal  axis,  the  second  along  the 
horizontal  axis,  and  the  third  along  the  vertical  axis.  The 
diagonal  fold  process  for  all  layers  was  the  same.  However, 
the  horizontal  and  vertical  folds  were  modified  depending  on 
which  layer  was  being  folded. 

After  all  the  folds  in  a  layer  for  the  (100)  orientation 
were  complete,  the  sputtering  count  for  each  atom  was 
divided  by  eight  to  correctly  scale  the  results. 


B.    (110)  ORIENTATION  FOLDING 


Figure  12  also  shows  the  (110)  orientation  reduced 
representative  impact  area  folding  geometry.  The  reduced  36 
impact  point  area,  when  folded,  represents  the  complete  144 
impact  point  area. 

Since  the  (110)  orientation  reduced  impact  area  is  a 
rectangle,  two  folds  were  required  to  complete  the  folding 
process.  The  first  fold  was  along  the  horizontal  axis  and 
the   second   fold  was  along  the  vertical  axis.   The  diagonal 
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fold  was  not  necessary.   The  horizontal  and   vertical   folds 
(exactly   the   same   as   for   the   (100)   orientation)   were 
modified  depending  on  which  layer  was  being  folded. 

After  all  the  folds  in  a  layer  for  the  (110)  orientation 
were  complete,  the  sputtering  count  for  each  atom  was 
divided  by  four  to  correctly  scale  the  results. 


C.    (111)  ORIENTATION  FOLDING 


The  (111)  orientation  was  by  far  the  most  difficult  to 
determine.  Figure  13  shows  the  (111)  orientation  reduced 
representative  impact  area  folding  geometry.  The  reduced  21 
impact  point  area,  when  folded,  represents  the  complete  126 
impact  point  area. 

Since  the  (111)  orientation  impact  area  is  a  30°-60° 
triangle  three  folds  along  the  diagonals  resulted  in  the 
overlap  shown  in  Figure  13.  However,  if  the  set  of  three 
folding  sequences  were  summed  the  results  could  be  scaled  to 
represent  the  complete  impact  area.  The  number  one  fold  was 
along  the  positive  slope  diagonal  axis,  the  number  two  fold 
was  along  the  negative  slope  diagonal  axis,  and  the  number 
three  fold  was  along  the  vertical  axis. 

Before  each  sequence  of  folds,  the  sputtering  count  for 
each  atom  was  reinitialized  to  the  original  value.  The 
first  sequence  (1,2,3)  was  performed  and  upon  completion  the 
sputtering  count  for  each  atom  was  divided  by  eight  to 
correctly  scale  the  results.  The  second  sequence  (1,2,1) 
was  performed  and  again  the  sputtering  count  for  each  atom 
was  divided  by  eight  for  scaling.  The  third  sequence 
(1,3,2)  was  performed  and  again  the  sputtering  count  for 
each  atom  was  divided  by  eight  for  scaling. 
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Finally,  the  results  obtained   from   each   sequence  was 

summed   and   the   final   sputtering   count  for  each  atom  was 

divided  by  four  to  correctly  scale  the  results   due   to  the 
summing  process. 

Since   the   geometry   of   the   (111)    orientation  was 

complicated  and  not  fully  developed,  the  folding  process  was 
performed  only  for  the  first  layer. 
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APPENDIX  C 


GLOSSARY 


AEUN1,  ABUN2 

The    abundanca   of    the    target's    isotopes. 

AC 

A   distance   measurement   used   in   impact    point 
generation. 

ALFA,  ALFA2 

The  cosine  of  the  angle  between  vectors   R1,R2, 
squared. 

ALPHA 

Input  Morse  potential  parameter. 

ANGLE 

The   angle   between   the  direction  of  motion  of  an 
ejected  atom  and  the  surface  normal  in  degrees. 

AVE 

The   sputtering   ratio. 

AIT,    AYT,    AZT 

Components  of  the  random  thermal  displacement. 

BENRGY 

The  lattice  binding  energy  in  ev. 
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BIND 

Negative  of  the  total  potential  energy  (TPOT)  at 
time  zero.  Retained  for  energy  conservation 
checks. 

BMAS 

The  mass  of  the  bullet  in  amu. 

BQZ 

A  fee  crystal  constant  -  (1/R00T (3) )  . 

BOL 

The  bullet  heading  Hollarith. 

BXr  BZ 

Ensealed  x,z  coordinates  of  the  impact  point. 

BXMAX,  BXflIN 

X  direction  boundaries  of  the  target  impact  area. 

BZMAX,  BZMIN 

Z  direction  boundaries  of  the  target  impact  area. 

CELS 

A  frictional  force  multiplier  when  frictional 
forces  are  included  in  the  program. 

CGB1,  CGB2 

Morse   potential    parameters. 

CGD1,    CGD2 

Morse    potential    parameters. 

CGF1,    CGF2 

Morse  force  function  parameters. 
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CFO,    CF1,    CF2 

Force  parameters  of  the  cubic  fit  between  Jlorse 
and   Born-Mayer    functions. 

CO,    CORF 

"Phonon"  displacement  correlation  parameter  used 
as  an  option  in  WARM. 

CORN 

1  -  CORF  . 

COX,  COY,  COZ 

Directional  cosines  of  the  incoming  ion's  velocity 
vector,    relative    to   the    surfaca    normal. 

COXI,  COYI,  COZI 

Negative  values  of  COX,  COY,  COZ. 

CPO,  CP1,  CP2,  CP3 

Potential  parameters  of  cubic  fit  between  Morse 
and   Born-Mayer    functions. 

CVBD 

A    ratio    used    to    avoid    repeated    division. 

CVD 

Distance  conversion  factor  -  angstroms  to  meters. 

CVE 

Energy  conversion  factor  -  ev  to  joules. 

CVED 

A  ratio  used  to  avoid  repeated  division,  CVE/CVD, 
used  to  convert  exponential  potential  functions  to 
forces . 
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CVM 

Mass  conversion  factor  -  amu  to  kg. 

CVR 

Conversion  factor  -  lattice  units  to  angstroms. 

DBXMAX,  DBZMAX 

Upper  boundary  limits  for  the  target  impact  area. 

DCON 

Input  Morse  potential  parameter. 

DFF 

The  distance  closer  than  ROE  that  an   atom   is   to 
the  bullet.  Dsed  when  forces  are  turned  on  or  off. 

DIST 

The  distance  between  any  two  atoms. 

DNN(I) 

Distance  to  the  Ith  ring  in  lattice  units. 

DRX,  DRY,  DRZ 

Components  of  DIST. 

DT 

Length  of  each  timestep  in  seconds. 

DTI 

Number  of  lattice  units  the   most   energetic   atom 
can  move  in  one  timestep. 

DTOD 

DT  over  distance  -  a  ratio  used  to  avoid   repeated 
division . 
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DT0M1,  DT0M2 

DT  over  target  isotope  mass   -   a   ratio   used   to 
avoid  repeated  division. 

DTOMB 

DT   over    bullet      mass      -      a      ratio      used      to      avoid 
repeated    division. 

DX(I)  ,    DY(I)  ,    DZ(I) 

Change  in  the  position  of  the  Ith  atom. 

EMAX 

The   velocity  squared  of  the  fastest  atom.  Used  to 
make  ETEST. 

ETEST 

The  kinetic  energy  of  the  most  energetic  atom. 

ETHR 

Energy    threshold   to    escape   the    surface. 

EV 

Primary  energy  of  the  bullet  in  ev. 

EVE 

Primary    energy   of    the   bullet    in    kev.      R       indicates 
a   read-in    and   saved   value. 

EXA,    EXB 

Input  Born-Mayer  potential  function  for  the  targst 
atom  -  target  atom  interaction. 

FA 

The  temporary   component   force   increment   on   an 
atom  . 
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FAC 


FINISH 


FLAT 


FM 


FMAX 


FORCE 


FPTC 


FRC 


FRAN 


FRANR 


The  decrement  size  used  in  SETUP  in   locating  the 
ion  impact  point. 


The  elapsed  program  running  time  in  seconds. 


A  lattice  unit  in  angstroms. 


A  small  number  used  in  checking   potential   energy 
zero  point. 


The   maximum   freguency   of   occurrence    in   the 
histogram  intervals. 


A  numerical  value  of   the   force   function   for   a 
specific  pair  of  atoms. 


The  corrective  force  at  ROE. 


Numerical   value   of   the  target  force  function  a 
ROE. 


A   real  variable  for  number  of  passes  cooiple-ed  at 
termination . 


The  reciprocal  of  FRAN. 
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FREQ1,  FREQ2 

Number   of   events   in   each    interval    of   the 
histogram. 

PSCAL 

A  histogram  scale  factor. 

FTDI 

A    parameter    used    to   determine    DT    by    maximum    energy 
method. 

FXA 

4 

A    Born-Mayer    force    function    parameter. 

FX(I)  ,    FI(I)  ,    FZ(I) 

Components  of  the  total  force  on  the  Ith  atom. 

HBHAS 

Half   Bullet  Mass  -  a  ratio  used  to  avoid  repeated 
division . 

HDTOD 

Half   DT   over   distance   -   a  ratio  used  to  avoid 
repeated  division. 

HDTOM 

Half   DT  over  most  abundant  isotope  mass  -  a  ratio 
used  to  avoid  repeated  division. 

HDTOMB 

Half   DT   over  bullet  mass  -  a  ratio  used  to  avoid 
repeated  division. 

HTMA$1,  HTMAS2 

Half   target   isotope  mass  -  a  ratio  usad  to  avoid 
repeated  division. 
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IB0X1,  IB0X2 

Ejected  atom  angle  tally  keeper. 

ICODNT 

Flag  used  in  sequencing. 


IDEEP 


IET 


IFFM 


IFLAG 


IFM 


IFT 


IH1 


IH2 


The    first    frozen    layer.    in   the    target 
microcrystal. 


Elapsed   time  in   2.6  microsecond  blocks.  Used  to 
determine  actual  elapsed  time. 


The   atom   number   of  any  target  atom  that  has  its 
mass  constrained  due  its  isotope  ratio. 


Flag  used  in  sequencing. 


The  isotope  type  (heavy  or  light)  of  any  target 
atom  that  has  its  mass  constrained  due  to  its 
isotope  ratio. 


A   number   used   as   a   constant   when   generating 
pseudo-random  numbers. 


Alphanumeric  array  of  the  program  heading 


Alphanumeric  array  of  the  Morse  parameters 
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IHB(1),  IHB(2),  IHB(3),  IHB  (4)  ,  IHB  (5) ,  IHB  (6) 

Alphanumeric  array  describing  the  bullet. 

IHBPOT 

The  type  of  ion-atom  potential. 

IHBULL 

The  bullet  elemental  symbol  or  name. 

IHS 

Alphanumeric  array  for  the  type  and  orientation  of 
the  crystal. 

IHT 

Alphanumeric  array  describing  the  target. 

IHT(1),  IHT  (2)  ,  IHT  (3),  IHT  (4)  ,  IHT  (5)  ,  IHT  (6) 

Alpha  numeric  array  describing  the  target. 

IHTARG 

The  target  elemental  symbol  or  name. 

IHTPOT 

The  type  of  target  atom  -  target  atom  potential. 

ILATT 

A    flag   used    in   the    lattice    generators. 

ILAYER 

A   variable   used  to  identify  which  layer  is  being 
folded . 

ILL 

The  number  of  the  last  atom  in  the  first  layer. 
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IMF 


IN 


INDEX 


INT 


IOR 


IOUT 


IPASS 


IPNUM 


IPNOMS 


IPROG 


IRAN 


The  mass  of  each   target   atom   with   the   correct 
isotope  ratio. 


The  number  of  intervals  in  the  histograms. 


A  variable  used  in  determining  phase  in  a  dynamics 
cycle  step. 


The  width  of  the  intervals  in  the  histograms 


The  target  plane  of  incidence. 


Flag  used  in  sequencing. 


The  number  of  passes  (shots)  completed 


The  impact  point  number. 


Scaled  impact  point  number. 


Input  variable  used  to  skip  attractive  potential 


Random   number   seed    used    in    NORMAL    subroutine    which 
generates    pseudo-random    numbers. 
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IRANGE 

The  upper  limit  in  the  histograms. 

IREAD 

The  number  of  atoms  which  have  the  isotope  type 
constrained  externally. 

ISEEDA,  ISEEDB,  ISEEDC 

Seeds  used  in  the  uniform  random  number  generator 
which  determines  the  isotope  ratio  and 
distribution . 

ISET 

Flag  used  in  sequencing. 

ISHOT 

A    parameter    used   to   shut   down   the    program. 

ISNDM 

The  atom  number  of  a  sputtered  atom. 

ISOGEN 

A  flag  used  as  an  option  to  generate  the  target 
atcm    masses    with    isotope    or    single    element    values. 

ISPX,    ISPZ 

Initial  start  point  x,z  coordinates. 

IT 

Unsealed  x  coordinate  used  in  lattice  generation. 

ITEMP 

The  lattice  temperature  in  degrees  Kelvin. 

ITIHE 

Flag  used  in  sequencing. 
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ITT 

An   odd-even   integer   used  to  determine  atom  site 
location . 

IX,  II,  IZ 

The  number  of  x,y,z,  planes  in  the  crystal. 

IIP,  IIP,  IZP 

Crystal  dimensions  in  x,y,z. 

JN 

The  numerical  value  of  the  histogram  intervals. 

JODT 

A  histogram  sequencing  variable. 

JSCAL 

A  histogram  scale  factor. 

JT 

Unsealed  y  coordinate  used  in  lattice  generation. 

JTS 

Variable  used  to  establish  a  torn  site. 

JTT 

Variable   used   to    establish   a  torn    site. 

K2S0M 

The   square   of  the  sputtering  count  which  is  used 
in  calculating  the  variance. 

KATOM 

The  number  of  atoms  in  a  ring. 
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KFIRST 

A  sequencing  variable. 

KEING 

The  highest  numbered  ring  in  a  layer. 

KSPDT 

The  number  of  sputs  from  all  atoms  in  a  ring. 

KSUM 

The   total    sputtering   count    before    folding. 

KT 

Unsealed  z  coordinate  used  in  lattice  generation. 

LCUT(I) 

Used  to  identify  an  I th  atom  which  can  be  omitted 
from  the  force  and  potential  calculations  by 
changing  LCUT(I)  from  zero  to  one. 

LD 

The  highest  numbered  atom  in  a  mobile  layer. 

LL 

The  highest  numbered  atom  in  the  entire  crystal. 

LL2,    LL3 

Twice  and  triple  the  value  of  LL. 

LNUM  (I) 

The  number  of  times  the  Ith  atom  sputters  after 
all  shots  in  the  target  impact  area. 

LNUHF(I) 

The  number  of  times  the  Ith  atom  spatters  af*:er 
all  shots  and  the  summary  is  complete  and  folded. 
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LNUMI(I) 

The   number   of   times  the  Ith  atom  sputters  after 
all  shots  but  before  folding  the  summary. 


LS 


LSS 


LSDR 


LTYPE 


MAX 


MCBO 


N1 


NCT 


Sum  of  the  Miller  index  integers. 


Variable  used  to  identify  the  type  of  surface. 


The  lattice  surface  type. 


The  type  of  lattice  (i.e.  fee  )  . 


A  height  scale  factor  used  in  the  histograms. 


The  number  of  simultaneous   equations   solved   by 
CEOSYM. 


NLAY  -  1,  a  sorting  variable. 


Increments  the  timestep  interval. 


NCTT,  ND 

Data  output  interval  in  number  of  timesteps.  Also 
the  timestep  interval  between  potential  energy 
calculat  ions . 
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NF1 

A  sequencing  variable. 

NFIRST 

The  number  of  the  first  atom  in  a  given  layer. 

NLAY 

The  highest  numbered  atom  in  a  given  layer. 

NLAY1 

The  highest  numbered  atom  in  the  first  layer. 

NLINE 

The  line  number  in  the  output  block. 

NOTH 

A    Hollarith   consisting    of   blanks. 

NPAGE 

A    page   numbering    variable. 

NRAN 

Variable  for  number  of  passes  (shots)  completed  at 
termination . 

NES 

A    restart    flag. 

NS 

The   next   timestep  number  at  which  output  will  be 
printed. 
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NSHOTS 

The  number  of  passes  (shots)  to  be  completed  in  a 
program.  Used  to  break  up  an  extremely  large, 
many-shot  program  into  a  few  small  programs  with 
manageable   computer  run   time. 

NSPOT(I) 

The  number  of  atoms  that  sputter  for  the  Ith 
impact  point . 


NSR 


NSS 


NT 


NTT 


PAC 


The  first  timestep  where  output  is  printed 


Flag  to  print  final  output  status  of  current  shot 


The  timestep  number. 


Timestep  number  limit  before  shut  down. 


Parameter  for  bullet  force  function  correction 


PBMAS 

Bullet  mass  in  kg. 

PEXA,  PEXB 

Input  Born-Mayer  potential  function  parameters  for 
the  bullet-target  interaction. 

PFIV 

Point  five  (0.  5)  . 
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PFPTC 

Bullet  force  function  evaluated  at  ROE. 

PPRC 

Bullet-atom      repulsive    force    interaction    evaluated 
at    ROE. 

PFXA 

Bullet  force  function  parameter. 

PKE(I) 

The  kinetic  energy  of  the  Ith  atom. 

PNEG 

Total  negative  potential  energy  in  the  lattice. 

PNOM 

Alphanumeric  variable  for  the  impact  point  number. 

POT 

The   potential  energy   between   two    atoms. 

PPE(I) 

The  potential  energy  of  the  Ith  atom. 

PPOS 

Total  positive  potential  energy  in  the  lattice. 

PPTC 

Bullet  potential  function  evaluated  at  ROE. 

PBOB 

The   probability   that  an  atom  will  sputter  from  a 
given  ring. 
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PTC 

Target   atom    potential   function    evaluated    at    ROE. 

PTE(I) 

The  total  energy  of  the  Ith  atom. 

QM 

A  small  number  used  in  checking  the  kinetic  energy 
zero  point. 

QUIT 

Coutoff  variable  checked  against  the  total 
potential  energy  which  can  be  used  as  a  shut  down 
criterion. 

R1,  R1S 

Distance  from  (0,0)  impact  point  to  the  initial 
bullet  position. 

R1,  R1Y,  R1Z 

The  x,y,z  components  of  R1S. 

R2 

Magnitude  of  vector  from  impact  point  to  the  first 
atom  hit  by  the  bullet. 

R3 

Magnitude  of  vector  from  impact  point  to  bullet 
start   position. 

R1,    R2 

Scalar  product  of  vectors  R1,  R2. 

R2SQ 

R2  squared. 
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RAD 

The  angle  between  the  direction  of  an  ejected  atom 
and  the  surface  normal  in  radians. 

BADDEG 

Conversion  factor  -  radians  to  degrees. 

RAN  (I) 

A  list  of  normally  distributed  random  numbers. 

RANGE 

A  scaling  variable  used  in  the  histogram. 

RBX,  RBZ 

Unsealed   x,z   coordinates   of   the   impact   area 
reference  point. 

RDIST 

The  distance  between  an  atom  and  the  center   point 
of  its  layer. 

RE 

Input  Morse  potential  parameter. 

RINCR 

The  RANGE  increments. 

RMIDPT 

The  midpoint  of  each  interval  in  the  histogram. 

ROE,  R0E2 

Nearest    neighbor    distance    -    squared. 

ROEA 

The    maximum    cutoff    for    the    Born-Mayer    potential. 
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ROEB 

The  minimum  cutoff  for  the  Morse  potential. 

BOEC,  ROEC2 

The  maximum   cutoff   for   the   Morse   potential 
squared. 

BOEM 

Defines   the   region   where   the   repulsive   force 
function  must  be  modified. 

EX  (I)  ,  BY  (I),  HZ  (I) 

X,y,z  coordinates  of  the  Ith  atom  at  any  time. 

BXBL,  BYBL,  BZBL 

Maximum   negative   x,y,z   coordinates   an  atom  may 
move  before  it  is  LCDT. 

BXBH,  BYBM,  BZBM 

Maximum   positive   x,y,z   coordinates   an  atom  may 
move  before  it  is  LCUT. 

BXF,  RYF,  BZF 

Final   x,y,z    coordinates   of   the   Ith   atom    when    it    is 
LCOT. 

BXI(I)  ,    BYI  (I)  ,    BZI(I) 

The   Ith    atom's    initial   x,y,z    coordinates. 

BXK(I)  ,  RYK  (I)  ,  BZK(I) 

Temporary  x,y,z  coordinates  of  the  Ith  atom  luring 
the  force  cycle. 

BXS,  BYS,  BZS 

Scaled  x,y,z  coordinates  of  the  impact  point. 
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SAVE 

One-half    POT. 

SBX,    SBZ 

Scaled    x,z    coordinates    of    the    impact    point. 

SCALE 

A   scaling  variable  for  the  frequency  of  events  in 
the  histogram. 

SCX,  SCY,  SCZ 

The   x,y,z   lattice  scale  factors  used  to  convert 
the  generated  lattice  to  match  a  real  crystal. 

SCXR,  SCZR 

The  inverse  of  SCX  and  SCZ  -  a  ratio  used  to  avoid 
repeated  division. 

SDIST 

The   square  of   the  distance  from  the  center  of  a 
layer  to  an  atom. 

SETIME 

A   library   subroutine   used   to  keep  track  of  run 
time  . 

SLOW 

A   cutoff  velocity  variable  checked  against  a  long 
DT. 

SPDX,  SPDZ 

The   fraction   of   the  target  impact  area  that  the 
impact  points  are  incremented. 
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SPX,  SP2 

The  x,z  components  of  the  distance  from  the  impact 
area  reference  point  to  the  impact  point. 

sscz 

A   z   scale   factor   used  in  the  (111)  orientation 
lattice  generation. 

STD 

The  standard  deviation  of  the  sputtering  ratio. 

STOPER 

This   variable   stops   calculating   and   saves   (SO 
seconds  for  compilation  and  output. 

SOR 

Alphanumeric  variable  describing  the  surface. 

TA,  TAHP 

The  thermal  amplitude  of  the  lattice. 

TAR 

Alphanumeric      variable        describing        the         target 
element. 

TARGET 

Alphanumeric      variable        describing        the         target 
element. 

TAS 

A    scale    factor    -    TAMP/ROOT  (3 .0)  . 

TASF,    TASN 

Random  thermal  amplitude  scale  constants. 
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TE 

The  total  kinetic  and  potential  energy  of  all   the 
atoms. 

TEMP 

The  temperature  of  the  lattice  in  degrees  Kelvin. 

TESTE 

A  termination  test  of  energy. 

TEAC 

A  time   factor   ratio   used   to   determine   DT   by 
maximum  force  methods. 

THERM 

Mean   thermal    energy   of    a      lattice      atom.       Used      in 
shut-down    tests. 

TIME 

The  elapsed  calculating  time  in  seconds. 

TITLS1,  TITLE2 

The  titles  for  the  histogram  outputs. 

TMAS1,  TMAS2 

The  target  isotope  masses  in  amu. 

TPKE 

The  total  kinetic  energy  of  all  the  atoms. 

TPOT 

The  total  potential  energy  of  all  the  atoms. 

TEON 

The  time  the  program  is  to  run  in  minutes. 
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VEL 

The   velocity    calculation    -    SQRT (ev/HBMAS) . 

vss 

A  temporary  storage  variable  for  the  velocity 
componen  ts. 

VTEST 

A  test  velocity  used  to  determine  if  the  atom  has 
sufficient  velocity  (energy)  to  escape  the 
surface. 

VTS 

A  scale  factor  -  0.8*VTEST. 

VX(I)  ,  VI  (I),  VZ(I) 

The    x,y,z    components    of    the    Ith    atom's    velocity. 

VXS,     VYS,    VZS 

The  x,y,z  components  of  the  bullet's  initial 
velocity . 

XDIST 

The  x  component  of  the  distance  from  the  center  of 
a  layer  to  the  atom. 

XMAX,  XMIN 

The  x  planes  to  determine  whether  two  atoms  are 
close  enough  to  interact. 

YLAX(I) 

The  crystal  relaxation  of  the  Ith  layer  in  the  -y 
direction. 
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TMAX,  YMIN 

The  y  planes  to  determine  whether   two   atoms   are 
close  enough  to  interact. 

ZDIST 

The  z  component  of  the  distance  from  the  center  of 
a  layer  to  the  atom. 

ZE 

Absolute  zero  (0.0). 

ZMAX,  ZMIN 

The   z   planes   to  determine  whether  two  atoms  are 
close  enough  to  interact. 

ZSTOP 

Used  to  print  out  a  warning  header  that  a  pass  was 
terminated  before  completion. 
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APPENDIX  D 


COMPUTER  PROGRAM 


C 

c 

c 
c 
c 
c 
c 
c 
c 
c 
c 


GLD3  IS  OLD1  MGDIFIEC  TG  INCREASE  SPEED 

ROLD  IS  GLD3  WITH  'WARMED*  LATTICE 

RING  IS  RPF  MODIFIED  TO  FGLD  LATTICE  AND 

STUDY  RING  FORMAT  I  ON/SPUTTER ING 
RMAS  IS  RING  MODIFIED  TG  STUDY  TARGET  ISOTCPES 
RANG  IS  RMAS  MOD.  TG  STUCY  ANGULAR  CISPERSIGN 


90 
90 
90 
90 
90 


J5 
1C 
2C 

30 

1   c 


C 
D 
C 
C 
C 
C 

c 

C 
X 

c 
c 

CC 
CC 

X  F 
CC 
CC 

X  P 
CC 
CC 
CC 
CC 
CC 
CC 
CC 
CC 
CC 
CG 
CC 
CC 
CC 
CC 

X 
CC 

X  L 
CC 

X  c 
CC 
CC 
CC 
CC 
CC 
CC 
FC 
FC 
FC 
FC 
FC 


MENSICN 

MENSICN 

MENSICN 

MENSICN 

MENSICN 

MENSICN 

MENSICN 

MENSICN 

ITLE2(24) 

MENSICN  H 

MENSICN  I 

MM0N/CCM1 

MMCN/CCM2 

m,pfiv,tp 

MMCN/C0M3 

MMCN/C0M4 

EX,RBZ 

MMCN/CGM4 

MMCN/C0M5 

MMCN/C0M6 

MMCN/CCM7 

MMCN/C0M8 

MMCN/C0M8 

MMCN/CGM1 

MMCN/C0M1 

MMCN/C0M1 

MMON/COMI 

MMCN/C0M1 

MMGN/C0M1 

MMON/COMI 

MM.CN/C0M2 

IHBULL(2) 

MMCN/C0M2 

SLR 

MMCN/C0M2 

F1,CF2,CG 

MMGN/C0M2 

MMCN/C0M2 

MM0N/C0M2 

MMCN/C0M2 

MMCN/C0M3 

MMCN/CCMA 

RMAT  (20A 

PMAT(20A4 

PMAT( 10A4 

PMAT  (2A4 

FMAT  (F5. 


X(600)  ,DY(600) ,DZ(60C) 

XK(600)  ,  RYM6GJ)  ,RZK(6Q0) 

KE(600)  ,PTE(600) 

SPUT(50) ,IPNUMS(50) 

SNUM(20) ,ISEQ(600) 

XF(600),RYF(600),RZF(6CO) 

BOXH20),  IBQX2(20),  RANNUM(250) 

PEQK80)  ,  FREQ2(80)  ,  RANGE  (20  )  ,  T  IT  LE  1  (24  )  , 


CT0M(2) 


FFM(600)  ,IFM(600) 

/RX(60C) ,RY(600) ,RZ( 

/RCE,RCE2,RGEM,AC,PA_ 

fT 
EXA,EXB,FXA,PEXA,PEXB,PFXA,TI- 
IX,IY,IZ,IXP,IYP,SCX,SCY,SCZ, 


/ 

CT 
/E 
/ 


*fcOC) ,LCUT(6 
AC,PPTC,PTC, 


THE 
I 


OC)  ,LL,LD 
PFPTCtFPTC, 

RM 
DEEP, 


A/  LNU 
/FX(60 
/COXI , 
/R1,LS 
/VX(60 
A/VXS( 
1/RXI ( 
3/  RAN 
4/  AXT 
6/BMAS 
6A/3X, 
7/TIME 
8/  ITE 
O/IHK 


M(6J 
OtF 
COYI 
S.SP 
C)  ,V 
600) 
600) 
(180 
(6JJ 
,TMA 
BY,B 
rNT, 
MP,T 
20)  , 


0) ,LNU 
Y(600) 
,CCZI  , 
X,SPZ, 
Y(600) 
,VYS(6 
TRYI(6 
0) »IRA 
) ,AYT( 
SltTMA 
ZtEDIS 
DT,DTI 
AMP, CO 
IHS( 10 


MI(6 

,FZ( 

RXS, 

COX, 

,VZ( 

00)  , 

00)  , 

N,TA 

60J) 

S2,A 

P,CV 

,TEM 

RF 

),IH 


00)  ,L 

600)  , 

RYS,R 

CCY,C 

6QC) 

VZS(6 

PZK6 

SN,TA 

,AZT( 

BUM, 

R 

P,CEL 


NUMF(60C) 
PPE(600  ) 
ZS,FAC 
OZ,YLAX  (2J) , ILL 


00) 
00) 

SF,LL2 
60  D) 
ABUN2, 


,LL3 
I  FNUM 
S,EV,NPAGE,NLINE 
T(6) ,IHB(6)  ,  IHTARG( 21 , 
1/IHBPCT(2),  IHTPCT(2)  ,  SUR ( 2 )  ,  IF2  ( 10 ) , I CR , 


CF2,CP3,CF0, 

CC- 


2/R0EA,RCEB,RCEC,RQEC2,CP0,CPl, 

C1,CGD2,CGB1,CG3  2,CGF1 ,CGF2, I  PR 

3/  SPDX,SPDZ,DBXMAX,DBZMAX,SSCZ 

4/  DNN (100) ,KATOM( 100  )  ,KSPUT ( 100  ),PR08  (  100  ) 

4A/  KSUM,KRING,RDIST(6CC) 

5/IMF(60G) , ISOGEN 

0/  ISEEDA,  ISEEDB,  ISEEDC 

/  A(4, 5) ,MCRG 

4) 

,I3,2A4,I2,3I3,2F3.0,2F5.2,I2) 
,A4,F8.3,2A4,2F10.5,F8.5,3F7.4, 12) 
2,15) 
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904C  FCRMAT(10A4,3F5.2,3F8.5,I1) 

905C  FCFMAT  (2X,  I3,15(  13,  1X,I1,2X) ) 

90  80  FCFMAT(5X,3F5.3,4I5,2F8.6,2X,A4,3F6.4) 

91C5  FCRMATt 11X,F5.0) 

96CC  FCFMATt//) 

96C1  FCFMATt  1H1) 

96 2C    FGBMAT(1X,3(I5,4X,I1,AX,3F6.2,8X) ) 

96  3  0    FCFMAT(105X,AHPAGE,I3,/,1H1 ) 

964C    FGFMAT(«     BINDING    ENERGY    =',F6.2,',    TEST    SPEED    =  «, 

X    F10.2,'    M/SEC  ,  '  , SAFETY    FACTOR    SPEED    =  ',F1C2, 

X     •     M/SEC  ,/,lHl) 
96  80    FCFMAT (2X,  '  ATOM'  ,  2X  ,  »  LCUT  '  ,9X  ,  •  VX  '  t9X  ,  •  V  Y  •  ,  9X  ,  '  VZ  •  , 

X    8Xf  •RXFS8X,  •RYF»,8X,,RZF*  ,8X,  fPPES8X,  'PKE'tSX, 

X  'FTE'/) 

9681  FCRMATt IX, 'ATOM    DX»  ,6X, • DY' , 6X  ,  •  DZ ' ,7X  ,  • VX »  ,7X ,  • VY •  , 
X   7X,'VZ',7X,'KE',8X,'PE',8X,'TE',7X,'RX',6X,'RY',6X, 
X   aRZ''f4Xf  »LCUT    ISEQ',/) 

9682  FCRMAT(I5,3F8.3,3F9.1,3F10.4,3F8.2,I4,3X,I4) 
96  85  FCFMATt2X,I3,3X,I3,5X,3tF8.1,3X),6(F8.3,3X) ) 

96SC  FCFMAT(/4X,F10.3, 'EV,  TCTAL  KINETIC  ENERGY  ,', F10 .3 , 

X   'EV,  TOTAL  POTENTIAL  ENERGY, • ,F10. 3,' EV,  TCTAL', 

X   '  ENEFGY,  •  ,/9X,  'TCTAL  POSITIVE  POTENTIAL  ENERGY  IS', 

X   F10.3,'TCTAL  NEGATIVE  POTENTIAL  ENERGY  IS', FLO. 3) 
9695  FCFMATt3X, ' ATOM  &  POTENTIAL  SUR F  •  ,5X , ' VX  •  ,  7X ,  • VY •  , 7X , 

X'\zZ'  ,  10X,  'PKE'  ,6X,  'PPE'  ,6X,  'PTE    KEV   CCX    COY',//) 
9701  FOPMAT<«  IFT  =',110,',    IRAN  =»,I10,»,     LL=',I5,//) 
98  0  0  FCFMATt  IX  ,  2  15  ,2A4  ,  14 , A4  ,  IX , 3F9 . 1 , 2F9 .4, F5 . 1 ,1X,2F6.2) 
98C5  FCRMATt 5X, • SPUTTERED  ATCM  /  ANGLE  WITH  NORMAL'/) 
9810  FCFMAT(//,'  AT', F7.lt'  EV,  THE  FOLLOWING  ATCMS  WERE', 

X  •  SPUTTEREC:'/,20I4,//) 
9815  FCRMATt  1X,5(5X, I3,3X,F7.4) ) 

982C  FCFMAT(  •  THERE  ARE  NO  SPUTS  AT  ',F8.1,'EV') 
984C  FCFMAT(25X,  'SUMMARY  OF  TARGET  IMPACT  AREA'/,25X, 

X   'FOLDED  TO  REPRESENT  TOTAL  COVERAGE'/) 
9845  FCFMAT  (IX,  'ATOM  NUMBER  /  TIMES  SPUTTERED  ') 
985C  FCRMATt///,'  SPUTTERING  RATIO  WAS',F6.3, 

X   ',   THE  VARIANCE  WAS',F8.3,',  FCR  ',14, 

X   '  IMPACT  POINTS,  THE  STD  DEV  CF  THE  MEAN  WAS', 

X    F6.3,//) 
9860  FCPMAT(2X,  ' STATISTICAL  SUMMARY  FCR ' , 14 ,  '  I M P ACT  POINTS' 

X  ,'  AT  ',F5.1,  »KEV  ,///) 
9865  FCRMATt'  NUMSER  OF  TIMES  EACH  ATCM  SPUTTERED'/) 
987C  FCFMAT  tlOt  16, 14)  ) 
9875  FCPMAT(//,»  BINDING  ENERGY  WAS»,F6.2,«  EV,  TEST  SPEED' 

X  ,»  WAS',F10.2,'  M/SEC,  SPUTTERED  ATOM  COLNT  IS', 15) 
9880  FCFMAT( 1X,3(»  RING  DISTANCE  ATCMS  SPUTS  PRGB  •)) 
9885  FCFMAT( IX, 3(  '  NN(  •  ,  I  2  ,  '  )  '  , 2X , F6 .3  ,3X , 12 , 5X ,  13  ,  3X , 

X  F6.3,1X)) 
989C  FCFNAT(25X,  'PROBAEILITY  THAT  A  RING  OF  ATCMS  SPUTTERS' 

X  //  ) 
9895  FCFMAT(25X, 'SPUTTERED  ATOMS  ANGULAR  DENSITY'/) 
99CC  FCRMAT(2X,'  INCREMENTS  ARE  IN',F4.2,'  DEGREE  INTERVALS' 

X  /) 
9905  FCFMAT(2X, 'ANGLE  (DEGREES)  /  NUMEER  OF  ATCMS'/) 
9910  FCRMAT  (2X  ,  19( I  3 , 3X ) , 13 ) 
9915  FCRMAT  (3X,I2,2(      IX, FA. 2)) 
9920  FCFMAT('  00-05  05-10  10-15  15-20  20-25  25-3C  30-35', 

X  •  35-40  40-45  45-50  50-55  55-60  60-65  65-7C  70-75', 

X  •  75-80  80-85  85-90' /) 
9925  FCFMAT(«  CO-10  10-20  20-30  30-40  40-50  5C-6C  60-70', 

X  •  70-8C  80-90'/) 
993C  FCPMAT('  NUMBER  OF  ATOMS  CHANGED  WAS', 14/,'  ATOMS', 

X  ■  CHANGEC  AND  ISCTOPE  TYPE  ARE') 
9935  FCRMAT  ( IX ,  13 ( I  3 , 2X , I  1 , 3X )  )  rn 

9940  FCRMAT  ( IX  ,  '  SCX  =  '  , F7 .4  ,  '    SCY=',F7.4,»    SCZ=»,F7.4) 
9945  FCRMAT(3( 110, 10X) ,15) 

995C  FCFMAT(«  ISEED  =  ',110)  _   ___  „  TTt 

9955  FCFMATt IX, •***  THIS  SHOT  WAS  TERMINATED  BEFCRE  IT', 

X    •  WAS  COMPLETED'/, 4X , '  ZSTGP=1   THE  CURRENT', 

X    '  TIMESTEP  EXCEEDED  THE  LIM  IT  '  / , 14X , • ANY  SPUTS', 

X    '  ARE  REMOVED  FROM  THE  SUMMARY') 
996C  FCFMAT( IX, • ***  THIS  SHOT  WAS  TERMINATED  3EFCRE  IT«, 
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9965 
9970 
9975 


6 
I 
L 

F 

P 

C 
<; 

W 

c 
I 
I 

1  c 

D 
I 

L 
R 
R 
R 
V 
V 
V 
P 
R 
R 
F 
F 
F 
F 
P 
P 
R 
R 

2  R 

R 

X 
R 
R 

X 
R 
5    R 

X 
R 
R 

P 
I 
I 
7 
A 
A 
2  C 
T 
C 
I 


•     WA 

14X, 

CPMAK 

CRMATC 

CFMAK 
ALL  CV 
ALL  SE 
EAD  (5 
EAC  (5 
EAC<5, 
T0PER= 
INISH= 
LATT=0 
CFC  =  4 
CTT=10 
ESTE=2 
EIND 
CVB  = 
ENRGY= 
SHLT=0 
C  =  0 

M1.0E 
FIV=0. 
MO. 01 
LCU=1. 
PITE(6 
C    1    1  = 

ecxid 

ECX2(I 
CMINU 
C  2  1  = 
NUMd  ) 
CLT(I) 
XK(I)= 
YK(I )= 
ZK(I)= 
X(I)=0 
(I)=0 
(I)=0 
(I)=0 
(I)  =  J 
(I)=0 
(I  )=0 
(I  )=0 
(I)  =  0 
TE(I)= 
KE(I )= 
PE(I)= 
XI(I)= 
VKI  )  = 
ZICI  )  = 
READ  { 
EAD(5, 
LTYPE 
EAD  (5 
EAC(5, 
,ABUM 
EAD  (  5 
£AG(5, 
,7A,CC 
EAC(5, 
EAD  (5 
EAD  (5 
F(EVR. 
F(  ISQG 
PAS2=6 
EUM  =  C 
ELN2=1 
CNTINU 
AMP=TA 
CPF=CC 
PASS=1 


S  CO 
•  AN 
4X,» 
4X,« 

ax, ' 

FLOW 
TIME 
,90  0 
,90C 
91C5 
TRUN 
0.0 


PPLETED'/,14X, •  ISHUT=1 
Y  SPUTS  ARE  REMOVED  FRCM 

IMPACT  POINT  • , 15) 
BX=«  ,F6.4,  »BZ=»  ,F6.4) 
NOT  IN  TARGET  APEA' ,/) 


5)  TITLE1 
5)  TITLE2 
)  TRUN 
*60. 0-60.0 


TIME  RAN  GUT'/, 
THE  SUMMARY* ) 


.0 

=  0.0 

0.774E-9 
2.0 


-1J 

c 


OE- 
,96 
1,2 
)=0 
)  =  0 
E 

1,6 
=  0 
=  0 
CO 
CO 
CO 
.0 
.0 
.0 
.0 
.0 
.0 
.0 
.0 
.0 
0.0 
CO 
CO 
Q.O 
CO 
CO 
5, 
902 


14 

01) 

0 


oc 


9C10)  I  HI 

0  )IHS, IQR,SUR,L5S,IX, IY,IZ,RBX,RBZ,DTI,RCE2, 


,903 
9030 
,A8U 
,904 
9080 
,ETH 
<5945 
,991 
,905 
LE.O 
EN.N 
3.54 
.000 
.000 
E 


C)  IHBULL,BUL,BMAS, IFEPGT,PEXA,P£X6,XXXX 

)IHTARG,TAR,TMAS1,IHTPGT,EXA,EXB,FLAT,TMAS2 

N2, ISOGEN 

G)  IH2  ,ROEA, ROES,  ROEC, DC CN, ALPHA ,RE, IPRCG 

CNSR,  NCI  DEEP  ,CO>,CGY,  ITE^P 


)EVR,SPX,SPZ,NTT 

R 

ISEEDA,ISEED3, 
IRANGEf RMIDPT 
IREAD,(  IFFM(J 
GO  TO  9999 
GO  TO  33 


) 

5) 

C) 

.0) 

E.O) 


ISEEDCNSHGTS 

,  R  I N  C  R 

)  ,IFW(J)  ,J  =  1, IREAD) 
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NCTT=ND 
IF(ILATT.NE.O) 


GO    TO    10 


MAKE    HEADING    HCLLARITH 


IFEQ) 
IFB(2) 
IFE(3) 
IHE(4) 
11-6(5) 
Il-e(6) 
IHT(L) 
IFT(2) 
IFT<3) 
IFT<4) 
IFT(5) 
IFT<6) 


IHBU 
IHBU 
IHTA 
IHTA 
IHBP 
IHBP 
IHB( 
IHB( 
IH8( 
IHB( 
IHTP 
IHTP 


LL(1) 

LL(2) 

RG(1) 

PG(2) 

CT(1) 

CT(2) 

2) 

4) 

3) 

4) 

CT(1) 

CT(2) 


MAKE    CONSTANTS 


RACCEG=57. 2957795 
RCE=SQRT(RCE2) 
EV  =  EVR*1C0C0 
NS=NSR 

FAC=0.1 
PCEP=RCE-CTI 
CVR=0.5*FLAT 
CVE=1.6C21E-19 
CVN=1.6604c-27 
CVC=CVR*1.0E-ld 
CVEC=CVE/CVD 
RXBL=-RCE 
RYBL=-J.8 
RZBL=-RCE 
IXP=(  IX  +  D/2 
IVF  =  (  IY  +  D/2 
IZP=(IZ+l)/2 
PEMAS=BMAS*CVM 
HEMAS=0.5*PBMAS/CVE 
PTMAS1=TMAS1*CVM 
PTNAS2=TMAS2*CVM 
HTMAS1=    0.5*PTMAS1/CVE 
HTMAS2=    0.5*PTMAS2/CVE 
VTEST=-SQRT(2.0*ETHR*CVE/PTMAS2) 
VTS=0.8*VTEST 

MAKE    REPULSIVE    POTENTIALS 

FXA=ALOG(-EXB*CVED)+EXA 

PFXA=ALCG(-PEXB*CVED)+PEXA 

AC=ALOG(CVEC)+EXA 

PAC=AC-EXA+FEXA 

PTC=EXP(EXA+EXB*RGE) 

FRC=EXP (FXA+EXB*RCE) 

FPTC=EXF(AC+EXB*RCE) 

PFPTC=EXP(PAC+PEXB*ROE) 

PFRC=EXP(PFXA+PEXE*ROE) 

PPTC=EXP(PEXA+PEXe*ROE ) 

GENERATE  THE  LATTICE 

IF  (I0R.EG.100)  GC  TO  6 
IF  (I0R.EQ.110)  GC  TO  7 

IF  (IOR.EQ.lll)  GO  TO  8 

CALL  KF100 

GC  TO  9 

CALL  KF110 

GC  TO  9 

CALL  KF1L1 

GC  TO  9 

CCNTINUE 

IF(  I0R.EQ.1C0)  THETA  =  1.0000 


98 


12 


65 


10 


25 


30 


F(IGR. 

F(IOR. 

F(  ISOG 

C    12    I 

=IFFM( 

VF(J)= 

CNTINU 

C    65    I 

=  1  +  1 

=  1+2 

RITE(6 

RXI(K) 

RITE     ( 

PITE 

RITE 

PITE 

RITE 

L2=2*L 

L2=3*L 

1S-R1 

LATT=1 

FAGE=1 

>=RBX* 

z=Rez* 

ex=bx 

ez=bz 

XBM=  (I 
YBP.=  (I 
Z8M-U 

CNTINU 
CZ-1.0 

CZ=ABS 
CZ=SQR 
CXI=-C 
CVI=-C 
CZI=-C 
FT=ITI 
RAN=2* 
F(  IRAN 
XMIN=B 
XMAX=e 
ZNIN=B 
ZMAX=E 
F(  IOR. 
F(SFX. 
ZMAX=B 
ZMIN=B 
F(ICR. 
F4SPX. 
ZMAX=8 

zmin=b 

CNTINU 
X=(RBX 
\=(RYI 
Z=(RBZ 
CXR=  1 
CZR=1. 
F(IOR. 


EQ.113)     THETA=C.OOOO 

EQ.lll)    THETA=(SQRT(  1 

EN. ECO)     GC    TO    12 

=1,IR5AD 

I) 

IFM(  I) 

E 

=2, LI, 3 


5)/3.0)/(SQRT(2.0)/2.0) 


,9620)1 

,RYI (K) 

6,9600) 

6,9930) 

6,9925) 

6,9600) 

6,9940) 

L 

L 


SCX 
SCZ 


INF(I)fRXKI)fRYKI)  tRZKDt  Kf  IMF(K) 

RZI(K),J,IMF(J),RXI(J),RYI(J),RZI(J) 

IREAD 

(  IFFM( J) ,IFM( J)  ,J  =  1, IREAD) 

SCX,SCY,SCZ 


X-l 
Y-l 
Z-l 

E 

-CC 

(CO 

T(C 

CX 

CY 

CZ 

NE( 

IFT 

.LT 

X 

X+D 

Z 

Z+D 

NE. 

EQ. 

ZMA 

ZMI 

NE. 

EQ. 

ZMA 

ZMI 

c 


)*SCX+RCE 
)*SCY+RCE 
)*SCZ+RCE 

X*COX-CCY*COY 

Z) 

CZ) 


II) 
+  1 

.0)     IRAN=-IRAN 

BXMAX 

8ZMAX 

100)     GO    TO    25 

0.0)    GO    TO    25 

X-10.0*SPX*SPCZ 

N+10.0*SPX*SPCZ 

111)    GO    TO    30 

0.0)    GO    TO    30 

X-5.0*SFX*SPDZ 

N+5.0*SFX*SPCZ 


FT=ITI 
RAN=2* 
FURAN 
RITE  ( 
PITE(6 
RITE(6 
F(TAMP 
£S=TAM 
CRN=1. 
ASN=TA 
ASF=TA 


+SPX)*SCX 

(4)  +  C001)*SCY 

+SPZ)*SCZ 

.O/SCX 

C/SCZ 

EQ.lll)     SCZR=1.0/SSCZ 

PREPARE  TC  WARM  LATTICE 

NE(II) 

IFT  +  1 

.LT.O) 

6,9600) 

,9701) 

,9640) 

.EQ. CO)  GC 

P/SQRT(3.0) 

C-CCPF 

S*CCRN 

S*CCRF 


IRAN=-IRAN 

IFT, IRAN, LL 
ETHR,VTEST,VTS 
TO  40 


99 


c 
c 

c 


40 


50 


60 


LCC 


2C0 


21C 


EACH  REPLICA  STARTS  HERE 

CCNTINUE 

IF(IPASS.GT.NSHQTS)  GO  TO  10JD 

IF(TAMP.GT.C.O)  CALL  WARM 

CCNTINUE 

TiyE=0.C 

NS=NSR 

NT  =  C 

NSS  =  0 

NCT=0 

INCEX=0 

ZSTCP=0.0 

SPX=(BX-S8X  )*SCXR 

SPZ=(EZ-SBZ)*SCZR 

ISPX=(SPX+0.001)*100.0 

ISPZ=(SPZ+0.001)*1CO.O 

IFNUN=1C0*ISPX+ISPZ 

IFNUMS(  IPASS)=IPNUM 

R1=R1S 


RECCNSTRUCT  LATTICE 


IS 
CC 
LC 
IS 
RX 
RY 
RZ 
FT 
PP 
VX 
VY 
VZ 
CC 
IF 
RX 
RY 
RZ 
Rl 
CA 
LC 
RX 
RY 
FZ 
PF 
VE 
VX 
VY 
VZ 
FC 
CT 


EG(1)=0 

60    I=2tLL 
LT( I) =0 
EC( I)=0 
(I)=RXI(I  ) 
(I)=RYI(  I  ) 
(  I)=RZI(  I  ) 
E  ( I )  =  0  .  0 
E(I)=C.O 
(I)     =    0.0 
(I)     =    0.0 
(I)    =    0.0 
NTINUE 
(TAMP  .NE 
S  =  EX 
S  =  BY 
S  =  EZ 
=R1+FAC 
LL    SETUP 
LT(1)=0 
I  (1)=RX(1  ) 
I (1)=RY(  1) 
I(1)=RZ(1) 
E(1)=Q. J 

L=SQRT(EV/HBMAS  ) 
(1)=VEL*CCX 
(1)=VEL*CCY 
<1)=VEL*CCZ 
T1=DTI*CVC 
=FDTI/VEL 


CO)     CALL    WARM 


BEGIN    TIMESTEP    CYCLE 


DTCD=DT/CVD 

FCTCD=0.5*QTGD 

CTCM1=DT/PTMAS1 

CTCM2=CT/PT^AS2 

HCTCM(  1)=0.5*DT0M1 

FCTCM(2)=0.5*DTOM2 

CTCMB=DT/P6NAS 

HCT0MB=C.5*DT0MB 

CALL  STEP 

IF{  INDEX. GT.O)  GO  TO 


250 


FIRST  HALF  OF  TIMESTEP 


INDEX=1 
1  =  1 


\ 


100 


IF(LCUT(I )  .EQ.2)  GC  TO  215 

IF(LCUTd)  .NE.O)  GG  TO  220 

RXK(I )=RX( I  ) 

RYKU )=RY( I ) 

RZK(I)=RZ( I ) 

RX(I)=RX( I )+DTOQ*(HDTOMB*FX (I )+VX(I ) ) 

RY( I)=RY(I )+DTOD*(HDTCME*FY (I)+VY(I)) 

RZ(I)=RZ(I)+DTOD*(hDTCMB*FZ<  I  )  +  VZ(I) ) 

GC  TO  220 
215  RX(I)=RX( I)+DTOD*VX(I  ) 

RV(I)=RY(I )+DTOD*VY(I ) 

RZ(I )=RZ(I )+DTOD*VZ(I ) 
220  DC  240  1=2, LL 

IF(LCUT(I )  .EQ.2)  GO  TO  235 

IF(LCUTU)  .GT.O)  GC  TO  240 

J  =  INF(I  ) 

RXK(I)  =  RX(I  ) 

RYKU  )  =  PY(  I  ) 

RZK(I )=RZ( I  ) 

RX(I)=RX(I )+OTOD*lHDTOM(J)*FX( I  )+VX(I  )) 

PY(I)=RY(  I )+OTOQ*(HDTCM(J)*FY( I  )+VY( I ) ) 

RZ(I)=RZ( I )+QTOD*(HDTCM(J) *FZ( I  )+VZ( I  ) ) 

GC  TO  2^0 
225  PX(I)=RX( I)+DTOD*VX(I  ) 

PY(I  )=RY(I  )+DTOD*VY(I  ) 

PZ(I)=RZ(I)+DTOD*VZ(I ) 
24C  CCNTINUE 

GC  TO  2C0 
C 

C  SECCND    HALF    OF    TIM5STEP 

C 

250    CCNTINUE 

II\CEX  =  0 

ENAX=0.0 

TINE=TINE+DT 

NT=NT+1 

1  =  1 

IF(LCUT(I )  .NE.O)     GO    TO    260 

\,SS=VX(  I) 

VX( I)=VSS+HCTOMB*FX(I  ) 

RXI)=RXK(  I  )+(VX(I)+VSS)*HDTOD 

^SS=VY(  I) 

VY(I)=VSS+HDT0M3*FY(I  ) 

RY(I)=RVK(I)+(VY(  I)+VSS)*HDTOD 

\,SS=VZ(  I) 

VZ( I)=VSS  +  HCTOMB*FZ(I  ) 

RZ(I)=RZK(I  )+(VZ(  D+VSS  )*HDTOD 

FX<1)=0.0 

FY(  1)=0.0 

FZU)=0.0 

PKE(I)  =  VX(  I  )*VX( I  )+VY(I )*VY( I ) +VZ(  I ) *VZ (  I  ) 

EPAX=PKE(  I) 

K  =  I 

PKE(I)  =  PKE(  I)*HBM4S 
26C    CC    290     1=2, LL 

IF(LCUT(I )  .NE.O)     GC    TO    29  0 

J=IfF(I  ) 

VSS=VX( I) 

VX(I)=VSS  +  HCTOM(J  )*FX<  I  ) 

RX( I)=RXK(I  )  +  (VX(  I  )+VSS)*HDTGD 

VSS  =  VY(  I  ) 

VY(I )=VSS+hDTOM( J )*FY( I ) 

RY(I )=RYK(I  )  +  ( VY(  I  )  +  VSS)*HDTOD 

VSS  =  VZ(  I) 

VZ(I)=VSS  +  HDTOM(J  )*FZ( I) 

RZ(I)=RZK(I  )+(VZ(  I  )+VSS)*HQTOD 

FX(I  )  =  0.0 

FY(  I  )=0.0 

FZ(I  )=0.0 

FKE( I)  =  VX(  I  )*VX( I  )+VY(I  )*VY(I )+VZ(  I  J*VZ(  I) 

IF(PKE(  I  )  .LT.EMAX  )    GO    TO    280 

EMX  =  PKE(  I  ) 


10L 


280 


C 

c 

c 


28 
29 


29 

29 


3C: 
3C" 


31C 


32C 


35C 


36C 


37 
38 


4CC 


K  =  I 

IF  (  IMF(  I)  .ECU     GC    TO 

PKEU  )=FTMAS2*PKE  (I) 

GC    TO    290 

PKE( I  )=FTMAS1*PKE  (I) 

CCNTINUE 


282 


CALL 
CALL 
FINIS 
1F<  F  I 

IF(  IM 
ETEST 
GC  TO 
ETEST 
CCMI 
IF(ET 
NCT  =  N 
IF(NC 
NCT  =  0 
CC  30 
IF(LC 
P  F  E  ( I 
CCMI 
C/JLL 
CC  35 
IF(LC 
IF(RY 
IF(RX 
IF(RZ 
IFIRY 
LCUT( 
P>F(  I 
RVFd 
RZF(I 
PTE(  I 
GC  TO 
I  F  <  I . 
IF(P 
PTEU 
IF  (PT 
LCLT( 
RXFU 
RYF(I 
RZF(I 

ccmi 


TEST  FOR  PRINT  OR  TERMINATION 

GETIME(IET) 

SETIME 

H=  IET  *2. 65-05  +FINISH 

NISH.GT.STOPER)  GO  TO  373 

F(K).EC.l)  GO  TO  293 

=HTMAS2*EMAX 

294 
=HTMAS1*EMAX 
NUE 

EST. LE. TESTE  )  GO  TO  380 
CT+1 
T.NE.NCTT)  GC  TO  360 


5  I 
UT( 
)=C 
NUE 
ENE 
0  I 
UT  ( 
) 


=  1,LL 

I) .GT.O)  GO  TO  305 

.0 


(I 
(I  ) 
(I  ) 
(I  ) 
I)  = 
)  =  R 
)=R 
)  =  P 
)  =  P 
35 
EQ. 
Pt( 
)  =  P 
E(I 
I)  = 
)=R 
)  =  R 
)=R 
NUE 


RGY 

=  1» 

I)  . 

.LT 

.LT 

.LT 

.LT 

1 

X(I 

Y(  I 

Z(  I 

KE( 

0 

1) 

I)  . 

KE< 

)  .L 

2 

X(I 

Yd 

Z(I 


GC  TO  350 
AND.VY(I) .LT.O.O) 
OR.RXd  )  .GT.RXBM) 
CR.RZ(I).GT.RZBM) 
GO  TO  350 


LL 

GT.O) 

.RYBL. 

.RXBL. 

.RZBL. 

.RYBM) 

) 

) 
) 

I)+PPE(I) 

GO  TO  310 

GT.O.OOOL)  GO  TO  350 

i)+ppe  (i) 

T.O.O)  GO  TO 


GO  TO  32C 
GO  TO  31C 
GO  TO  31C 


350 


TEST  FOR  CUTFUT  OF  CURRENT  STATUS 


CT=FDTI/SQRT(EMAX ) 
IF(NT.GE.NTT)  GO  TO  380 

GO  TO  NEXT  TIMSSTEP 

IFiNS.NE.NT  )  GO  TC  103 

GC  TO  400 

ISFUT=1 

GC  TO  4C0 

NSS  =  1 

IF(NT.GE.NTT)  ZSTCP=1.0 

CALCULATE  POTENTIAL  ENERGY  AND  PRINT  STATUS 
OF  THE  CURRENT  TIMESTEP 

CCNTINUE 

TFCT=0.0 

PFCS=O.C 

FNEG=0.0 

CC  450  1  =  1, LL 

IF(ISEQ(I  )  .EQ.l)  ISEQ(I)=2 


102 


IF(LCUTd)  .NE.O)  GC  TO  450 

PFEd)  =  C.O 

PTE(I)=0.0 
45C  CCNTINUE 

CALL  ENERGY 
595  PTE(1)=PKE(1)+PPE(1) 

TFKE=PKE( 1 ) 
6CC  CC  625  1=2, LL 

TFKE=TFKE+PKE( I ) 

IF(PPE(  I)  .GT.O.O)  GO  TO  615 

FNEG=FNEG+PPE( I) 

GC  TC  620 
615  FFCS=PPCS+PPE( I ) 
62C  PTE(  I)=PK5(  D+PPE  (I) 
625  CCNTINUE 

TE=TPKS+TFOT+BIND 

IFdNSS.NE.O)  .ANC.  ( ZSTOP . EQ .0 . 0  )  )  GO  TG  95J 

CALL  KPFINT 

WRITE  (6,9681) 
730  CC  750  1=1, LL 

IF  (PTE(  I)  .LT.0.1  )  GO  TC  750 

IFdSEQ(i).EQ.O)  ISEQ(I)=1 

CXd)=RX(  D-RXI  (I  ) 

CV(I  )=RY(  D-RYI  (  I  ) 

CZ(I)=RZ(  D-RZI  (I  ) 

WRITE (6,  96 82)  I,DX(I),DY(I),DZ(I),VX(I),VY(  I),VZ(  I), 
XPKE(I),PPE(I),PTE(I),RX(I),RY(I),PZ(I),LCUT(I),ISEC(I) 
75C  CCNTINUE 

WFITE(6,96  9C)  TPKE , TP CT ,TE , PPOS  ,  PNEG 

IF(  ISHUT.NE.O)  GO  TO  972 

IF(ZSTOP.NE.  0.0)  GO  TC  973 

WRITE(6,963C)  NPAGE 

NFAGE=NFAGE+1 

NS=NS+NC 
C 

C  GO  TO  NEXT  TIMESTEP 

C 

GC  TO  ICO 
95C  CCNTINUE 
C 

C  TERMINATE  REPLICA  ANC  TFEN 

C  PRINT  FINAL  STATUS  OF  CURRENT  PASS 

C 

J  =  l 

K=C 

IK  =  0 

CC  960  1=1, LL 

IF  (LCUT(  I)  .LT.2)  GO  TC  960 

IF(VYd  KGT.VTEST  )  GO  TC  96C 

K=K  +  1 

I5MJM(J)*I 

LNLMd  )=LNUMd  )+l 

J=J  +  1 
96C  CCNTINUE 

K  K  =  K 

NSFUTdPA3S)=KK 

CALL    KPPINT 

IF(KK.EC.O)    GO    TO    972 

WRITE(6  ,9813)    E V ,  ( I  SNUM  ( I ) , I  =  1 , KK  ) 

WFITE(6,960C) 

WFITE<6,9695) 

DO  963  11=1, KK 

I  =  ISNUM(II  ) 

W P  ITE( 6,980 C) IPNUM,I ,BUL,TAR, I  OR  ,SUR(1)  , VX (  I )  ,  VY (  I  ) , VZ 
X(I),PKE(I),PPE(I),PTE(I),EVR,COX,CCZ 

968  CCNTINUE 

CC  969  1=2, LL 

I  F(LCUT(I ) .EQ.2   .AND.   VY( I ) • LT.VTS )  WR ITE (6 ,9300 ) 
X   I  PNUM  ,1  ,BUL,  TAR  ,  IOR,SUR(  1)  ,VX(  I  )  ,V  Yd  )  ,VZd  )  ,PK  E(I  )  , 
X   FPE(  I)  ,PTE( I)  ,EVR 

969  CCNTINUE 


97C 


972 


972 


974 
975 


99C 


10QC 


IOC 


DETERMINE  SPUTTERED  ATOM  ANGLE 

F(KK.EQ.O)  GO  TO  972 

RITE  (6,9600) 

FI7E<6,9805) 

C  970  1  =  1, KK 

=  ISMIM(I) 

£C  =  A7AMSGR7(VX(L)**2  +  VZ(L)**2)/(-VY(L) 

NGLE=PAC*RADDEG 

=IFIX((ANGLE/RINCR)+1.0) 

F(IMF(  L)  .EC.L)  IEOXK  J)  =  IBCX1(  J)  +1 

F(IMF(L)  .EC. 2)  IB0X2( J)=I3GX2(J)  +1 

RITE  (6,9815)  L, ANGLE 


)  ) 


CMINUE 
FUISHLT 


EQ.l) 


RITE(6,963C)  NPAGE 

PAGE=NFAGE+1 

C  TC  990 

F(  ISHUT.EQ.O)  GO  TO 

FITE(6,9965)  IPNUM 

FITE(6,960G) 

RI7E(6  ,9960) 

SFLT( IFASS  )=0 

F4SS=I FASS-1 

5hCTS=NSH07S-l 

C  70  975 

F(ZS7CF.EC.O)  GO  70 

RI7E(6,9965)  IPNUM 

FI7E(6  ,96JC) 

RI7E(6,9955) 

SFL7(IPASS  )=0 

FASS=IFASS-1 

<FC7S=NSHQ7S-1 

C    7C    975 

RI7E(6,982C)    EV 

PI7E(6,963C)    NPAGE 

FAGE=NPAGE+1 


CR.        (ZS70P.EQ.1.0) )     GO    7C    972 


973 


974 


DEFINE    NEW    REPLICA    OR    7ERMNA7E 


IPASS=IPASS+1 

IF(ISHUT.NE.O) 

E2=BZ+SPDZ 

IF  (EZ.LE.BZNAX) 

EX=EX+SPDX 

IF(BX.G7.BXMAX) 

IF(  I0R.EQ.1C0) 

IF(  IGR.EC.lll) 

EZ=EZMIN 

IF(  I0R.EQ.1C0) 

IF  (IOR.EQ.lll) 

GC    70    4  0 


GO    70    1000 

GO    70    40 

GO    70    1000 
BZMIN=BZMIN+SPDZ 
3ZMIN=EZMIN+SPDZ/2 

BZMAX=BZMAX-SPDZ 
BZMAX=EZMAX-SPDZ/2 


NRAN,EVR 
( ( IPNUMS( 


I)  ,NSPU7(  I)  )  ,1  =  1  ,NRAM 


PRIN7  SUMMARY  OF  ALL  SHC7S 

CCN7INUE 

CALL  KPFIN7 

WRI7E  (6,9600) 

NPAIS  =  IPASS-1 

WPI7E(6,986C) 

WP  I7E(6,987C) 

KSLM=0 

K2SUM=Q 

DC  1005  I=1,NRAN 

KSLN=KSLM+NSPU7(I ) 

K2SLM=K2SUN+NSPU7(I )**2 

CCN7INUE 

FRAN=NRAN 

FRANR=1.0/FRAN 

AVE=FL0A7(KSUM)*FRANR 

VAR=FLCA7(K2SUM)*FRANR-AVE**2 

S70=VAR*FRANR 


134 


ETHR,VTEST,KSUM 
AVE  ,VAR,NRAN,STD 


(  I,  LNUMd  )  ,1=1. LL) 


10C7 
1008 
IOCS 


1011 


STC=SQRT(STC) 
WRITE  (6,9600) 
WFITE(6,9875) 
WPITE(6,985C) 
WFITE  (6,96C0) 
WFITE  (6,9865) 
WFITE(6,987C) 
WRITE  (6,9600) 
WRITE  (6,9920) 
WRITE  (6,9935) 
IF(  IOR.NE.100) 
WRITE(6 ,995C)ISEEDA 
IF(IOR.NE.llO)  GO  TC 
WRITE(6,995G)  ISEECB 
IF( IOR.NE.lll)  GO  TO 
WFITE(6,995C)  ISEECC 
CONTINUE 

WFITE(6,963C)  NPAGE 
NPAGE=NPAGE+1 
CALL  KPNT 


FOLD  FINAL  SUMMARY 

IF  (IOR.EQ.100)  GO  TO  1011 

IF  (I0R.EQ.110)  GO  TO  1012 

IF  (IOR.EQ.lll)  GO  TO  1013 
C^LL  KFCLDA 
GC  TO  1C14 


IREAD 

(  IFFM(J) , IFM(J) , J  =  1,IREA0  ) 

GO  TO  1007 


1008 


1009 


1012  C*LL  KFCLOB 


1012 
1014 


GC  TO  1014 
C4LL  KFCLDC 
WFITE  (6,9600) 
WFITE  (6,9840) 
WFITE  (6,9845) 
WPITE(6  ,9870) 
WRITE  (6,96C0) 


( I, LNUM(I) ,1=1, LL) 


DETERMINE  RING  SPUTTERING  PROBABILITY 


101 


CALL  KPPCB 
WFITE(6 ,9890) 
WFITE(6,98SC) 
CC  1020 
J=  1  +  1 
K  =  I+2 

WFITE(6,98  85)I,DNN(I) ,KATOM ( I ) , KSPUT( I ) , PROE  ( I ) , J , 
DNN( J) ,KATON( J)  ,KSPUT(J)  ,PROB( J)  ,  K , DNN ( K  )  ,KATOM ( K  )  , 


I=1,KRING,3 


X 

X   KSPUT(K)  ,FROB (K) 

1020  ccraiNUE 

WRITE  (6,96C0) 

WRITE  (6,9630)  NPAGE 

NFAGE=NPAGE+1 

C/i  LL  KPNT 

kFITE  (6,96C0) 


DRAW  ISOTOPE  HISTOGRAMS 

1025  CC  1030  I=1,IRANGE 

FFECKI  )=IBCX1(  I  ) 

FFEG2(I  )  =  IBCX2( I  ) 
103C    CCNTINUE 

PANGE(1  )=RMIDPT 

CC  1040  I=2,IRANGE 

J  =  1-1 

RANGE(I)=    RANGE(J  )  +  RINCR 
104C    CCNTINUE 

IN=-IRANGE 

CALL  HI  STO  (  IN,  FRE  CI,  RANGE,  TITLED 


WFITE 
WRITE 
WFITE 


(6,96C0) 
(6,9895) 
(6,9900) 


RINCR 


IF  ( IRANGE.EQ.9)   WRITE  (6,9925) 


105 


IF  (IRANGE.EQ.18)  WR  I  TE ( 6 , 9920 ) 

WRITE  (6,9910)  < I BCX1 ( J ) , J= 1  ,  IR ANGE ) 

WRITE    <6,96C0) 

WFITE    (6,96C0) 

C4LL    HIST0(  IN,FREC2,R£NGE,TITLE2) 

WRITE  (6,9600) 

WRITE  (6,9895) 

WRITE  (6,9900)  RINCR 

IF  (  IRANGE.EQ.9)   WRITE  (6,9925) 

IF  (  IRANGE.EQ.18)  WRITE (6 , 9920 ) 

WRITE  (6,9910)  ( I 8GX2 ( J ) , J= 1 ,  IRANGE ) 

WRITE  (6,96C0) 

WRITE  (6,9630)  NPAGE 


NFAGE=NPAGE+1 
IF(ISHUT.EG.O) 
9999  STCP 
EfsC 


GO  TO  5 


106 


SLEROUTINE    SETUP 


C 

c 
c 
c 
c 
c 


THIS  S/R  FINDS  THE  STARTING  POINT  FCR  EACH 
IMPACT  POINT  IN  THE  TARGET  AREA 


1GQ 


CCM 
CCN 
X  FN 
CCM 
CCM 
COM 
1  =  1 
Pl  = 


oc 


3 
97 


CC 


cc 


F( 

IX 

iy 

12 
XM 
YM 
ZM 
XNA 
YNA 
ZNA 
CC 
IF< 
IF( 
IF( 
IF( 
IF( 
IF( 
IF( 
CFX 
CRY 
CR2 
OIS 
IF( 
GC 
CCN 
GC 
CCN 
R2S 
R2  = 
RIP 
X  ( 
ALF 
ALF 
R2  = 
RX( 
RY( 
RZ( 
FET 
PX( 
RY  ( 
R2( 
FCR 

RET 
ENC 


MCN/COM1 
MCN/CCM2 
,PFIV,TP 
MCN/C0M6 
yCN/CCM7 
MCN/COMl 

Rl-FAC 
Rl  .LT.O. 
=R1*CCXI 
=R1*CCYI 
=R1*CCZI 


N 
N 
N 
X 
X 
X 
195 


R1X- 
P1Y- 
R12- 
R1X  + 
P1Y  + 
P1Z  + 
J  =  2, 


LE 
GT 
LE 

GE 

LE 


LCUT  (J) 

RY(J).GE 

RY(  J  ) 

RX(J  J 

RX(  J) 

RZ(J) 

RZ(J  ) 

=RX t J)-P 

=RY(J  )-R 

=R2( J)-P 

T=DRX*#2 

CIST.GE. 

TO  2C0 

TINUE 

TO  1GO 

TINUE 

Q=(RXI(  J 

SCRT  (R2S 

2=( R1X-R 

R12-R2S ) 

A=RlF2/( 

A2=ALFA* 

R2*ALFA+ 

1)=R3*CC 

1)=R3*CC 

1)=R3*CC 

LRN 

1)=RXS 

1)=RYS 

1)=R2S 

MAT(41h 

ITE  (  6, 

LRN 


/RX(60C)  ,RY(600)  ,RZ(60C) ,LCUT  (  60 C  )  , LL , LD 
/RCE,RCE2,RQEM,AC,PAC,FPTC,PTC,PFPTC  ,FPTC, 
CT 

/CQXItCOYI,CCZI,RXS,RYS,RZSfFAC 
/R1,LSS,SPX,SPZ,C0X,CCY,CCZtYLAX  (23)tTLL 
1/RXI (60C) tRYI (600) fRZI  (600) 


C)  GO  TO  300 

+RXS 

+RYS 

+RZS 

ROE 

ROE 

ROE 

ROE 

ROE 

ROE 

LL 

GT.O)  GC  TO  195 

.YMAX)  GC  TO  195 

.YMINJ  GO  TO  195 

.XMAX)  GO  TO  195 

.XMIN)  GC  TO  195 

.ZMAX)  GO  TO  195 

.ZMIN)  GO  TO  195 

IX 

1Y 

1Z 

+DRY**2+DRZ**2 

R0E2)  GO  TO  195 


)-RXS )**2+(RYI ( J)-RYS)**2+(RZI( J  )-RZS) **2 

Q) 

XS)*(RXI(J)-RXS)+(R1Y-RYS)*(RYI ( J)-RYS  )  + 

*(RZI ( JJ-RZS) 

R1*R2) 

ALFA 

SQRT(R2SQ*ALFA2-R2SQ+RCE2) 

XI+RXS 

YI+RYS 

ZI+RZS 


THIS  DIC  NCT  LOCATE  A  PROPER  START  POINT) 
9700) 


SLEROUTINE  VvARM 


THIS  IS  A  MGDEL  3  WARMER 

EACH  ATOM  IS  DISPLACED  AN  AMOUNT  CETERMINEC 
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C  BY     ITS    THERMAL    AMPLITUDE 

C 


c 


CCNMGN/COM1/RX(600),RY(60J)  ,RZ ( 60  0 ) , LCUT (60 C  )  ,LL,LD 
CCPMCN/COM13/  RAN (1800) , I  RAN, TA SN , TASF , L L2  ,  LL3 
CCNM0N/C0M14/  AXT (6G0 ) , AYT ( 600 )  ,AZT(600) 
AXT(1)=0 
AYT(1)=C 
A2T(1)=C 

CALL  NCFMAL(IFAN,RAN,LL3) 
CC  1100  1=2, LL 
£XT(I)=RAN(  I) 
AVT(I)=PAN(  I  +  LL) 
AZT(I )  =  RAN(  I+LL2) 

FX(I)=RX(I )+TASN*AXT( I )+TASF*AXT<  1-1) 
RV(I)-RY(I  )+TASN*AYT(I)+T4SF*AYT(I-l) 
R2(I)=RZ(I )+TASN*AZT( I )+TASF*AZT(  1-1) 
1100  CCNTINUE 
RETURN 
END 


SieRCUTINE  KF10 

c 

c 

c 

LATTICE 

c 

LATTICE 

c 

THE  Z  D 

c 

c 

GENERATCR  FCR  A  FCC  (100)  ORIENTATION 

IS  DEVELQPEC  IN  THE  X  QIREC7ICN,  THEN 

IRECTION,  THEN  FINALLY  THE  Y  DIRECTICN 

CCVMCN/CCM1/RX(60C) ,RY(600)  ,RZ ( 60 C ) , LCUT ( 6  0 C  )  , LL , LD 
CCyMCN/CCM2/RCE,RCE2,RQEM,AC,PAC,PPTC,PTC,PFPTC,FPTC, 

X  FN,PFIV,TPCT 
CCMfCN/C0M4/IX,IY,IZ,  IXP, I YP , SCX  ,  SCY , SC Z  ,  ID EEP , 

X  F  P  X  R  P  2 

CCNMCN/CCM7/R1,LSS,SPX,SPZ,CCX,CCY,CCZ,YLAX(20) ,  ill 

CCNMCN/C0M11/RXI (6JO),RYI(600),RZI(60  0) 

CCMMCN/C0M22/  SPDX,SPDZ,03XNAX,CeZMAX,SSCZ 

CCMMCN/CCM2  5/IMF(600) , ISOGEN 

CCMMCN/C0M3C/  ISEECA,  ISEEDB,  ISEECC 

C1NENSICN  R(600) 
9001  FCFMAT(«  I  S EED= ' ,  1 10/ /  ) 

9690  FCFMAT(25XT 'SUMMARY  OF  ATOMS    FCC  (100)  ORIENTATION' 
X   /,35X,'MASS  1  IS  THE  HEAVIER/LEAST  ABUNDANT  ISOTOPE' 
X   //,3('    ATOM   MASS        POSITION  '),/) 

LC  =  C 

N  =  600 

CC  3  1=2, N 

INF(I)=2 
2  CCNTINUE 

ISEED=ISEEDA 

WRITE  (6,90C1)  I5EED 

WRITE  (6,9690) 

IF  (ISOGEN. EQ.O)  GO  TO  4 
C 

C  GENERATE  RANDOM  MASSES 

C 

CALL  GGU3( ISEED,N,R) 

CC  4  1=2, N 

IF(R(I)  .LE.C.309)  IMF (  I  )  =  1 

4  CCNTINUE 

CC  5  1  =  3, IY 

5  YLAX(I)=0.0 
C 

C  DETERMINE  ORIENTATION  PARAMETERS 

C 

SCX=1.0 

SCY=1.J 

SCZ=1.0 

SFCX  =  0.  1 
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c 
c 
c 


20 


57 

56 
5<5 


6C 


ac 


1C 

n 

12 


SFCZ  =  0.  1 

CEXNAX=0.501 

CEZMAX=1-001 

YL£X(1)=-0.12 

YLAX(2)=-0.03 


M  =  2 
J7  =  0 
Y=-SCY 
CC  60  J 
Y=V+SCY 
K7  =  C 
Z=-SCZ 
CC  59  K 
Z=Z+SCZ 
IT  =  0 
X=-SCX 
CC  58  I 

x=>+scx 
iit=it+ 

IFUTT- 
RXI(M)= 
RYI (M)= 
RZI(M)= 
V=N  +  1 
CCNTINL 

n  =it+ 

CCNTINL 

KT=KT+1 

CCNTINU 

IF(JT.E 

Jl  =  JT 

IF(JT.5 

CCNTINU 

LL=M-1 

IF(LD.E 

IF{CCY. 

C-C  TC  ( 

R1=RCE/ 

GC  TO  1 

ISFX=RB 

ISFZ=RE 

N=(IX/2 

IF(  IX-( 

N=N+ISP 

LCLT(N) 

Pl=2.0 

LCLT(2) 

GC  TC  1 

V  =  2 

RXI(M)= 

RVI(M)= 

RZI (M)= 

P1=(R0E 

CCNTINU 

RETURN 

Rl=1.5 

RETURN 

ENC 


GENERATE  THE  LATTICE 


=  1,  IY 


=  1,IZ 


=  1,  IX 

JT  +  KT 

(ITT/2J*2)  57,30,57 

X 

Y  +  YLAXU) 

Z 


CO)  ILL=M-1 
+  l 

C.IDEEP)  LC=M-l 
E 


)  LD=LL-1 
0)  GO  TC  120 
,80,90)  ,LSS 
(COY) 


SFZ+ISPX/2+3 
2)*2)  83,85,83 


CO 

EQ. 
70 
A8S 
CO 
X 

z 
)*I 

IX/ 
X/2 
=  1 

=  -1 
00 


RBX*SCX 
-SCY+YLAX( 1) 
IPBZ+1.0)*SCZ 
+SCY)/A3S(C0Y) 


C 

C 

c 

C 

c 

C 


SLEPOUTINE  KF110 


LATTICE  GENERATOR  FOR  A  FCC  (110)  ORIENTATICN 
LATTICE  IS  DEVELOPED  IN  THE  X  DIRECTICN,  THEN 
THE  Z  DIRECTION,  THEN  FINALLY  THE  Y  DIRECTICN 
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90 
96 


01 


CC 
CC 

X  F 
CC 

X  F 
CC 
CC 
CC 
CC 
CC 
CI 
FC 
FC 


11 
12 
21 
22 


LC 
M 

CC 
IN 
CC 
IS 

fcP 
*F 
IF 


MMCN/ 
MMCN/ 
N,PFI 
NMCN/ 

ex, Re 

NMCN/ 

NMCN/ 

MMCN/ 

NMON/ 

MMCN/ 

MENSI 

FMAT( 

PMAT( 

/,35X 

//,3( 

=  C 

6C0 

3  1  = 
F(I)  = 
MINU 
EED=I 
ITE  ( 
ITE     ( 

(ISO 


C0M1/RX(60C),RY(600),RZ(6CC),LCUT(6QC),LL,LD 

CCN2/RCE,RCE2,RaEM,AC,PAC,PPTC,PTC,PFPTC,FPTC, 

V,TPCT 

C0M4/IX,IY,IZ,IXP,IYP,SCX,SCY,SCZ,IDEEP, 

Z 

C0M7/RlfLSS,SPX,SPZ,C0X,CCY,C0ZfYLAX(20) ,ILL 

C0M11/RXI (60  J) ,RYI(600) ,RZI (633) 

CQM23/  SPDX,SPDZ,DBXNAX,DBZNAX,SSCZ 

C0M25/IMF(6CC) , ISOGEN 


COM30/  ISEECA, 
CN  R(600) 
»  ISEED=» ,  110//) 
25X,  'SUMMARY  OF  ATOMS 
IS  THE 
MASS 


ISEEOBt  ISEEDC 


,«MASS  1 
'    ATOM 


FCC 

HEAVIER/LEGIST 
POSITIGN 


(1L0)  CRIENTATICN' 
ABUNDANT  ISOTCPE' 
'),/) 


2,N 

2 

E 

SEEDS 

6,9001)  ISEED 

6,9690) 

GEN.EQ.O)  GC  TG  4 

GENERATE  RANDOM  MASSES 


CALL  GGLB(  ISEED, N,P) 

CC  4  1=2, N 

IF(R(I)  .LE.C.309)  IMF(  I  )=1 

CCNTINUE 

CC  5  1  =  3, IY 

YLAXd  )=0.C 

CETERMINE  ORIENTATION  PARAMETERS 

PC=1.0/SQRT(2.0) 

SCX=RO 

SC>=RO 

SCZ=1.0 

SPCX=0.2*SCX 

SFCZ=0.2*SCZ 

CEXNAX=SCX+3.01 

CEZMAX=SCZ  +  C01 

YLAX(1)=-0.15 

YL/5X(2)=-3.Q3 

GENERATE  THE  LATTICE 


M=2 

JT=0 

Y=-SCY 

CC  60  J  =  l,  IY 

Y=Y+SCY 

KT  =  C 

2=-SCZ 

CC  59  K-l,  IZ 

Z=Z+SCZ 

IT  =  0 

X=-SCX 

CC  58  1  =  1, IX 

X=X+SCX 

IF(IT-(IT/2)*2)  21,11 

IF(JT-{ jT/2)*2)  57,12 

IF (KT-(KT/2)*2)  57, 3C 

IF<JT-< JT/2)*2)  22,57 

IF(KT-(KT/2)*2)  30,57 

PXI (M)=X 

RYI (M)=Y+YLAX( J) 

RZI(M)=Z 

N=N  +  1 


21 
57 
57 
22 

3C 


110 


57 


60 


7C 

80 


90 


ICC 
11C 
12C 


CC 
IT 
CC 
KT 
CC 
JT 
IF 
CC 
LI 
IF 
IF 
GC 
Rl 
GC 
15 
IS 
N  = 
LC 
Rl 
LC 
GC 
N  = 
RX 
RY 
RZ 
Rl 
CC 
RE 
PI 

EN 


NTINUE 

=  IT+1 
NTINUE 
=  KT  +  1 
MINUE 

=  JT  +  1 
(JT.EC.IDEEP)  LC=M-1 
NTINUE 
=  N-1 

(LD.EC.O)  LD=LL-1 
(COY.EQ.O)  GO  TO  120 

TO  (  70,80,93) ,LSS 
=RCE/ABS(COY) 

TO  100 
FX=SPX 
FZ=SPZ 

(IX  +  U/2*  ISPZ/2+ISPX+2 
LT(N)=1 
=RCE/A8S  (COY) 
LT(2)=-1 

TO  100 


I  (M) 
KM) 

I  (M) 
=  (PG 
NTIN 
TURN 
=  1.5 
RET 
C 


=(SPX+1.0)*SCX 
=-SCY+YLAX(  1) 
=(SPZ+1.0)*SCZ 
E+SCY)/ASS(COY) 
UE 


URN 


SLERCUTINE  KF111 


LATTICE  GENERATCR  FCR  A  FCC  (111)  ORIENTATION 
LATTICE  IS  DEVELOPED  IN  THE  X  DIRECTION,  THEN 
THE  Z  DIRECTION,  THEN  FINALLY  THE  Y  DIRECTICN 


900 
969 


CC 
CC 

X  F 
CC 

X  R 
CC 
CC 
CC 
CC 
CC 
DI 
FC 
FC 


LC 

N  = 
DC 
IV 
CC 
IS 

WP 

WR 
IF 


MM0N/C0M1 

NMCN/CCM2 

t  ,PFIV,TP 

NMCN/C0M4 

EX,RB2 

NMCN/CCM7 

MMCN/C0M1 

fMCN/C0M2 

PMCN/CGM2 

NMCN/CCM3 

MENSICN  R 

RMATC  IS 

FMAT(25X, 

/,35X,'MA 

//,3(  • 

=  C 

600 

3  1=2 
F(I)=2 
NTINUE 
EEC=ISEEC 
ITE  (6,90 
ITE  (6,96 

(ISOGEN. 


N 


/RX(60C)  ,RY(600)  ,RZ(6CC) ,LCUT ( 60 C ) , LL , LD 

/PCE,RCE2,R0EM,AC,PAC,PPTC,PTC,PFPTC,FPTC, 

CT 

/IX,IY,IZ,IXF,IYP,SCX,SCY,SCZ,IDEEP, 

/R1,LSS,SPX,SPZ,CGX,CCY,CGZ,YLAX(2C),ILL 

1/RXI(600),RYI(600),RZI(60J) 

3/  SPDX,SPDZ,DBXMAX,Ce2yAX,SSCZ 

5/IMF(6CC) , ISOGEN 

0/  ISEEDA,  ISEEDB,  ISEECC 

(600) 

EED=« , 110//) 

•SUMMARY  OF  ATOMS    FCC  (111)  CRIENTATICN' 

SS  1  IS  THE  HEAVIER/LEAST  ABUNDANT  ISOTCPE' 

ATOM   MASS        POSITION  •),/) 


C 

01)  ISEED 

90) 

EQ.O)  GG  TO  4 


GENERATE  RANDOM  MASSES 


C£LL  GGLB( ISEED, N , R ) 
CC  4  1=2, N 


i  l  i 


30 


43 


58 

EC 

60 

7C 
8C 

90 


I  F  CR ( I ) .LE.  0.309)  IMF(I)=1 
CONTINUE 
DC  5  1  =  3,  IV 
YIAX( I)=0.0 

DETERMINE  ORIENTATION  PARAMETERS 

SCX=1.0/SQRT(2.0) 

SCY=2.0/SQRT(3.0) 

SC2=SQRT(1.5) 

SSCZ=2.0*SCZ/3.0 

SPCX=0.2*SCX 

SFCZ=0.2*SSCZ 

CE>MAX=SCX  +  C01 

C£ZMAX=SSCZ40.01 

YLAX( 1) =-0.C5 

YL*X(2)=-0.C1 


N=2 
JT  =  C 
Y=-SCY 
CC  6J  J 
Y=Y+SCY 
J1S=JT+ 
Z=-SCZ 
KT  =  C 
CC  59  K 

z=z+scz 

IT  =  C 

x=-scx 

DC  58  I 

X=X+SCX 

IN=IT+J 

IF(IN-( 

F*I(M)= 

RYI  (M)  = 

IFUT-3 

JTT=JT 

J1T=JTT 

IF(JTT) 

J7T=JTT 

ZF=JTT 

RZI(M)= 

GC  TO  5 

RZI(M)= 

CCNTINU 

IT=IT+L 

CCNTINL 

KT=KT+1 

CCNTINL 

JT=JT+1 

IF(JT.E 

CCNTINU 

Ll=P-l 

IF(LD.E 

IFtCOY. 

GC  TC  ( 

R1=RCE/ 

LCtT(2) 

GC  TO  1 

I<FX=SP 

ISFZ=SP 

N=IX*(I 

LCLT(N) 

R1=RCE/ 

GC  TC  1 

l*  =  2 

PXI(y)= 

RYI  (M)  = 


GENERATE  THE  LATTICE 

=  1,IY 
JT/3 

=  1,IZ 


=  1,  IX 

TS  +  KT 

IN/2)*2)  57T3C,57 

X 

Y+YLAX( J ) 

*(JT/3)1  41t45,41 

-3 

43,45,42 
+3 

Z+ZP*SSCZ 

C 

Z 


C.IDEEP)  LD=M-1 
E 

CO)  LD  =  LL-L 
EQ.O)  GO  TO  120 

70,80,9:)) ,LSS 
ABS(COY) 
=  -L 
00 
X 

z 

SPZ/2)+ISPX+2 

=  1 

ABS(COY) 

00 

SPX*SCX 
-SCY+YLAX( 1) 
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RZI (M)=(SPZ+1.0)*SCZ 

R1=(R0E  +  SCY  )/ABS(COY) 
ICC    CCMINUE 
11C    RETURN 
12C    Rl«l. 5 

RETURN 

end 


c 
c 
c 
c 
c 
c 
c 


c 
c 
c 
c 


SIERCUTINE  ENERGY 


THIS  S/R  CALCULATES  THE  MUTUAL  POTENTIAL  ENERGIES 
THIS  S/R  ALLOWS  ATTRACTIVE  POTENTIALS 
THIS  IS  THE  MUNDKUR  HIGH  SPEED  VERSION 

CCMfCN/CGMl/RX(60C) ,RY(600)  ,RZ ( 600  )  , LCUT (60C),LL,LC 
CCMMCN/CCM2/RCE,RCE2,R0EM,ACtPACtPPTC,PTCtPFPTC,FPTC, 

X   FN,PFIV,TFOT 
CCMMCN/CaMS/EXATEXBtFXAtPEXAtPEXEtPFXATThERN 
CCMMCN/COM5/FX(60  0)  , FY (600)  ,FZ(6QC),PPE(600) 
CCNMCN/C0M22/R0EA,RGEBf R0EC,R0EC2fCP3,CPl»CF2,CP3,CF0, 

X   CFlfCF2fCGDl,CGD2tCGBl,CG62fCGF,  I  PROG 
1  =  1 

IFUCUT(I)  .NE.O)  GG  TO  200 
IF=I+1 

INITIAL  POTENTIAL  ENERGY  INTERACTIGN 
MAXIMUM  DISTANCE  BOUNDARIES 

XNAX  =  RX(I)+ROE 
YNAX  =  RY(I)+ROE 
ZNAX  =  RZ(I)+ROE 
XfIN  =  RX(I  )-ROc 
YMN  =  RY  (I  J-ROE 
ZMN  =  RZ(I  J-RQE 

TEST  TO  DETERMINE  IF  BULLET  AND  ATONS 
ARE  CLOSE  ENOLGH  TO  INTERACT 


DC  195   J=IP,LL 
IF(LCUT(J) .GT.O)  GC  TO 
IF(RY(J  ).GE.YMAX)  GO  T0 
IF(RY(J  ).LE.YMIN) 


IF(RX(J) 
IF(RX(J  ) 
IF(RZ(J) 
IF(RZ(J) 


GT.XMAX) 
LE.XMIN) 
GE.ZMAX) 
LE.ZMIN) 


GO 
GC 
GO 
GC 
GO 


TO 
TO 
TO 
TO 
TO 


195 
195 
195 
195 
195 
195 
195 


195 

2C0 


TEST  TO  DETERMINE  IF  THE  ATOMS  ARE  INSIDE  A 
CIRCLE  OF  RADIUS  ROE  FRCM  THE  BULLET 

CPX=RX( J)-RX(  I) 
CRV=RY( J)-RY( I) 
CRZ  =  RZ( J)-RZ(  I) 

C  IST  =  CRX**2+0RY**2  +  DRZ**2 
IF(DIST.GE.R0E2)  GC  TO  195 
CIST=SQPT(CIST) 

CALCULATION  OF  THE  POTENTIAL  ENERGY 

PCT  =  EXP (PEXA+PEXB*OIST  )-PPTC 

IF(PCT.LE.C  .00001  )  POT  =  0.0 

TFCT=TPCT+PCT 

SAVE=PFIV*PCT 

PPE(  I)  =  PPE(  D+SAVE 

PPE(J)=PFE( Ji+SAVE 

CONTINUE 

CCNTINUE 


113 


28C 


295 

3CC 


310 


CC  300  1  =  2, LD 
IF(LCUTd)  .GT.O) 
IP=I+L 

XNAX  =  RX(  I  )+ROE 
XNIN  =  RX(I  )-ROE 
YNAX  =  RY(I)+ROE 
YMN  =  RY  (I  )-ROE 
ZNAX  =  RZ(I)+ROE 
ZMN  =  RZ(I)-R0E 


GC  TO  300 


TEST  TO  DETERMINE  IF  THE  ATOMS 
ARE  CLOSE  ENOUGH  TO  INTERACT 


CC  295 

IF(LCUT 
IF(RY(J 
IF(RY(J 
IFCRXU 
IF(RX(J 
IF(RZ(J 
IF(RZ(J 
CRX=RX< 
DRY=RY( 
CRZ=RZ( 
CIST=DR 
IF(CIST 
CIST=SQ 


J=IPtLL 

(J) .GT.O) 

).GE.YMAX) 

).LE.YMIN) 

).GE  .XMAX) 

).LE.XMIN) 

) .GE.ZMAX) 

J.LE.ZMIN) 

J)-RX<  I) 

J)-RY( I) 

J)-RZ(I) 

X**2+DRY**2+CRZ**2 

.GE.R0E2)  GC  TO  295 

RT(DIST) 

CALCULATICN  CF  THE  POTENTIAL  ENERGY 


GC  TO  295 
GO  TO  295 
TO 
TO 
TO 
TO 
TO 


GG 
GC 

GC 
GC 
GO 


295 
295 
295 
295 
295 


PCT=EXP(EXA+EXB*DIST)-PTC 

IF(F0T.LE.Q.00C01 )    POT=0.0 

TFCT=TPCT+PCT 

SAVE=PFIV*PCT 

PPE(  I)  =  PPE(  D+SAVE 

PPE(J)=PPE( J)+SAVE 

CCNTINUE 

CCNTINUE 

DC  310  1=1, LL 

IF(LCUT( I) .EQ.O)  GC  TO  310 

TFCT  =  TPCT  +  PFE(I  ) 

CCNTINUE 

RETLRN 

ENC 


C 
C 
C 
C 
C 

c 
c 


c 

c 

C 


SLERCUTINE  STEP 

THIS  S/R  CALCULATES  THE  FCRCES  BETWEEN  ATOMS 
ATTRACTIVE  FORCES  ARE  ALLCWED  BETWEEN  ATOMS 
THIS  IS  THE  MLNOKUR  HIGH  SPEED  VERSION 

CCMMGN/C0M1/RX(60C) ,RY(600)  ,RZ(60C) ,LCUT(60C)  ,LL,LD 
CCNMCN/CCM2/RCE,RCE2,RCEM,AC,PAC,FFTC,PTC,PFPTC,FPTC, 

X   Ff,PFIV,TPCT 
CCNMCN/C0M3/EXA,EXB,FXA,PEXA,PEXe,PFXA,TFERi* 
CCNMCN/C0M5/FX(60G) , FY (600) ,FZ(600),PPE(600) 
CCMN0N/C0M2  2/R0EA,R0EB,R0EC,R0EC2,CP0,CP1,CF2,CP3  ,CFC, 

X   CF1,CF2,CGD1,CGC2,CGB1,CGE2,CGF,IPR0G 
1  =  1 

IF(LCUT (I)  .NE.O)  GO  TO  200 
IF=I+1 

INTIAL  FORCE  DISTANCE  BOUNDARIES 

XNAX  =  RX(I  )+ROE 
YN£X  =  RY (I  )  +  ROE 


114 


c 
c 

c 

c 
c 


c 
c 
c 

c 


165 


18C 


195 
2CC 


ZNAX  =  RZ(I  )+ROE 

XMN  =  RX(I)-ROE 

YNIN  =  RY(I  )-ROE 

If  IN  =  RZ(I  )-RQE 

TEST  TO  DETERMINE  IF  BULLET  AND  ATQNS  ARE 
CLOSE  ENOUGH  TO  INTERACT  AND  WHICH  FCRCE 
FUNCTION  TC  USE 


CC  195  J=IP,LL 
IF(LCUTU)  .GT.O) 
IF(RY(J  ).GE.YMAX) 
IF  (RY( J  J.LE.YMIN) 
IF(RX(J).GT.XMAX) 
IF(RX(J  ).LE.XMIN) 
IF(RZ(J).GE.ZMAX  ) 
IF(RZ(J ) .LE.ZMIN) 


GC  TO  195 

GO  TO  195 

GO  TO  195 

GO  TO  195 

GO  TO  195 

GO  TO  195 

GC  TO  195 


A  CIRCLE 


TEST  TO  DETERMINE  IF  ATCM  IS  INSIOE 
OF  RADIUS  RCE  FROM  THE  EULLET 

CRX=RX( J)-RX( I) 

CPY=RY( J)-RY( I ) 

CPZ=RZ( J  )-RZ( I) 

DIST=DRX**2+DRY**2+DRZ**2 

IF(DIST.GE.R0E2)     GC    TO    195 

CIST=SQRT(CIST) 

IF(DIST.GT.RCEM)     GC    TO    165 

ACTUAL     FORCE    FUNCTION    DETERMINED    FRCM    THE 
BCRM-MEYER    BULLET-ATCM    PCTENTIAL 

FCRCE=EXP(PFXA+PEXB*DIST) 
GC    TC    180 
CFF=ROE-DIST 
IF(DFF.LT.l.OE-lO)    GO    TO    195      f 

AN    AVERAGE    FORCE    FUNCTION    FOR    THE    BLLLET-ATCM 
INTERACTICN    AT    TURNCN 

FCPCE=(EXP( FAC+PEXB*DIST)-PFPTC)/CFF 
IF(FM.GT.FCPCE)     GC    TO    195 
FCC    =FGRCE/CIST 


FCO*DPX 
=    FX(J)+FA 
=    FX(  D-FA 

FCC*DRY 
=    FY(J)+FA 
=    FY(  D-FA 

FCD*DRZ 
=    FZ(J)+FA 
=FZ( I )-FA 


GO    TO    300 


F/S    = 

FX(J) 

FMI) 

FA    = 

FY(  J) 

FY(I) 

FA    = 

FZ(J) 

FZ(  I) 

CCNTINUE 

CCNTINUE 

CC    200    I=2tLD 

IFUCUT  (I  )  .GT.O) 

IF=I+1 

XMN    =    RX(  I  )-ROE 

XVAX    =    RX(  I  )+ROE 

YMN    =    RY  (I  )-ROE 

YNAX    =    RY(  I  )  +  ROE 

ZMN    =    RZ(I)-ROE 

ZNAX    =    RZID  +  ROE 

TEST  TO  DETERMINE  WHETHER  ATOMS  ARE  CLOSE 
ENOUGH  TO  INTERACT  AND  WHICH  FORCE  FUNCTION 
TO  USE. 

DC  295  J=IPtLL 
IFUCUT(J)  .GT.O)  GC  TO  295 
IF(RY{J) .GE.YMAX)  GC  TO  295 
IF(RY( J J.LE.YMIN)  GO  TO  295 
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c 
c 

c 
c 


28C 


295 
3CC 


IF(RX(J) .GE.XMAX)  GC  TO  295 
IF(RX(J ).LE.XMIN)  GO  TO  295 
IF(RZ(J).GE.ZMAX)  GO  TO  295 
IF(RZU).LE.ZMIN)  GO  TO  295 
CRX  =  RXU)~RX(  I  ) 
DPY  =  RY( J)-RY(  I) 
CRZ=RZ( J)-RZ(I) 
DIST=DRX**2+DRY**2+DRZ**2 
IF(DIST.GE.R0E2)  GO  TO  295 
DIST=SQRT(DIST) 
IF(DIST.GT.ROEM)  GC  TO  265 

ACTUAL  FORCE  FUNCTION  DETERMINED  FRCM  THE 
CU-CU  (GIB-IIJ  POTENTIAL  FUNCTION 


FCRCE=E>P(FXA+EX3*DIST) 
GC  TO  280 
265  CFF=RCE-DIST 

IF(CFF.LT  .1.0  E-1C)  GO  TO 


295 


AN    AVERAGE    FORCE    FUNCTION    FOR    THE    ATCM-ATCM 
INTERACTION    AT    TURN    ON 

FCRCE=(EXP( AC+EXB*CIST)-FPTC)/DFF 

IF(FM.GT. FORCE)     GC    TO    295 

FCC=FORCE/DIST 

FA=FOD*DRX 

FX(J)=FX( J)+FA 

FX(I)=FX( I) -FA 

FA=FCD*CRY 

FY(J)=FY( J)+FA 

FV(I)=FY(I )-FA 

FA=FOD*ORZ 

FZ(J)=FZ(J )+FA 

FZ(I)=FZ(I )-FA 

CCNTINUE 

CCNTINUE 

RETURN 

END 


SLEROUTINE  KPRINT 


96  1 
964 


964 

965 


CCMMCN/ 

CCMMCN/ 
CCMMCN/ 

X  PEX,RE 
CCMMCN/ 
CCMMCN/ 
CCMMCN/ 
CCNMCN/ 
CCMMCN/ 
CCMMCN/ 
CCMMCN/ 

X  IFBUL 
CCMMCN/ 

X  LSLR 
CCMMCN/ 

X  CF1,CF 
C  FCRMAT( 
C  FCPMAT( 

X  •  ENER 

X  •  X  * 
5  FCPMAT< 

X  'LATT 
C  FCRMAT( 


THIS  SUBROUTINE  PRINTS  THE  PAGE  FEACING 


C0M1/RX<60  0) , RY(600) ,RZ(6CC),LCUT(60C),LL,LD 

C0M3/EXAT  EXB,FXA,PEXA,PEXE,PFXA,TFER^ 

C0M4/IX,IY,IZ,  IXF,IYP,SCX,SCY,SCZ,IDEEP, 

2 

C0M7/R1 ,LSSfSPX,SPZfCCX,CCY,C0ZTYLAX<2J)»ILL 

C0M11/RXI  (6C0),RYI(6  00),RZI  (600) 

C0M16/BMAS,TMAS1,TMAS2,  ABUM,A8UN2,  I  PNUM 

C0M16A/BX,BY,BZ,EDISF,CVP 

C0M17/TIME,NTtDTtDTI ,TEMF,CELS,EV,NP4GE,NLINE 

C0M18/  ITEMP, TAMP, CORF 

COM2  0/IH1(20),IHS(10),IHT(6),IHB(6),  IHTARG(2) , 

L(2) 

C0M21/IHBPCT(2) , IHTPOT ( 2 ) , SUR ( 2 ) ,  Ih2  ( 10 ) , I  OR , 

C0M22/RaEA,R0EB,R0EC,R0EC2,CPO,CPl,CF2,CP3,CF0, 

2fCGCl?CG02iCGBltCGB2,CGFl,CGF2,  IPRCG 

2CX,20A4, /,4CX,10A4,///) 

'  (SIB,1)  PLANE  ',2A4,'  SURFACE 

GY  =',F7.1,'  EV,  CRYSTAL  SIZE  (' 


12,'  )  ,  IMPACT  POINT  ' 
•  TARGET  -• ,2A4,20X, ' 


PR  IMARY' 
12,  •  X  ', 


,19/) 

BULLET  -'  ,2A4,7X 


12, 


ICE  UNIT  =«  ,F7.4,» 
•  MASS1  ='  ,F7.2,2X 


ANG'  ) 

•ABUNDANCE  =  • 


F5.3,6X, 
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X   'MASS  =•  ,F7.2t9X,  'LATTICE  TEMPERATURE  = »  t A4 , 
X   •  THERMAL  AMPLITUDE  =',F7.4) 
9652  FCRMAT('  MASS2  =  '  ,  F7  .2 ,2X , ' ABUNCANCE  =',F5.3, 

X   29X, 'CORRELATION  =',F5.3,//) 
96  5  5  FCFMAT(12X,2(2A4,1X)  ,2A4, 3X ,5 HP EXA= , F9 . 5 , 2X ,5HPEXB  =  , 

X       F9.5,2X,5FPFXA=,F9.5/) 
966C    FCFMAT(12H    POTENTIAL       ,  2( 2A 4  ,  IX  )  ,  2A4 ,3X , 5FEXA    =,F9.5, 

X   2X,5HEXB  =,F9.5,2X,5HFXA  =,F9.5) 
9665  FCPMAT(30H  TARGET  POINT  ON  CRYSTAL    X  =,F7.3,«,  Y  =  ', 
X   F7.3,',  Z  =',F7.3,13H,  COS  TO  X  =,F8.5, 
X   12H,  CCS  TO  Z  =,F8.5) 
967C  FGPMAT(30H  PRIMARY  START  POINT  (LU)   X  =,F7.3,5H,  Y  =, 

X   F7.3»5H,  Z  =,F7.3/) 
9675  FCPMAT(/«  TIMESTEF  ',I4,40X,'   ELAPSED  TIME  (SEC)  =', 

X   E10.4,',  LAST  TIMESTEP  WAS  =',E10.4/) 
9725  FCFMATC  CUT-OFF  AT',F5.2,',  WHEN  R  =',F6.3,'  LU, MORSE 
X   POTENTIAL  PARAMETERS  ARE '  ,  L0A4  ,  // ,10X ,  •    CGD1  =  ', 
X   F8.4,',  CGD2  =',F8.4,',  CGB1  =',F8.4,»t  CC-B2  =  ', 
X   F8.4,',  CGF1  =«,F8.4,»,  CGF2  =',F8.4,/) 
973C  FCFMAT('  WHEN »  , F8 .4 , '  <R<',F8.4,«  THE  MATCHING  POTENTI 
XAL  PARAMETER  ARE',//,3X,'      CPC  =',E12.4,',  CP1  =  », 
X   E12.4,'8  CP2  =',E12.4,»,  CP3  =  '  ,E12 .4 , /7X  ,  •  CFO  =  ', 
X   E12.4,',  CF1  =',E12.4,',  CF2  =',E12.4,//) 
974C  FCPMATQlh  THE  FIRST  ,I3,22F  ATOMS  ARE  BLCCKEDt    ,13, 
X   24H  LAYERS  ARE  FREE  TC  MOVE,/) 
WRITE  (  6,9610)  Ihl,IHS 

WRITE  (  6,9640)  I CR , SUR,EV, IX , I Y , IZ, IPNUM 
WRITE  (  6,9645)  IHTARG , IHBULL , CVR 
WRITE  (6,9650)  TMAS  1 , ABUN 1 , BMAS , I  TEMP , TAMP 
WRITE  (6,9652)  TM AS2 , ABLN2 , COR F 
WRITE  (  6,9660)  I HT t E XA  ,  EXB , FXA 
WRITE  (  6,9655)  I HB , PEXA, PEXB , P FXA 
WRITE  (  6,9665)  BX , BY , B Z , CO  X ,CCY  ,CCZ 
WRITE  (  6,9670)  RX  I  ( 1  )  , RYI ( 1 )  , RZ  I  ( 1 ) 
IFUPRCG.NE.4)  GO  TO  100 

WRITE  (  6,9730)  RCEA , ROEB , C FO , CP 1  ,CP2 , C P3  ,  C FO  ,CF1 , CF2 
WRITE(6,972  5)R0EC,RCEB,IH2,CGD1,CGC2,CG31,CGB2,CGF1, 
X   C  G  F  2 
IOC  WRITE  (  6,9675)  NT, TIME, DT 
NLINE=1 
RETLRN 
END 


SLBROLTINE  KPNT 
C 
C 

C  THIS  SUBRCUTINE  PRINTS  THE  HEADING 

C  FOR  OUTPUT  THAT  IS  FILED 


CCMMCN/C0M4/IX,IY,IZ, IXP , IYP  ,3CX  ,  SCY ,SC Z , I CEEP , 
X  R  P  X  R  R  7 
CCMMCN/C0M16/BMAS  ,TMAS1  ,TMAS2 , A BUM , ABUN2 , I FNUM 
CCNMCN/C0M16A/BX,BY,8Z,EDISP,CVR 

CCMMCN/C0M17/TIME,NT,DT,DTI ,TEMP,CELS,EV ,NP AGE ,NL INE 
CC^MCN/C0M20/IH1(20),IHS(10) , I HT ( 6 ) , IHB ( 6 ) ,  IHTARG ( 2  ) 
X   IHBULL(2) 

CCNPCN/C0M21/IHBPCT(2) , IHTPOT (2 ) , SUR(2) , IH2 ( 10 ) , I  OR, 
X  LSUR 
964C  FCRMAT('  (',13,')  PLANE ', 10X ,'  ION  ENERGY  =',F7.L,10X 

X   'CRYSTAL  SIZE  (',12,'  X',I2,'  X ',  12 ,')',// ) 
9645  FCPMATC  TARGET  -',2A4,20X,'  BULLET  -',2A4,6X, 

X  'PCTENTIAL   '  ,2 ( 2A4  ,  IX ) , 2A4) 
965C  FCRMAT(  »   MASS  1  ='  ,F7 . 2 , 2X  ,' ABUNCANCE  =',F5.3,4X, 

X   'MASS  =  '  ,F7.2,16X,2(2A4,1X)  ,2A4) 
9652  FCPMAT  (•   MASS  2  =', F7 .2 , 2X ,' ABUNDANCE  =',F5.3,//) 
WRITE  (6,9640)  I  OR  ,E V , I  X , I Y  ,  I Z 
WRITE  (6,9645)  I HTARG, I HBUL L , I HT 
WRITE  (6,9650)  TMAS1 , ABUN1 , EMAS , IHB 


117 


WRITE  (6,9652) 

RETURN 

ENC 


TMAS2,ABUN2 


C 

c 
c 

c 

c 
c 


SUBROUTINE  HISTOl IN , FREG1 , RANGE  ,  T  ITLE1 ) 

THIS  SUBROUTINE  CRAWS  THE  FISTORAM  (FREQUENCY 
VS.  EJECTEC  ANGLE)  FOR  THE  TYPE  1  HEAVY  ATCMS 

DIMENSICN  FREQK  20),  RANGE  (20),  JOLT  (20)  ,  TITLE  1(24) 
CINENSICN  PCUT(20) 
CATA  NQTH/'     »/,K/»****»/ 
3  FCPMAT( •0INTERVAL,»2X,lS(F5.1t IX) ,F5.1) 


10 
12 


7C 

£0 


16 
51 

15 

ICC 


FCFNAT(//,8X,24A4/) 
FCRNAT(  'OFREQUENC V 
FCPMAT( 'OINTERVAL • , 
-«  ,128( »-' 
TITLE1 


FCFNAT( 

WRITE  (6,4) 

JN=IN 

IM  =  0 

IF(JN.GT.O)  GO  TO  10 

INT=1 

J!^=-JN 

CC  12  1=1, JN 

JCLT(I) =FREC1( I ) 

WRITE  (6,5)  (JOUT(I), 

WRITE  (6,7) 


,2016) 

2X,19(F5.1,1X)  ,F5.1) 
)) 


1=1, JN) 


FIND  LARGEST  FREQUENCY 


FNAX=0. 
CC  20  1=1, JN 
IF(FREQ1 (I ) .GT 
2C  CCNTINUE 

SCALE 


FMAX)     FMAX=FPEQ1(I ) 


JSCAL=1 

I  F(FMAX.GT.6J.)  J SCAL= ( FMAX+59 . )/60. 

FSCAL=JSCAL 

CC  50  1=1, JN 

JCLT( I) =NOTH 

yAX=FMAX/FSCAL 

CC  8J  1=1, MAX 

X=MAX-( 1-1) 

CC  70  J=1,JN 

I  F(  FREQK  JJ/FSCAL.GE.X  )    JOUT(  J)=K 

CCNTINUE 

IX=X*FSCAL 

WRITE  (6,2)  IX, ( JCUT( J)  ,J  =  1  ,JN) 

FCFMAT( I6,4X,20(2X,A4) ) 

WRITE  (6,7) 

IFdNT.EQ.l  )  GO  TO  16 

WRITE  (6,3)  (RANGEU  ),  J  =  l,  JN) 

GC  TO  15 

CC  51  1=1, JN 

FCLT(I)=RANGE( I ) 

WRITE  ( 6,6)  (POUT  (I  )  ,1=1, JN) 

CCNTINUE 

WRITE  (6,130) 

FCPMAT(  IX, 'MIDPOINT1/) 

RETURN 

ENC 
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c 
c 
c 
c 
c 
c 


1C 
12 


5C 


7C 

80 


16 

51 

15 

1GQ 


SLEPGUTINE  HIS  TO ( IN , FREC2 ,R ANGE , T  ITLE2 ) 


THIS  SUBROUTINE  CRAWS  THE  HISTORAM  (FREQUENCY 
VS.  EJECTEC  ANGLE)  FGR  THE  TYPE2  LIGHT  ATOMS 


CIMENSICN  FREQ2(2J),RANGE<20), J OUT (20) 

CINENSICN  PCUT(20) 

CATA  NGTH/'     •/»!</•****■/ 

FCFNAT(  I6,4X,20(2X,A<t)  ) 


TITLE2(24) 


FCPMAT( 'OINTERVAL- , 
FCFMAT(//,8X,24A4/) 
FCFNAT(  'OFREQUENCY' ,2016) 
FCPMAT( 'OINTERVAL •  ,2X,19(F5 
FCFMAT( '  -•  ,128( •-•  )) 
WRITE  (6,4)  TITLE2 
J  N  =  I N 
INT  =  J 

if(jn.gt.o)  go  to  10 
int=i 

ji^=-jn 

DC  12  1=1, JN 

JCLT(I) =FREC2(  I) 

WRITE  (6,5)  (  JOUTd  ),I  =  1,JN) 

WRITE  (6,7) 

FINC  LARGEST  FREQLENCY 


2X,19(F5.1,1X)  ,F5.1) 


1,1X)  ,F5.1) 


FNAX=0. 
CC  20  1=1, JN 
IF(FREQ2(I ) .GT.FMAX) 
2C  CONTINUE 


FMAX=FREQ2(I) 


SCALE 


)  JSCAL=(FMAX+59. )/60 


JSCAL=1 

IF(FMAX.GT.60 

FSCAL=JSCAL 

DC  50  1=1, JN 

JCLT(I) =NOTH 

MAX=FMAX/FSCAL 

CC  80  1=1, MAX 

X  =  PAX-(  1-1) 

CC  70  J=1,JN 

IF(FREQ2(J)/FSCAL. GE.X  ) 

CCNTINUE 

IX=X*FSCAL 

WRITE  (6,2) 

WRITE  (6,7) 

IF( INT.EQ.l  ) 

WPITE  (6,3) 

GC  TO  15 

CC  51  1=1, JN 

PCLTd  )=RANGE(  I  ) 

WRITE  (6,6)  (P0UT(I),I=1,JN) 

CCNTINUE 

WRITE  (6, ICC) 

FCPMAT( IX, 'NIDPOINT' 

RETURN 

END 


JOUT( J)=K 


IX,  ( JCUT( J) ,J  =  1 ,JN) 

GO  TO  16 
(RANGE( J  ),J  =  1,JN) 


/) 


C 

c 
c 
c 
c 
c 


SLEPOUTINE  KFOLDA 


THIS  SUBROUTINE  STATISTICALLY  FOLDS  THE 
TARGET  IMPACT  AREA  FGR  THE  FCC  (100)  CRYSTAL 
LATTICE  TO  REPRESENT  FULL  COVERAGE 
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CCMMON/CCM1/RX(60  0)  ,RY(600)  ,RZ ( 600 ) ,LCUT ( 60 Z ) ,LL,LD 
CCMM0N/C0M4/IX,I Y,IZ,  IXPTIYP,SCX,SCY,SCZ,IDEEP, 
v  c  P  x  R  B  Z 
COMCN/C0M4A/  LNUM6J0)  ,LNUMI  (  6  0  0  )  ,LNUMF  (  6  0  C  ) 
CC^M0N/C0M21/IHBPCT(2) , IHTPCT ( 2 )  ,  SUR ( 2 )  ,  I  h2  ( 10 )  ,  I  OR 
X  LSLR 
ILAYER=1 
NFIRST=2 
C 

C  DETERMINE  THE  FIRST  ATOM  AND  THE  LAFGEST 

C  NUMBERED  ATOM  IN  EACH  L£YER 

C 

10  GC  TO  (20,3C,40,5C,6C,70) ,  ILAYER 
20  MAY  =(  (LL-l)/(  IY/2)  )/2+2 
MNFIRST 
NLAY1=NLAY 
GC  TO  8C 
20  NLAY=((LL-l)/(IY/2))+l 

K=NFIRST  +  IX/2 
rr     tq  or 

40  MAY  =  NLAY1  +  NLAY  -1 
K=NFIRST 
GC  TO  80 
5C  NLAY  =  NLAY1  +  NLAY  -2 
K=NFIRST  +  IX/2 
GC  TO  80 
6C  NLAY  =  NLAY1  +  NLAY  -3 
K=NFIRST 
GC  TO  80 
7C  NLAY  =  NLAY1  +  NLAY  -4 
K=fSFIRST  +  IX/2 
GC  TO  8C 
C 

C  FIRST  (DIAGCNAL)  FOLD  -  ALL  LAYERS 

C 

80  DC  90  J=NFIRST,NLAY 

LNUM  J)  =LNUN(  J)  +    LNUM(K) 
IF(J.EQ.K)  LNUM(J)=2*LNUM(J  ) 
LNLMK)=LNUV(J ) 
K  =  K  +  I  X 

IF(K.GT  .NLAY)  K=K-NLAYl+2 
9C  CCNTINUE 

DC     100    J=NFIRST,NLAY 
LNUM  J)=LNUM(J ) /2 
IOC    CCNTINUE 
C 

C  SECCND    (HORIZONTAL)     FOLC    -    ODD    LAYERS 

r 

105    ICCUNT=1 

ISET=0 

GC    TO    {  110,160,  110,160,  110,  160)  ,     ILAYER 
110    K=NLAY-IX/2 

CC    150    J=NFIRSTtNLAY 

LNLM  J)=LNUM<  J)+LNUM(K) 

IF(J.EQ.K)     LNUM< J)=2#LNUM(J ) 

LNLM(K)=LNUN(J) 

ICCUNT=ICOLNT+l 

IF(ISET.EQ.l)    GO    TO    130 

IF(  ICCUNT.GT.    IX/2    +1)     GO    TO    120 

K=K  +  1 

GC  TO  150 
12C  ISET=1 

K=K-IX+1 

ICCUNT=1 

GC  TO  150 
13C  IF  (ICCUNT.GT,  IX/2)  GO  TO  140 

K=K  +  1 

GO  TO  150 
L4C  ISET=0 

K-K-IX+1 

ICCLNT=1 
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c 
c 
c 


170 


18J 


190 


2JC 
210 

22C 


15C  CCNTINUE 
GC  TC  210 

SECCNO  (HCRIZQNTAL)  FOLC  -  EVEN  LAYERS 

16C  K=NLAY-IX/2+l 

DO  200  J=NFIRST,NLAY 

LNUM  J)=LNUM  J)+LNUM(K) 

IF(J.EQ.K)  LNUM(J)=2*LNUM(J) 

LNLM(K)=LNUN( J) 

ICCUNT=ICGUNT+1 

IF(ISET.EQ.l)  GO  TO  180 

IF(  ICGUNT.GT.IX/2)  GO  TO  170 

K  =  K  +  1 

GC  TO  2C0 

ISET=1 

K=K-IX+1 

ICCLNT=1 

GC  TC  2G0 

IF  ( ICCLNT.GT.IX/2+1)  GO  TO  190 

K=K  +  1 

GC  TO  200 

I<ET=0 

K=K-IX+1 

ICCUNT=1 

CCNTINUE 

DC    220    J=NFIRST,NLAY 

LNLM  J)=LNUM(  J)/2 

CCNTINUE 

THIRD    (VERTICAL)     FCLD    -    CDC    LAYERS 

ICCUNT=1    « 

ISET=0 

GC    TO    ( 23C,280,23C,28C,230,280) ,     ILAYER 
23C    K=NFIRST+IX/2 

DC    270     J=NFIRST,NLAY 

LNLM J)=LNUM( J)+LNUM(K) 

IF(J.EQ.K)     LNUM( J  )=2*LNUM(J  ) 

LMjM(K)  =LNUM(J) 

ICCUNT=ICCUNT+1 

IF(  ISET.EQ.l)  GO  TC  250 

IF(ICCUNT.GT.IX/2  +1)  GO  TO  240 

K=K-1 

GC  TO  270 
24C  ISET=1 

K=K+IX-1 

ICCLNT=1 

GC  TO  270 
250  IF  (ICCUNT.GT.IX/2)  GO  TO  260 

K=K-1 

GC  TO  270 
26C  ISET=0 

K=K+IX-1 

ICCUNT=1 
27C  CCNTINUE 

GC  TO  330 

THIRD  (VERTICAL)  FQLC  -  EVEN  LAYERS 

280  K=NFIRST+IX/2-l 

DC  220  J=NFIRSTTNLAY 

LNUN ( J)=LNUM(J)+LNUM(K) 

IF(J.EQ.K)  LNUM( J)=2*LNUM(J) 

LNUMi K)=LNUN(J) 

ICCLNT=ICCUNT+1 

IF(  ISET.EQ.l)  GO  TC  300 

IF(  ICOUNT.GT.IX/2)    GO    TC    290 

K=K-1 

GC    TO    2  20 
290     ISET=1 

K=K+IX-1 
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3QC 


31C 


32C 
33C 

34C 


35C 


ICCLNT=1 

GC  TO  320 

IF<  ICOUNT.GT.IX/2+1)  GG  TO  310 

K=K-1 

GC  TO  320 

ISET=0 

K=K+IX-1 

ICCUNT=1 

CCNTINUE 

CC  340  J=NFIRST,NLAY 

LNLM(  J)=LNUM(J)/2 

CCNTINUE 

ILAYER=ILAYER+1 

NFIRST  =  N'LAY+1 

IF(  ILAYER.GT.IY)  GC  TO  350 

GC  TO  1C 

RETURN 

ENC 


C 

c 
c 
c 
c 

c 

c 


SIEROUTINE  KF0LD3 


THIS  SUBROUTINE  STATISTICALLY  FOLCS  THE 
TARGET  IMPACT  AREA  FOR  THE  FCC  (110)  CRYSTAL 
LATTICE  TC  REPRESENT  FULL  COVERAGE 

CC^MCN/CCM1/RX(60C)  ,RY(600),RZ(600),LCUT(6QC),LL,LC 
CCMMCN/C0M4/IX, IY,IZtIXP,IYP,SCXtSCY,SCZ,IDEEPt 

X  REX  R  B  Z 
CCfMCN/C0M4A/  LNUM600),LNUMI(600)  ,LNUMF(60C) 
CC^MCN/C0M21/IHBP0T(2) t  IHTPQT(2)  ,  SUR(2)  ,IF2(10),I0R, 

X  LSLR 
ILAYER=1 
NFIRST=2 


DETERMINE  THE 
NUMBERED  ATOM 


FIRST  ATOM  AND  THE  LARGEST 
IN  EACH  LAYER 
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NUMBERED    ATOM    IN    EACH    LAYER    FOR     110    LATTICE 
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NUMBERED    ATGM    IN    EACH    LAYER    FOR    111    LATTICE 

2CC    GC    TO    (210, 220, 230, 240), ILAYER 
21C    M_AY=((LL-2)/( IY/2) )/2+2 

MAY1  =  NLAY 

GC    TO    8C 
22C    NLAY=( ( LL-2)/( IY/2) )+l 

GC    TO    8  C 
22C    MAY  =  M_AY1  +  NLAY-1 

C-  c    t  r    8  c 
24C    MAY=MAY1  +  NLAY-1 
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Figure    1    -      Sputtering    Summary   -     (100)     2keV    Cu/Ne+    KSE-B 


Figure    2    -      Sputtering    Summary    -     (100)     2keV    Cu/Au+    KSS-5 
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Figure  3  -   9x4x9,  (100)  Orientation  Microcrystallite 
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Figure  U   -   9x4x9,  (110)  Orientation  Hicrocrystallite 
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Figure  7  -   intrinsic  Geometry  of  Target  Impact  Areas 
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Figure  8  -   Reduced  Representative  Impact  Areas 
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Figure    12    -       (100)     and    (110)    Orientation    Folding    Geometry 
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Figure  13  -   (111)  Orientation  Folding  Geometry 
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Fiaure    14   - 


Sputtering    Ratio    vs.     Ion    Atomic    Number 
(100)  ,     100-eV,     9x4x9 
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(a)Gibson-II 


(b)KSE-B 


Figure    15    -      Sputtering    Summary    -     (100) ,    100-97,    Cu/ff=v 
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Figure    16    -      Ring   Probability    -     (100)  ,     100-eV,    Cu/ i 
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Figure    17   -      Angular   Dispersion   -     (100),    100-aV,    Cu/He+ 
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Fiaure    18   -      Sputtering    Summary   -     (100),    100--V,    Cu/Na+ 
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Figare    19    -      Ring   Probability   -     (100)  ,    100-eV,    Cu/Na 
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Figure    20    -      Angular    Dispersion    -     (100) ,    100-eV,    Ca/Ne+ 
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Figure    21    -      Sputtering    Summary    -     (100),     100-eV,    Cu/Ar+ 
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Figure  22  -   Ring  Probability  -  (100)  ,  100-eV,  Cu/Ar- 
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Figure    23   -      Angular    Dispersion    -     (100)  ,    100-eV,    Cu/Ar  + 
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Figure    24    -      Angular    Dispersion    -     (100)  ,    100-eV,    Cu/Ar  + 
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Figure    25    -      Sputtering    Mechanism    Trac9 
(100),    100-eV,    Cu/Ar+,    Gibson-IT 
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.gure    26    -      Sputtering    Summary    -     (100),     100-eV,    Cu/Cu 


153 


100  + 


0.75-  • 


—   0.50- 

CO 


(a)Gibson-II 


«-H 1- 


co 

o 

< 


CO 

cc 

CO 


a. 

CO 


1.00  + 


0.75-- 


0.50  •• 


0.25-- 


0.00 


(b)KSE-B 

2                  3                 4  5  6 

H 1 1  I  I 


RING     DISTANCE     FROM     REFERENCE     POINT     [LU] 


Figure    27    -       Ring    Probability    -     (100),     100-<?V,    Cu/Cir 
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Figure    28    -      Angular    Dispersion    -     (100),     100-97,    Cu/Cu+ 
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Figure    29    -      Sputtering    Summary    -     (100)  ,    100-eV,    Ci/*-  + 
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Figure    30    -       Ring    Probability    -     (100)  ,     100-eV,    Cu/Kr  + 
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Figure    31    -      Angular    Dispersion    -     (100),     100-eV,    Cu/Kr+ 
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Figure    32    -      Sputtering    Summary    -     (100),    100-eV,    Cu/X 
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Figure    33    -       Ring    Probability    -     (100),     100-e7,    Cu/X 
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Figure    34    -      Angular    Dispersion    -     (100),    100-eV,    Cu/Xe+ 
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Figurs    35    -      Sputtering    Summary   -     (100)  ,    100-eV,    Cu/Au  + 
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Figure    31    -      Angular    Dispersion   -     (100),    100-eV,    Cu/At 
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Figure    38    -      Sputtering    Summary    -     (100),     100-eV,    Cu/Hg 
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Figure    40    -      Angular    Dispersion    -     (100) ,    100-sV,    Cu/Ha* 
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Figure    41    -      Angular    Dispersion    -     (100),    100-27,    Cu/Kg+ 
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Figure  43  -   Spattering  Summary  -  (100),  1-kaV,  Cu/H« 
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Figure    44    -      Ring   Probability   -     (100) ,    1-keV,    Cu/He+ 
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Figure    45   -      Angular    Dispersion   -     (100)  ,    1-ksV,    Cu/He"1- 
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Figure  46  -   Sputtering  Summary  -  (100),  1-ksV,  Cu/Ne+ 
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Figure  47  -   Ring  Probability  -  (100) ,  1-keV,  Cu/N9- 
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Figure    48   -       Angular    Dispersion    -     (100)  ,    1-kaV,    Cu/Ne-  + 
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Figure  49  -   Spattering  Summary  -  (100),  1-keV,  Cu/Ar+ 
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Figure    50    -      Ring   Probability    -     (1u0) ,    1-kaV,    Cu/Aj 
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Figure    51    -       Angular    Dispersion    -     (100) ,    1-k=V,    Cu/Az+ 
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Figure  52  -   Spattering  Summary  -  (100),  1-ksV,  Cu/Cu+ 
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Figure    53    -      Ring   Probability    -     (100)  ,    1-keV,    Cu/Cu" 
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Figure    54    -       Angular    Dispersion    -     (100),    1-keV,    Cu/Cu+ 
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Figure  55  -   Sputtering  Summary  -  (100),  1-JcaV,  Cu/Kr+ 
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Figure    56    -      Ring    Probability    -     (100) ,    1-keV,    Cu/Xr+ 
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Figure    57   -       Angular    Dispersion    -     (100),    1-kaV,    Ci/Kr+ 
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Figure  58  -   Sputtering  Summary  -  (100)  ,  1-ksV,  Cu/X9  + 


185 


0.4+ 


0.3-- 


, , 

0.2 

co 

t- 

3 

a. 

CO 

— 1 

0.1 

< 

H- 

O 

H- 

0.0 

(a)Gibson-II 


CO 

5 

o 

< 


CD 


CO 


Ou 
CO 


(b)KSE-B 

0                 12                 3  4  5  6 

I 1  t III* 

RING     DISTANCE     FROM     REFERENCE     POINT     [LU] 


Figure    59    -      Ring   Probability    -     (100)  ,    1-ke7,    Cu/X^ 
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Figure    60   -      Angular   Dispersion   -     (100),    1-keV,    Ci/Xe*- 
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Figure    61    -       Sputtering    Summary    -     (100),    1-ksV,    Cu/Au+ 
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Figure    62    -      Ring   Probasilit y    -     (100),    1-keV,    Cu/Au+ 
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Figure    63  -      Angular   Dispersion    -     (100),    1-kaV,    Cu/Au+ 
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Figure  64  -   Sputtering  Summary  -  (100),  1-kaV,  Cu/Hj+ 
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Figure    65    -      Ring    Probability    -     (100),    1-keV,    Cu/Hcr 
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Figure    66    -       Angular    Dispersion    -     (100),    1-keV,    Cu/Ha+ 
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Figure  96  -   Angular  Dispersion  -  (100),  100-97,  Cu/C'i  + 
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Figure  105  -   Angular  Dispersion  -  (100),  1-keV,  Cu/Ci 
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Figure  107  -   Atom  Probability  -  (100),  100-aV,  Cu/Hg 
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Figure  108  -   Angular  Dispersion  -  (100),  100-eV,  Cu/"jj  + 
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Figure    110    -      Atom    Probability    -     (100),    200-9V,    Cu/Hg+ 
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Figure  111  -   Angular  Dispersion  -  (100),  200-eV,  Cu/Hg+ 
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Figure    113    -      Atom    Probability    -     (100),    500-sV,    Cu/Hg+ 
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Figure    114   -       Angular    Dispersion    -     (100),    500-eVf    Cu/H< 
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Figare  116  -   Atom  ProDability  -  (100),  1-keV,  Cu/Hg 
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Figure  118  -   Sputtering  Mechanism  Trace 
(100),  100-eV,  Cu/Cu+,  KSE-B 
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'igure  119  -   Sputtering  Mechanism  Tr^ce 
(100),  100-eV,  Cu/Cu+,  KSE-3 
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Figure    122    -      Placement    of    144    Random    Shots 
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Figure  123  -   Sputtering  Summary  -  Random  Points 
(100),  100-eV,  Cu/Hg+,  KSE-3,  11xUx11 
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Figure  124  -   Sputtering  Summary  -  Grid  Points 
(100),  100-eV,  Cu/Hg+,  KSE-3,  11x4x11 
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Figure    125    -      Atom    Spattering    Probability    -    Random    Points 
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Figure    126    -       Atom    Sputtering    Probaoility   -    ^rid    Poinds 
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figure  127  -   Angular  Dispersion  -  Random  Points 
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Figure    128    -      Angular   Dispersion    -    Grid    Point: 
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